diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/Makefile agcm5.3-g77/src/Makefile *** agcm5.3/src/Makefile Thu Mar 5 20:00:00 1998 --- agcm5.3-g77/src/Makefile Tue Feb 23 22:42:18 1999 *************** *** 146,151 **** --- 146,153 ---- cd .. ; \ done + TAGS=ctags + tags: tags.local tags.dir tagssm: tagssm.local tagssm.dir diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/Mkinclude agcm5.3-g77/src/Mkinclude *** agcm5.3/src/Mkinclude Thu Mar 5 22:26:23 1998 --- agcm5.3-g77/src/Mkinclude Mon Jul 5 18:10:32 1999 *************** *** 8,21 **** ############## 分解能の指定 ######################## ###RESOLUTION = t10l01 ###RESOLUTION = t10l02 ###RESOLUTION = t10l12 ###RESOLUTION = t10l45 ###RESOLUTION = t10l12xz ###RESOLUTION = t10l45xz ! ###RESOLUTION = t10l5 ! RESOLUTION = t10l5sm ###RESOLUTION = t21l02 ###RESOLUTION = t21l12 --- 8,23 ---- ############## 分解能の指定 ######################## + RESOLUTION = t42l16 + ###RESOLUTION = t10l01 ###RESOLUTION = t10l02 ###RESOLUTION = t10l12 ###RESOLUTION = t10l45 ###RESOLUTION = t10l12xz ###RESOLUTION = t10l45xz ! ## RESOLUTION = t10l5 ! ###RESOLUTION = t10l5sm ###RESOLUTION = t21l02 ###RESOLUTION = t21l12 *************** *** 24,30 **** ###RESOLUTION = t42l01 ###RESOLUTION = t42l02 ###RESOLUTION = t42l12 ! ###RESOLUTION = t42l16 ###RESOLUTION = t42l16yz ###RESOLUTION = t42l23 --- 26,32 ---- ###RESOLUTION = t42l01 ###RESOLUTION = t42l02 ###RESOLUTION = t42l12 ! ## RESOLUTION = t42l16 ###RESOLUTION = t42l16yz ###RESOLUTION = t42l23 *************** *** 38,44 **** ############# 実行ファイル名指定 agcm5.$(RESOLUTUON)$(MODELNAME) ##### ! MODELNAME = test # 適当な名前をつけましょう. ############## usr/ (利用者設定ルーチン) の指定 ########### --- 40,46 ---- ############# 実行ファイル名指定 agcm5.$(RESOLUTUON)$(MODELNAME) ##### ! MODELNAME = adj # 適当な名前をつけましょう. ############## usr/ (利用者設定ルーチン) の指定 ########### *************** *** 48,81 **** ############## nonstd/ (非標準ルーチン) の指定 ########### ###### 積雲パラメタリゼーション ##### ! CUMLUS = # 標準 Kuo 方式 ! ###CUMLUS = nonstd/p2cuma.o # 対流調節方式 ! ###CUMLUS = nonstd/p2cum0.o # 積雲パラメタリゼーションなし ###### 乾燥対流調節 ###### DRYADJ = # 標準 あり ! ###DRYADJ = nonstd/p2dad0.o # 乾燥対流調節なし ###### 浅い積雲対流 ###### SHLCNV = # 標準 Tiedtke ! ###SHLCNV = nonstd/p2shl0.o # 浅い積雲なし ##### FFT ####### FFT = # 標準 Numerical Receipes, 2^n のみ ! ###FFT = nonstd/ufftp.o # NCAR/Nakamura, 汎用 ##### 放射 ####### RADIAT = # 標準 ! ###RADIAT = nonstd/p2rad0.o # 放射なし ##### 鉛直拡散 ####### VDFUS = # 標準 ###VDFUS = nonstd/p2vdf0.o # 拡散係数一定(バルク法はそのまま) ! ###VDFUS = nonstd/p2vdf0.o \ ! ### nonstd/p2sfc0.o # 拡散係数一定 ##### 重力波抵抗 ####### ! ###GRVD = # 標準 GRVD = nonstd/p2gra0.o # 重力波抵抗なし NONSTD = $(CUMLUS) $(DRYADJ) $(SHLCNV) $(FFT) \ --- 50,83 ---- ############## nonstd/ (非標準ルーチン) の指定 ########### ###### 積雲パラメタリゼーション ##### ! ### CUMLUS = # 標準 Kuo 方式 ! CUMLUS = nonstd/p2cuma.o # 対流調節方式 ! #CUMLUS = nonstd/p2cum0.o # 積雲パラメタリゼーションなし ###### 乾燥対流調節 ###### DRYADJ = # 標準 あり ! #DRYADJ = nonstd/p2dad0.o # 乾燥対流調節なし ###### 浅い積雲対流 ###### SHLCNV = # 標準 Tiedtke ! #SHLCNV = nonstd/p2shl0.o # 浅い積雲なし ##### FFT ####### FFT = # 標準 Numerical Receipes, 2^n のみ ! #FFT = nonstd/ufftp.o # NCAR/Nakamura, 汎用 ##### 放射 ####### RADIAT = # 標準 ! # RADIAT = nonstd/p2rad0.o # 放射なし ##### 鉛直拡散 ####### VDFUS = # 標準 ###VDFUS = nonstd/p2vdf0.o # 拡散係数一定(バルク法はそのまま) ! #VDFUS = nonstd/p2vdf0.o \ ! nonstd/p2sfc0.o # 拡散係数一定 ##### 重力波抵抗 ####### ! # GRVD = # 標準 GRVD = nonstd/p2gra0.o # 重力波抵抗なし NONSTD = $(CUMLUS) $(DRYADJ) $(SHLCNV) $(FFT) \ *************** *** 92,116 **** ############## システム指定(SYSMAKE) ################### ! SYSTEM = Sun ###SYSTEM = HITAC ###SYSTEM = SX3 ############## コンパイルオプションの指定 ########## ! STDFFLAGS = $(SYSFFLAGS) -O ###STDFFLAGS = $(SYSFFLAGS) -g -u ! STDCFLAGS = $(SYSCFLAGS) -O ###STDCFLAGS = $(SYSCFLAGS) -g ! STDLDFLAGS = $(SYSLDFLAGS) ###STDLDFLAGS = $(SYSLDFLAGS) -g ! STDCPPFLAGS = $(SYSCPPFLAGS) ! AUTODBL = $(SYSAUTODBL) AUTODBLDEF = -DAUTO_DBL ############## AGCM のソースツリー原点 ############ ! DENNOUDIR = $(HOME)/dennou-com ! AGCMDIR = $(HOME)/dennou-com/agcm5 ############## インストールディレクトリの指定 ########## --- 94,123 ---- ############## システム指定(SYSMAKE) ################### ! SYSTEM = Linux ###SYSTEM = HITAC ###SYSTEM = SX3 ############## コンパイルオプションの指定 ########## ! # STDFFLAGS = $(SYSFFLAGS) -O ###STDFFLAGS = $(SYSFFLAGS) -g -u ! # STDCFLAGS = $(SYSCFLAGS) -O ###STDCFLAGS = $(SYSCFLAGS) -g ! # STDLDFLAGS = $(SYSLDFLAGS) ###STDLDFLAGS = $(SYSLDFLAGS) -g ! STDFFLAGS = -g -u ! STDCFLAGS = -g ! STDLDFLAGS = -g ! STDCPPFLAGS = $(SYSCPPFLAGS) -DREAL='ReaL*8' -DReal=DBLE ! #AUTODBL = $(SYSAUTODBL) ! #AUTODBLDEF = -DAUTO_DBL AUTODBLDEF = -DAUTO_DBL + AUTODBL = ############## AGCM のソースツリー原点 ############ ! DENNOUDIR = $(HOME)/e ! AGCMDIR = $(HOME)/e/agcm5 ############## インストールディレクトリの指定 ########## *************** *** 126,139 **** ############## 他のライブラリのディレクトリの指定 ###### ######## dcl ライブラリのありか ! DCLLIBDIR = /usr/local/lib/dcl-5.0 ######## symake ライブラリのありか ! SYSMAKEDIR = $(DENNOUDIR)/sysmake ######## gtool のライブラリのありか ! GTOOLLIBDIR = $(DENNOUDIR)/gt3-dcl5/lib ! GTOOLINCDIR = $(DENNOUDIR)/gt3-dcl5/include ######## ispack smpack ライブラリのありか ISPACKLIBDIR = $(DENNOUDIR)/lib --- 133,147 ---- ############## 他のライブラリのディレクトリの指定 ###### ######## dcl ライブラリのありか ! DCLLIBDIR = /usr/lib/dcl-5.0.1 ######## symake ライブラリのありか ! SYSMAKEDIR = /usr/lib/sysmake ######## gtool のライブラリのありか ! GTOOLLIBDIR = /home/toyoda/e/gt3-dcl5/lib-g ! #GTOOLLIBDIR = /usr/lib/gt3-dcl5 ! GTOOLINCDIR = /usr/include/gt3-dcl5 ######## ispack smpack ライブラリのありか ISPACKLIBDIR = $(DENNOUDIR)/lib diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/Mkinclude.rotty agcm5.3-g77/src/Mkinclude.rotty *** agcm5.3/src/Mkinclude.rotty Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/Mkinclude.rotty Mon Jul 5 18:08:06 1999 *************** *** 0 **** --- 1,208 ---- + # include file for Makefile for $(DENNOUDIR)/model/agcm/agcm5 + # + # You need dcl, sysmake, gtool libraries to install agcm5. + # + ####################################################################### + ######################## 実験設定部分ここから ######################### + ####################################################################### + + ############## 分解能の指定 ######################## + + RESOLUTION = t10l5 + + ###RESOLUTION = t10l01 + ###RESOLUTION = t10l02 + ###RESOLUTION = t10l12 + ###RESOLUTION = t10l45 + ###RESOLUTION = t10l12xz + ###RESOLUTION = t10l45xz + ## RESOLUTION = t10l5 + ###RESOLUTION = t10l5sm + + ###RESOLUTION = t21l02 + ###RESOLUTION = t21l12 + ###RESOLUTION = t21l12yz + + ###RESOLUTION = t42l01 + ###RESOLUTION = t42l02 + ###RESOLUTION = t42l12 + ## RESOLUTION = t42l16 + ###RESOLUTION = t42l16yz + ###RESOLUTION = t42l23 + + ###RESOLUTION = t84l01 + ###RESOLUTION = t84l02 + ###RESOLUTION = t84l16 + + ###RESOLUTION = v00l12 + ###RESOLUTION = v00l16 + ###RESOLUTION = v00l45 + + ############# 実行ファイル名指定 agcm5.$(RESOLUTUON)$(MODELNAME) ##### + + MODELNAME = rotty # 適当な名前をつけましょう. + + ############## usr/ (利用者設定ルーチン) の指定 ########### + + ###SPECIALDIR = usr/pochaka # ルーチンのおき場所の例です. + ###SPECIAL = usr/pochaka/p2tokubetu1.o usr/pochaka/p2tokubetu2.o + + ############## nonstd/ (非標準ルーチン) の指定 ########### + + ###### 積雲パラメタリゼーション ##### + CUMLUS = # 標準 Kuo 方式 + #CUMLUS = nonstd/p2cuma.o # 対流調節方式 + #CUMLUS = nonstd/p2cum0.o # 積雲パラメタリゼーションなし + + ###### 乾燥対流調節 ###### + DRYADJ = # 標準 あり + #DRYADJ = nonstd/p2dad0.o # 乾燥対流調節なし + + ###### 浅い積雲対流 ###### + SHLCNV = # 標準 Tiedtke + #SHLCNV = nonstd/p2shl0.o # 浅い積雲なし + + ##### FFT ####### + FFT = # 標準 Numerical Receipes, 2^n のみ + #FFT = nonstd/ufftp.o # NCAR/Nakamura, 汎用 + + ##### 放射 ####### + RADIAT = # 標準 + # RADIAT = nonstd/p2rad0.o # 放射なし + + ##### 鉛直拡散 ####### + VDFUS = # 標準 + ###VDFUS = nonstd/p2vdf0.o # 拡散係数一定(バルク法はそのまま) + #VDFUS = nonstd/p2vdf0.o \ + nonstd/p2sfc0.o # 拡散係数一定 + + ##### 重力波抵抗 ####### + # GRVD = # 標準 + GRVD = nonstd/p2gra0.o # 重力波抵抗なし + + NONSTD = $(CUMLUS) $(DRYADJ) $(SHLCNV) $(FFT) \ + $(RADIAT) $(VDFUS) $(GRVD) + + + ####################################################################### + ######################## 実験設定部分ここまで ######################### + ####################################################################### + + ####################################################################### + ######################## サイト依存部分ここから ####################### + ####################################################################### + + ############## システム指定(SYSMAKE) ################### + + SYSTEM = Linux + ###SYSTEM = HITAC + ###SYSTEM = SX3 + + ############## コンパイルオプションの指定 ########## + + # STDFFLAGS = $(SYSFFLAGS) -O + ###STDFFLAGS = $(SYSFFLAGS) -g -u + # STDCFLAGS = $(SYSCFLAGS) -O + ###STDCFLAGS = $(SYSCFLAGS) -g + # STDLDFLAGS = $(SYSLDFLAGS) + ###STDLDFLAGS = $(SYSLDFLAGS) -g + STDFFLAGS = -g -u + STDCFLAGS = -g + STDLDFLAGS = -g + STDCPPFLAGS = $(SYSCPPFLAGS) -DREAL='ReaL*8' -DReal=DBLE + #AUTODBL = $(SYSAUTODBL) + #AUTODBLDEF = -DAUTO_DBL + AUTODBLDEF = -DAUTO_DBL + AUTODBL = + + ############## AGCM のソースツリー原点 ############ + DENNOUDIR = $(HOME)/e + AGCMDIR = $(HOME)/e/agcm5 + + ############## インストールディレクトリの指定 ########## + + ######## AGCM のインストールディレクトリ原点 + DESTDIR = $(AGCMDIR) + + ######## AGCM ライブラリのインストール先 + LIBDEST = $(DESTDIR)/lib + + ######## AGCM 実行ファイルインストール先 + BINDEST = $(DESTDIR)/bin + + ############## 他のライブラリのディレクトリの指定 ###### + + ######## dcl ライブラリのありか + DCLLIBDIR = /usr/lib/dcl-5.0.1 + + ######## symake ライブラリのありか + SYSMAKEDIR = /usr/lib/sysmake + + ######## gtool のライブラリのありか + GTOOLLIBDIR = /home/toyoda/e/gt3-dcl5/lib-g + #GTOOLLIBDIR = /usr/lib/gt3-dcl5 + GTOOLINCDIR = /usr/include/gt3-dcl5 + + ######## ispack smpack ライブラリのありか + ISPACKLIBDIR = $(DENNOUDIR)/lib + ISPACKLIBNAME = isp + ###SMPACKLIBNAME = sm + + ############## dcl が無い時はコメントを外す ######### + + ###NODCL = dclcompat/dclsubf.o # dcl がないとき用 + + ####################################################################### + ######################## サイト依存部分ここまで ####################### + ####################################################################### + + + ##################################################### + + include $(SYSMAKEDIR)/Makedef.$(SYSTEM) + + LIBDIR = $(AGCMDIR)/lib + BINDIR = $(AGCMDIR)/bin + + SRCDIR = $(AGCMDIR)/src + INCDIR = $(AGCMDIR)/src/include + + LIBNAME = agcm5$(RESOLUTION) + LIB = lib$(LIBNAME).a + LIBRARY = $(LIBDIR)/$(LIB) + LIBS = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAME) + + LIBNAMESM = agcm5sm$(RESOLUTION) + LIBSM = lib$(LIBNAMESM).a + LIBRARYSM = $(LIBDIR)/$(LIBSM) + LIBSSM = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAMESM) + + GTOOLLIBNAME = gtool3 + GTOOLLIB = lib$(GTOOLLIBNAME).a + GTOOLLIBFILE = $(GTOOLLIBDIR)/$(GTOOLLIB) + + DCLLIBNAME = dcl + DCLLIB = lib$(DCLLIBNAME).a + DCLLIBFILE = $(DCLLIBDIR)/$(DCLLIB) + + DCLDBASENAME = dcldbase + DCLDBASEDIR = $(DCLLIBDIR)/$(DCLDBASENAME) + + STDLIBS = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) + + STDLIBSSM = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) \ + $(LINKOPT)-L$(ISPACKLIBDIR) \ + $(LINKOPT)-l$(ISPACKLIBNAME) + + ### $(LINKOPT)-l$(SMPACKLIBNAME) \ + ### $(LINKOPT)-l$(ISPACKLIBNAME) + + MKINCLUDE = $(SRCDIR)/Mkinclude diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/Mkinclude.t42l16.adj agcm5.3-g77/src/Mkinclude.t42l16.adj *** agcm5.3/src/Mkinclude.t42l16.adj Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/Mkinclude.t42l16.adj Mon Jul 5 18:11:53 1999 *************** *** 0 **** --- 1,208 ---- + # include file for Makefile for $(DENNOUDIR)/model/agcm/agcm5 + # + # You need dcl, sysmake, gtool libraries to install agcm5. + # + ####################################################################### + ######################## 実験設定部分ここから ######################### + ####################################################################### + + ############## 分解能の指定 ######################## + + RESOLUTION = t42l16 + + ###RESOLUTION = t10l01 + ###RESOLUTION = t10l02 + ###RESOLUTION = t10l12 + ###RESOLUTION = t10l45 + ###RESOLUTION = t10l12xz + ###RESOLUTION = t10l45xz + ## RESOLUTION = t10l5 + ###RESOLUTION = t10l5sm + + ###RESOLUTION = t21l02 + ###RESOLUTION = t21l12 + ###RESOLUTION = t21l12yz + + ###RESOLUTION = t42l01 + ###RESOLUTION = t42l02 + ###RESOLUTION = t42l12 + ## RESOLUTION = t42l16 + ###RESOLUTION = t42l16yz + ###RESOLUTION = t42l23 + + ###RESOLUTION = t84l01 + ###RESOLUTION = t84l02 + ###RESOLUTION = t84l16 + + ###RESOLUTION = v00l12 + ###RESOLUTION = v00l16 + ###RESOLUTION = v00l45 + + ############# 実行ファイル名指定 agcm5.$(RESOLUTUON)$(MODELNAME) ##### + + MODELNAME = adj # 適当な名前をつけましょう. + + ############## usr/ (利用者設定ルーチン) の指定 ########### + + ###SPECIALDIR = usr/pochaka # ルーチンのおき場所の例です. + ###SPECIAL = usr/pochaka/p2tokubetu1.o usr/pochaka/p2tokubetu2.o + + ############## nonstd/ (非標準ルーチン) の指定 ########### + + ###### 積雲パラメタリゼーション ##### + ### CUMLUS = # 標準 Kuo 方式 + CUMLUS = nonstd/p2cuma.o # 対流調節方式 + #CUMLUS = nonstd/p2cum0.o # 積雲パラメタリゼーションなし + + ###### 乾燥対流調節 ###### + DRYADJ = # 標準 あり + #DRYADJ = nonstd/p2dad0.o # 乾燥対流調節なし + + ###### 浅い積雲対流 ###### + SHLCNV = # 標準 Tiedtke + #SHLCNV = nonstd/p2shl0.o # 浅い積雲なし + + ##### FFT ####### + FFT = # 標準 Numerical Receipes, 2^n のみ + #FFT = nonstd/ufftp.o # NCAR/Nakamura, 汎用 + + ##### 放射 ####### + RADIAT = # 標準 + # RADIAT = nonstd/p2rad0.o # 放射なし + + ##### 鉛直拡散 ####### + VDFUS = # 標準 + ###VDFUS = nonstd/p2vdf0.o # 拡散係数一定(バルク法はそのまま) + #VDFUS = nonstd/p2vdf0.o \ + nonstd/p2sfc0.o # 拡散係数一定 + + ##### 重力波抵抗 ####### + # GRVD = # 標準 + GRVD = nonstd/p2gra0.o # 重力波抵抗なし + + NONSTD = $(CUMLUS) $(DRYADJ) $(SHLCNV) $(FFT) \ + $(RADIAT) $(VDFUS) $(GRVD) + + + ####################################################################### + ######################## 実験設定部分ここまで ######################### + ####################################################################### + + ####################################################################### + ######################## サイト依存部分ここから ####################### + ####################################################################### + + ############## システム指定(SYSMAKE) ################### + + SYSTEM = Linux + ###SYSTEM = HITAC + ###SYSTEM = SX3 + + ############## コンパイルオプションの指定 ########## + + # STDFFLAGS = $(SYSFFLAGS) -O + ###STDFFLAGS = $(SYSFFLAGS) -g -u + # STDCFLAGS = $(SYSCFLAGS) -O + ###STDCFLAGS = $(SYSCFLAGS) -g + # STDLDFLAGS = $(SYSLDFLAGS) + ###STDLDFLAGS = $(SYSLDFLAGS) -g + STDFFLAGS = -g -u + STDCFLAGS = -g + STDLDFLAGS = -g + STDCPPFLAGS = $(SYSCPPFLAGS) -DREAL='ReaL*8' -DReal=DBLE + #AUTODBL = $(SYSAUTODBL) + #AUTODBLDEF = -DAUTO_DBL + AUTODBLDEF = -DAUTO_DBL + AUTODBL = + + ############## AGCM のソースツリー原点 ############ + DENNOUDIR = $(HOME)/e + AGCMDIR = $(HOME)/e/agcm5 + + ############## インストールディレクトリの指定 ########## + + ######## AGCM のインストールディレクトリ原点 + DESTDIR = $(AGCMDIR) + + ######## AGCM ライブラリのインストール先 + LIBDEST = $(DESTDIR)/lib + + ######## AGCM 実行ファイルインストール先 + BINDEST = $(DESTDIR)/bin + + ############## 他のライブラリのディレクトリの指定 ###### + + ######## dcl ライブラリのありか + DCLLIBDIR = /usr/lib/dcl-5.0.1 + + ######## symake ライブラリのありか + SYSMAKEDIR = /usr/lib/sysmake + + ######## gtool のライブラリのありか + GTOOLLIBDIR = /home/toyoda/e/gt3-dcl5/lib-g + #GTOOLLIBDIR = /usr/lib/gt3-dcl5 + GTOOLINCDIR = /usr/include/gt3-dcl5 + + ######## ispack smpack ライブラリのありか + ISPACKLIBDIR = $(DENNOUDIR)/lib + ISPACKLIBNAME = isp + ###SMPACKLIBNAME = sm + + ############## dcl が無い時はコメントを外す ######### + + ###NODCL = dclcompat/dclsubf.o # dcl がないとき用 + + ####################################################################### + ######################## サイト依存部分ここまで ####################### + ####################################################################### + + + ##################################################### + + include $(SYSMAKEDIR)/Makedef.$(SYSTEM) + + LIBDIR = $(AGCMDIR)/lib + BINDIR = $(AGCMDIR)/bin + + SRCDIR = $(AGCMDIR)/src + INCDIR = $(AGCMDIR)/src/include + + LIBNAME = agcm5$(RESOLUTION) + LIB = lib$(LIBNAME).a + LIBRARY = $(LIBDIR)/$(LIB) + LIBS = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAME) + + LIBNAMESM = agcm5sm$(RESOLUTION) + LIBSM = lib$(LIBNAMESM).a + LIBRARYSM = $(LIBDIR)/$(LIBSM) + LIBSSM = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAMESM) + + GTOOLLIBNAME = gtool3 + GTOOLLIB = lib$(GTOOLLIBNAME).a + GTOOLLIBFILE = $(GTOOLLIBDIR)/$(GTOOLLIB) + + DCLLIBNAME = dcl + DCLLIB = lib$(DCLLIBNAME).a + DCLLIBFILE = $(DCLLIBDIR)/$(DCLLIB) + + DCLDBASENAME = dcldbase + DCLDBASEDIR = $(DCLLIBDIR)/$(DCLDBASENAME) + + STDLIBS = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) + + STDLIBSSM = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) \ + $(LINKOPT)-L$(ISPACKLIBDIR) \ + $(LINKOPT)-l$(ISPACKLIBNAME) + + ### $(LINKOPT)-l$(SMPACKLIBNAME) \ + ### $(LINKOPT)-l$(ISPACKLIBNAME) + + MKINCLUDE = $(SRCDIR)/Mkinclude diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/Mkinclude.t42l16.adj.old agcm5.3-g77/src/Mkinclude.t42l16.adj.old *** agcm5.3/src/Mkinclude.t42l16.adj.old Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/Mkinclude.t42l16.adj.old Sun Jul 4 18:25:04 1999 *************** *** 0 **** --- 1,206 ---- + # include file for Makefile for $(DENNOUDIR)/model/agcm/agcm5 + # + # You need dcl, sysmake, gtool libraries to install agcm5. + # + ####################################################################### + ######################## 実験設定部分ここから ######################### + ####################################################################### + + ############## 分解能の指定 ######################## + + RESOLUTION = t42l16 + + ###RESOLUTION = t10l01 + ###RESOLUTION = t10l02 + ###RESOLUTION = t10l12 + ###RESOLUTION = t10l45 + ###RESOLUTION = t10l12xz + ###RESOLUTION = t10l45xz + ## RESOLUTION = t10l5 + ###RESOLUTION = t10l5sm + + ###RESOLUTION = t21l02 + ###RESOLUTION = t21l12 + ###RESOLUTION = t21l12yz + + ###RESOLUTION = t42l01 + ###RESOLUTION = t42l02 + ###RESOLUTION = t42l12 + ## RESOLUTION = t42l16 + ###RESOLUTION = t42l16yz + ###RESOLUTION = t42l23 + + ###RESOLUTION = t84l01 + ###RESOLUTION = t84l02 + ###RESOLUTION = t84l16 + + ###RESOLUTION = v00l12 + ###RESOLUTION = v00l16 + ###RESOLUTION = v00l45 + + ############# 実行ファイル名指定 agcm5.$(RESOLUTUON)$(MODELNAME) ##### + + MODELNAME = adj # 適当な名前をつけましょう. + + ############## usr/ (利用者設定ルーチン) の指定 ########### + + ###SPECIALDIR = usr/pochaka # ルーチンのおき場所の例です. + ###SPECIAL = usr/pochaka/p2tokubetu1.o usr/pochaka/p2tokubetu2.o + + ############## nonstd/ (非標準ルーチン) の指定 ########### + + ###### 積雲パラメタリゼーション ##### + ### CUMLUS = # 標準 Kuo 方式 + CUMLUS = nonstd/p2cuma.o # 対流調節方式 + #CUMLUS = nonstd/p2cum0.o # 積雲パラメタリゼーションなし + + ###### 乾燥対流調節 ###### + DRYADJ = # 標準 あり + #DRYADJ = nonstd/p2dad0.o # 乾燥対流調節なし + + ###### 浅い積雲対流 ###### + SHLCNV = # 標準 Tiedtke + #SHLCNV = nonstd/p2shl0.o # 浅い積雲なし + + ##### FFT ####### + FFT = # 標準 Numerical Receipes, 2^n のみ + #FFT = nonstd/ufftp.o # NCAR/Nakamura, 汎用 + + ##### 放射 ####### + RADIAT = # 標準 + # RADIAT = nonstd/p2rad0.o # 放射なし + + ##### 鉛直拡散 ####### + VDFUS = # 標準 + ###VDFUS = nonstd/p2vdf0.o # 拡散係数一定(バルク法はそのまま) + #VDFUS = nonstd/p2vdf0.o \ + nonstd/p2sfc0.o # 拡散係数一定 + + ##### 重力波抵抗 ####### + # GRVD = # 標準 + GRVD = nonstd/p2gra0.o # 重力波抵抗なし + + NONSTD = $(CUMLUS) $(DRYADJ) $(SHLCNV) $(FFT) \ + $(RADIAT) $(VDFUS) $(GRVD) + + + ####################################################################### + ######################## 実験設定部分ここまで ######################### + ####################################################################### + + ####################################################################### + ######################## サイト依存部分ここから ####################### + ####################################################################### + + ############## システム指定(SYSMAKE) ################### + + SYSTEM = Linux + ###SYSTEM = HITAC + ###SYSTEM = SX3 + + ############## コンパイルオプションの指定 ########## + + # STDFFLAGS = $(SYSFFLAGS) -O + ###STDFFLAGS = $(SYSFFLAGS) -g -u + # STDCFLAGS = $(SYSCFLAGS) -O + ###STDCFLAGS = $(SYSCFLAGS) -g + # STDLDFLAGS = $(SYSLDFLAGS) + ###STDLDFLAGS = $(SYSLDFLAGS) -g + STDFFLAGS = -g -u + STDCFLAGS = -g + STDLDFLAGS = -g + STDCPPFLAGS = $(SYSCPPFLAGS) + AUTODBL = $(SYSAUTODBL) + #AUTODBLDEF = -DAUTO_DBL + AUTODBLDEF = + + ############## AGCM のソースツリー原点 ############ + DENNOUDIR = $(HOME)/e + AGCMDIR = $(HOME)/e/agcm5 + + ############## インストールディレクトリの指定 ########## + + ######## AGCM のインストールディレクトリ原点 + DESTDIR = $(AGCMDIR) + + ######## AGCM ライブラリのインストール先 + LIBDEST = $(DESTDIR)/lib + + ######## AGCM 実行ファイルインストール先 + BINDEST = $(DESTDIR)/bin + + ############## 他のライブラリのディレクトリの指定 ###### + + ######## dcl ライブラリのありか + DCLLIBDIR = /usr/lib/dcl-5.0.1 + + ######## symake ライブラリのありか + SYSMAKEDIR = /usr/lib/sysmake + + ######## gtool のライブラリのありか + GTOOLLIBDIR = /usr/lib/gt3-dcl5 + GTOOLINCDIR = /usr/include/gt3-dcl5 + + ######## ispack smpack ライブラリのありか + ISPACKLIBDIR = $(DENNOUDIR)/lib + ISPACKLIBNAME = isp + ###SMPACKLIBNAME = sm + + ############## dcl が無い時はコメントを外す ######### + + ###NODCL = dclcompat/dclsubf.o # dcl がないとき用 + + ####################################################################### + ######################## サイト依存部分ここまで ####################### + ####################################################################### + + + ##################################################### + + include $(SYSMAKEDIR)/Makedef.$(SYSTEM) + + LIBDIR = $(AGCMDIR)/lib + BINDIR = $(AGCMDIR)/bin + + SRCDIR = $(AGCMDIR)/src + INCDIR = $(AGCMDIR)/src/include + + LIBNAME = agcm5$(RESOLUTION) + LIB = lib$(LIBNAME).a + LIBRARY = $(LIBDIR)/$(LIB) + LIBS = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAME) + + LIBNAMESM = agcm5sm$(RESOLUTION) + LIBSM = lib$(LIBNAMESM).a + LIBRARYSM = $(LIBDIR)/$(LIBSM) + LIBSSM = $(LINKOPT)-L$(LIBDIR) \ + $(LINKOPT)-l$(LIBNAMESM) + + GTOOLLIBNAME = gtool3 + GTOOLLIB = lib$(GTOOLLIBNAME).a + GTOOLLIBFILE = $(GTOOLLIBDIR)/$(GTOOLLIB) + + DCLLIBNAME = dcl + DCLLIB = lib$(DCLLIBNAME).a + DCLLIBFILE = $(DCLLIBDIR)/$(DCLLIB) + + DCLDBASENAME = dcldbase + DCLDBASEDIR = $(DCLLIBDIR)/$(DCLDBASENAME) + + STDLIBS = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) + + STDLIBSSM = $(LINKOPT)-L$(GTOOLLIBDIR) \ + $(LINKOPT)-l$(GTOOLLIBNAME) \ + $(LINKOPT)-L$(DCLLIBDIR) \ + $(LINKOPT)-l$(DCLLIBNAME) \ + $(LINKOPT)-L$(ISPACKLIBDIR) \ + $(LINKOPT)-l$(ISPACKLIBNAME) + + ### $(LINKOPT)-l$(SMPACKLIBNAME) \ + ### $(LINKOPT)-l$(ISPACKLIBNAME) + + MKINCLUDE = $(SRCDIR)/Mkinclude diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/aadmn.F agcm5.3-g77/src/aadmn.F *** agcm5.3/src/aadmn.F Tue Sep 23 04:56:23 1997 --- agcm5.3-g77/src/aadmn.F Mon Jul 5 15:47:27 1999 *************** *** 286,333 **** * [INTERN PARAM] INTEGER ISTYPL DATA ISTYPL / 1 / ! REAL VMISS * CALL GZDBGT ( 'MISS', VMISS ) !" 欠損値 * CALL HISTRG !" 出力の登録 I ( 'U ', 'u-velocity ' ,'m/s ', 'ALEV', ! I VMISS , VMISS , 2.5 , 10. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'V ', 'v-velocity ' ,'m/s ', 'ALEV', ! I VMISS , VMISS , 2.5 , 10. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'T ', 'temperature ' ,'K ', 'ALEV', ! I VMISS , VMISS , 2.5 , 10. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'PS ', 'surface pressure ' ,'mb ', 'ASFC', ! I VMISS , VMISS , 2.5 , 10. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'Q ', 'specific humidity ' ,'g/g ', 'ALEV', ! I VMISS , VMISS , 1.E-3 , 5.E-3 , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(3PF12.3)' ) * CALL HISTRG I ( 'VOR ', 'vorticity ' ,'1/s ', 'ALEV', ! I VMISS , VMISS , 1.E-6 , 5.E-6 , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(1PE12.3)' ) * CALL HISTRG I ( 'DIV ', 'divergence ' ,'1/s ', 'ALEV', ! I VMISS , VMISS , 1.E-6 , 5.E-6 , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(1PE12.3)' ) * CALL HISTRG I ( 'TG ', 'ground temperature ' ,'K ', 'ASFC', ! I VMISS , VMISS , 2.5 , 10. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * RETURN --- 286,339 ---- * [INTERN PARAM] INTEGER ISTYPL DATA ISTYPL / 1 / ! REAL VMISS, DIVSX, DIVLX * CALL GZDBGT ( 'MISS', VMISS ) !" 欠損値 * + DIVSX = 2.5 + DIVLX = 10. CALL HISTRG !" 出力の登録 I ( 'U ', 'u-velocity ' ,'m/s ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'V ', 'v-velocity ' ,'m/s ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'T ', 'temperature ' ,'K ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'PS ', 'surface pressure ' ,'mb ', 'ASFC', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * + DIVSX = 1.E-6 + DIVLX = 5.E-6 CALL HISTRG I ( 'Q ', 'specific humidity ' ,'g/g ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(3PF12.3)' ) * CALL HISTRG I ( 'VOR ', 'vorticity ' ,'1/s ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(1PE12.3)' ) * CALL HISTRG I ( 'DIV ', 'divergence ' ,'1/s ', 'ALEV', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(1PE12.3)' ) * + DIVSX = 2.5 + DIVLX = 10. CALL HISTRG I ( 'TG ', 'ground temperature ' ,'K ', 'ASFC', ! I VMISS, VMISS, DIVSX, DIVLX, ISTYPL, I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * RETURN diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/admin/asetl.F agcm5.3-g77/src/admin/asetl.F *** agcm5.3/src/admin/asetl.F Sat Sep 27 03:21:06 1997 --- agcm5.3-g77/src/admin/asetl.F Sun Jul 4 22:09:33 1999 *************** *** 16,23 **** * INTEGER JDIMD PARAMETER ( JDIMD=200 ) ! REAL * 8 QSINLA( JDIMD ) !" sin（緯度）：倍精度 ! REAL * 8 QGW ( JDIMD ) !" ガウス荷重 ：倍精度 * #ifdef SYS_IBMS INCLUDE (ZHDIM) --- 16,23 ---- * INTEGER JDIMD PARAMETER ( JDIMD=200 ) ! ReaL * 8 QSINLA( JDIMD ) !" sin（緯度）：倍精度 ! ReaL * 8 QGW ( JDIMD ) !" ガウス荷重 ：倍精度 * #ifdef SYS_IBMS INCLUDE (ZHDIM) *************** *** 125,132 **** * PI = ATAN( 1. ) *4. DO 1100 I = 1, IDIM ! ALON ( I ) = REAL( I-1 ) / REAL( IMAX ) * 2. * PI ! DLON ( I ) = 1. / REAL( IMAX ) 1100 CONTINUE * *" < 2. 軸の名称 > --- 125,132 ---- * PI = ATAN( 1. ) *4. DO 1100 I = 1, IDIM ! ALON ( I ) = Real( I-1 ) / Real( IMAX ) * 2. * PI ! DLON ( I ) = 1. / Real( IMAX ) 1100 CONTINUE * *" < 2. 軸の名称 > diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/autodbl.pl agcm5.3-g77/src/autodbl.pl *** agcm5.3/src/autodbl.pl Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/autodbl.pl Tue Jul 6 17:23:38 1999 *************** *** 0 **** --- 1,12 ---- + #!/usr/bin/perl + + while (<>) { + s/REAL\s*\(/Real\(/g; + next unless /^\s+REAL/i; + next unless /REAL\s*\*/i; + s/REAL/ReaL/; + } continue { + print; + } + + exit 0; diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics/dg2wd.F agcm5.3-g77/src/dynamics/dg2wd.F *** agcm5.3/src/dynamics/dg2wd.F Sun Apr 10 20:01:30 1994 --- agcm5.3-g77/src/dynamics/dg2wd.F Sun Jul 4 21:51:48 1999 *************** *** 40,46 **** * * [INTERNAL WORK] COMMON /COMWRK/ ! & WTKE REAL WTKE ( NMDIM , KMAX ) !" 運動エネルギー項 ＫＥ * INTEGER K, NM --- 40,49 ---- * * [INTERNAL WORK] COMMON /COMWRK/ ! & GDUU , GDVV ! REAL GDUU ( IDIM*JDIM, KMAX ) !" 西風 Ｕ ! REAL GDVV ( IDIM*JDIM, KMAX ) !" 南風 Ｖ ! EQUIVALENCE (GDUU, WTKE) REAL WTKE ( NMDIM , KMAX ) !" 運動エネルギー項 ＫＥ * INTEGER K, NM *************** *** 192,199 **** REAL GDQ ( IDIM*JDIM, KMAX ) !" 比湿 ｑ * * [INTERNAL WORK] COMMON /COMWRK/ ! & GDPI REAL GDPI ( IDIM*JDIM ) !" log Ps ＝π * INTEGER IJ --- 195,207 ---- REAL GDQ ( IDIM*JDIM, KMAX ) !" 比湿 ｑ * * [INTERNAL WORK] + * g77 が文句をいわないようにするため, + * サブルーチン GUV2DV と COMWRK の大きさを同じにしている COMMON /COMWRK/ ! & GDUU , GDVV ! EQUIVALENCE (GDUU, GDPI) ! REAL GDUU ( IDIM*JDIM, KMAX ) !" 西風 Ｕ ! REAL GDVV ( IDIM*JDIM, KMAX ) !" 南風 Ｖ REAL GDPI ( IDIM*JDIM ) !" log Ps ＝π * INTEGER IJ *************** *** 254,259 **** --- 262,269 ---- REAL UVFACT( IDIM*JDIM ) !" u→U のファクター * * [INTERNAL WORK] + * g77 が文句をいわないようにするため, + * サブルーチン GUV2DV と COMWRK の大きさを同じにしている COMMON /COMWRK/ & GDUU , GDVV REAL GDUU ( IDIM*JDIM, KMAX ) !" 西風 Ｕ diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics/dmtrx.F agcm5.3-g77/src/dynamics/dmtrx.F *** agcm5.3/src/dynamics/dmtrx.F Sun Apr 10 20:01:31 1994 --- agcm5.3-g77/src/dynamics/dmtrx.F Sun Jul 4 22:09:40 1999 *************** *** 50,56 **** * *" < 1. semi-implicit の行列 > * ! FLAPLA = - REAL( N ) * REAL( N+1 ) / ER**2 * DO 1000 M = 0, MIN(N,MMAX), MINT IF ( N-M .LE. LMAX ) THEN --- 50,56 ---- * *" < 1. semi-implicit の行列 > * ! FLAPLA = - Real( N ) * Real( N+1 ) / ER**2 * DO 1000 M = 0, MIN(N,MMAX), MINT IF ( N-M .LE. LMAX ) THEN diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics/dsetc.F agcm5.3-g77/src/dynamics/dsetc.F *** agcm5.3/src/dynamics/dsetc.F Sun Apr 10 20:01:31 1994 --- agcm5.3-g77/src/dynamics/dsetc.F Sun Jul 4 22:21:23 1999 *************** *** 43,50 **** INTEGER IEFOLD !" 最大波数減衰時間 INTEGER IDORDR !" 水平拡散のオーダー * ! REAL * 8 QSINLA( JDIM ) !" sin（緯度）：倍精度 ! REAL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 REAL PI !" 円周率π REAL OMEGA !" 自転角速度Ω * --- 43,50 ---- INTEGER IEFOLD !" 最大波数減衰時間 INTEGER IDORDR !" 水平拡散のオーダー * ! ReaL * 8 QSINLA( JDIM ) !" sin（緯度）：倍精度 ! ReaL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 REAL PI !" 円周率π REAL OMEGA !" 自転角速度Ω * *************** *** 78,84 **** *" < 3. コリオリ因子, Uファクター ) * PI = ATAN ( 1. )*4. ! OMEGA = 2. *PI / REAL( IDAYL ) * DO 3100 I = 1, IDIM DO 3110 J = 1, JDIM --- 78,84 ---- *" < 3. コリオリ因子, Uファクター ) * PI = ATAN ( 1. )*4. ! OMEGA = 2. *PI / Real( IDAYL ) * DO 3100 I = 1, IDIM DO 3110 J = 1, JDIM *************** *** 243,249 **** DO 1110 L = 0 , LEND AN = L + M IF ( AN .GT. 0.1 ) THEN ! DIFTN = ( AN*( AN+1. )/REAL(NMAX) )**(IDORDR/2) DIFMAX = MAX( DIFMAX, DIFTN ) DIFT( NMO(1,M,L) ) = DIFTN DIFT( NMO(2,M,L) ) = DIFTN --- 243,249 ---- DO 1110 L = 0 , LEND AN = L + M IF ( AN .GT. 0.1 ) THEN ! DIFTN = ( AN*( AN+1. )/Real(NMAX) )**(IDORDR/2) DIFMAX = MAX( DIFMAX, DIFTN ) DIFT( NMO(1,M,L) ) = DIFTN DIFT( NMO(2,M,L) ) = DIFTN *************** *** 253,266 **** * IF( NMAX .GT. 0 )THEN DO 1200 NM = 1 , NMDIM ! DIFV( NM ) = DIFT( NM ) - ( 2. /REAL(NMAX) )**(IDORDR/2) 1200 CONTINUE ENDIF *" !" 竹広 92/07/20 上記 IF 節追加 * *" < 2. 拡散係数を決める > * ! DIFVAL = 1. / REAL(IEFOLD) * *" < 3. 拡散のファクター > * --- 253,266 ---- * IF( NMAX .GT. 0 )THEN DO 1200 NM = 1 , NMDIM ! DIFV( NM ) = DIFT( NM ) - ( 2. /Real(NMAX) )**(IDORDR/2) 1200 CONTINUE ENDIF *" !" 竹広 92/07/20 上記 IF 節追加 * *" < 2. 拡散係数を決める > * ! DIFVAL = 1. / Real(IEFOLD) * *" < 3. 拡散のファクター > * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics/dsphe.F agcm5.3-g77/src/dynamics/dsphe.F *** agcm5.3/src/dynamics/dsphe.F Sun Apr 10 20:01:31 1994 --- agcm5.3-g77/src/dynamics/dsphe.F Sun Jul 4 21:51:48 1999 *************** *** 40,52 **** REAL ZDATA ( IDIM*JDIM, KMAX ) !" 東西スペクトル REAL WORK ( IDIM*JDIM, KMAX ) !" ワーク * ! REAL * 8 QSINLA( JDIM ) !" sin(緯度) (倍精度 ! REAL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 REAL ER !" 地球半径 INTEGER NMO ( 2, 0:MMAX, 0:LMAX ) !" スペクトルの添字順番 * ! REAL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! REAL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * INTEGER J, L, M, NMH * --- 40,52 ---- REAL ZDATA ( IDIM*JDIM, KMAX ) !" 東西スペクトル REAL WORK ( IDIM*JDIM, KMAX ) !" ワーク * ! ReaL * 8 QSINLA( JDIM ) !" sin(緯度) (倍精度 ! ReaL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 REAL ER !" 地球半径 INTEGER NMO ( 2, 0:MMAX, 0:LMAX ) !" スペクトルの添字順番 * ! ReaL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! ReaL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * INTEGER J, L, M, NMH * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics_sm/ds1.F agcm5.3-g77/src/dynamics_sm/ds1.F *** agcm5.3/src/dynamics_sm/ds1.F Sat Mar 7 23:23:29 1998 --- agcm5.3-g77/src/dynamics_sm/ds1.F Sun Jul 4 21:51:51 1999 *************** *** 18,43 **** ! [INTERNAL WORK FOR SMPACK] ! REAL*8 X (IMAX) !" 南北格子点座標値(radian) ! REAL*8 Y (JMAX) !" 東西格子点座標値(radian) ! REAL*8 W (JMAX) !" ガウス荷重 INTEGER IT (5) !" SMT???が用いる配列 ! REAL*8 T (IMAX*2) !" SMT???が用いる配列 INTEGER IP (2*MMAX+1) !" SMT???が用いる配列 ! REAL*8 P ((MMAX+1)*(MMAX+3)*JMAX/2) !" SMT???が用いる配列 ! REAL*8 QSG (JMAX/2) !" SMTS2Gが用いる配列 ! REAL*8 QGS (JMAX/2) !" SMTG2Sが用いる配列 ! REAL*8 QSV (JMAX/2) !" SMTS?Vが用いる配列 ! REAL*8 QVS (JMAX/2) !" SMTV?Sが用いる配列 ! REAL*8 QSU (JMAX/2) !" SMTS?Vが用いる配列 ! REAL*8 QUS (JMAX/2) !" SMTV?Sが用いる配列 ! REAL*8 R (MMAX*(MMAX-1)*2) !" SMTV?S, SMTS?Vが用いる配列 INTEGER ML (NMDIM) !" 東西波数 m のテーブル !" SMTV?S, SMTS?Vが用いる INTEGER NL (NMDIM) !" 全波数 n のテーブル ! REAL*8 RF (NMDIM) !" SMCLAPが用いる配列 ! REAL*8 RB (NMDIM) !" SMCLAPが用いる配列 SAVE X,Y,W,IT,T,IP,P SAVE QSG,QGS,QSV,QVS,QSU,QUS --- 18,43 ---- ! [INTERNAL WORK FOR SMPACK] ! ReaL*8 X (IMAX) !" 南北格子点座標値(radian) ! ReaL*8 Y (JMAX) !" 東西格子点座標値(radian) ! ReaL*8 W (JMAX) !" ガウス荷重 INTEGER IT (5) !" SMT???が用いる配列 ! ReaL*8 T (IMAX*2) !" SMT???が用いる配列 INTEGER IP (2*MMAX+1) !" SMT???が用いる配列 ! ReaL*8 P ((MMAX+1)*(MMAX+3)*JMAX/2) !" SMT???が用いる配列 ! ReaL*8 QSG (JMAX/2) !" SMTS2Gが用いる配列 ! ReaL*8 QGS (JMAX/2) !" SMTG2Sが用いる配列 ! ReaL*8 QSV (JMAX/2) !" SMTS?Vが用いる配列 ! ReaL*8 QVS (JMAX/2) !" SMTV?Sが用いる配列 ! ReaL*8 QSU (JMAX/2) !" SMTS?Vが用いる配列 ! ReaL*8 QUS (JMAX/2) !" SMTV?Sが用いる配列 ! ReaL*8 R (MMAX*(MMAX-1)*2) !" SMTV?S, SMTS?Vが用いる配列 INTEGER ML (NMDIM) !" 東西波数 m のテーブル !" SMTV?S, SMTS?Vが用いる INTEGER NL (NMDIM) !" 全波数 n のテーブル ! ReaL*8 RF (NMDIM) !" SMCLAPが用いる配列 ! ReaL*8 RB (NMDIM) !" SMCLAPが用いる配列 SAVE X,Y,W,IT,T,IP,P SAVE QSG,QGS,QSV,QVS,QSU,QUS *************** *** 51,59 **** ! [INTERNAL WORK] CHARACTER HFUNC*4 ! REAL*8 WORK ( IDIM*JDIM*3 ) !" ワーク ! REAL*8 GW ( JMAX ) !" GW ( ガウス荷重 ) ! REAL*8 Y0 ( JMAX ) !" 東西格子点座標値(radian) LOGICAL OFIRST DATA OFIRST /.TRUE./ --- 51,59 ---- ! [INTERNAL WORK] CHARACTER HFUNC*4 ! ReaL*8 WORK ( IDIM*JDIM*3 ) !" ワーク ! ReaL*8 GW ( JMAX ) !" GW ( ガウス荷重 ) ! ReaL*8 Y0 ( JMAX ) !" 東西格子点座標値(radian) LOGICAL OFIRST DATA OFIRST /.TRUE./ diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/dynamics_sm/dsm.F agcm5.3-g77/src/dynamics_sm/dsm.F *** agcm5.3/src/dynamics_sm/dsm.F Sat Mar 7 23:22:58 1998 --- agcm5.3-g77/src/dynamics_sm/dsm.F Sun Jul 4 21:51:51 1999 *************** *** 18,52 **** ! [INTERNAL WORK FOR SMPACK] ! REAL*8 X (IMAX) !" 南北格子点座標値(radian) ! REAL*8 Y (JMAX) !" 東西格子点座標値(radian) ! REAL*8 W (JMAX) !" ガウス荷重 INTEGER IT (5) !" SMT???が用いる配列 ! REAL*8 T (IMAX*2) !" SMT???が用いる配列 INTEGER IP (KMAX*(2*MMAX+1)) !" SMT???が用いる配列 ! REAL*8 P (KMAX*(MMAX+1)*(MMAX+3)*JMAX/2) !" SMT???が用いる配列 ! REAL*8 QSG (JMAX/2) !" SMTS2Gが用いる配列 ! REAL*8 QGS (JMAX/2) !" SMTG2Sが用いる配列 ! REAL*8 QSV (JMAX/2) !" SMTS?Vが用いる配列 ! REAL*8 QVS (JMAX/2) !" SMTV?Sが用いる配列 ! REAL*8 QSU (JMAX/2) !" SMTS?Vが用いる配列 ! REAL*8 QUS (JMAX/2) !" SMTV?Sが用いる配列 ! REAL*8 R (KMAX*MMAX*(MMAX-1)*2)!" SMTV?S,SMTS?Vが用いる配列 INTEGER ML (NMDIM) !" 東西波数 m のテーブル !" SMTV?S, SMTS?Vが用いる INTEGER NL (NMDIM) !" 全波数 n のテーブル ! REAL*8 RF (NMDIM) !" SMCLAPが用いる配列 ! REAL*8 RB (NMDIM) !" SMCLAPが用いる配列 SAVE X,Y,W,IT,T,IP,P SAVE QSG,QGS,QSV,QVS,QSU,QUS SAVE R,ML,NL,RF,RB ! REAL*8 QSVMER(JMAX/2) !" QSV*ER ! REAL*8 QVSDER(JMAX/2) !" QVS/ER ! REAL*8 QSUMER(JMAX/2) !" QSU*ER ! REAL*8 QUSDER(JMAX/2) !" QUS/ER SAVE QSVMER, QVSDER, QSUMER, QUSDER ! [INPUT or OUTPUT] --- 18,52 ---- ! [INTERNAL WORK FOR SMPACK] ! ReaL*8 X (IMAX) !" 南北格子点座標値(radian) ! ReaL*8 Y (JMAX) !" 東西格子点座標値(radian) ! ReaL*8 W (JMAX) !" ガウス荷重 INTEGER IT (5) !" SMT???が用いる配列 ! ReaL*8 T (IMAX*2) !" SMT???が用いる配列 INTEGER IP (KMAX*(2*MMAX+1)) !" SMT???が用いる配列 ! ReaL*8 P (KMAX*(MMAX+1)*(MMAX+3)*JMAX/2) !" SMT???が用いる配列 ! ReaL*8 QSG (JMAX/2) !" SMTS2Gが用いる配列 ! ReaL*8 QGS (JMAX/2) !" SMTG2Sが用いる配列 ! ReaL*8 QSV (JMAX/2) !" SMTS?Vが用いる配列 ! ReaL*8 QVS (JMAX/2) !" SMTV?Sが用いる配列 ! ReaL*8 QSU (JMAX/2) !" SMTS?Vが用いる配列 ! ReaL*8 QUS (JMAX/2) !" SMTV?Sが用いる配列 ! ReaL*8 R (KMAX*MMAX*(MMAX-1)*2)!" SMTV?S,SMTS?Vが用いる配列 INTEGER ML (NMDIM) !" 東西波数 m のテーブル !" SMTV?S, SMTS?Vが用いる INTEGER NL (NMDIM) !" 全波数 n のテーブル ! ReaL*8 RF (NMDIM) !" SMCLAPが用いる配列 ! ReaL*8 RB (NMDIM) !" SMCLAPが用いる配列 SAVE X,Y,W,IT,T,IP,P SAVE QSG,QGS,QSV,QVS,QSU,QUS SAVE R,ML,NL,RF,RB ! ReaL*8 QSVMER(JMAX/2) !" QSV*ER ! ReaL*8 QVSDER(JMAX/2) !" QVS/ER ! ReaL*8 QSUMER(JMAX/2) !" QSU*ER ! ReaL*8 QUSDER(JMAX/2) !" QUS/ER SAVE QSVMER, QVSDER, QSUMER, QUSDER ! [INPUT or OUTPUT] *************** *** 59,69 **** INTEGER J, K, NM CHARACTER HFUNC*4 ! REAL*8 WORKA( NMDIM,KMAX ) !" ワーク ! REAL*8 WORKB( NMDIM,KMAX ) !" ワーク ! REAL*8 WORK ( IDIM*JDIM*KMAX*3 ) !" ワーク ! REAL*8 GW ( JMAX ) !" GW ( ガウス荷重 ) ! REAL*8 Y0 ( JMAX ) !" 東西格子点座標値(radian) LOGICAL OFIRST DATA OFIRST /.TRUE./ --- 59,69 ---- INTEGER J, K, NM CHARACTER HFUNC*4 ! ReaL*8 WORKA( NMDIM,KMAX ) !" ワーク ! ReaL*8 WORKB( NMDIM,KMAX ) !" ワーク ! ReaL*8 WORK ( IDIM*JDIM*KMAX*3 ) !" ワーク ! ReaL*8 GW ( JMAX ) !" GW ( ガウス荷重 ) ! ReaL*8 Y0 ( JMAX ) !" 東西格子点座標値(radian) LOGICAL OFIRST DATA OFIRST /.TRUE./ diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/etc/txt2bin.F agcm5.3-g77/src/etc/txt2bin.F *** agcm5.3/src/etc/txt2bin.F Sun Apr 10 20:01:41 1994 --- agcm5.3-g77/src/etc/txt2bin.F Sun Jul 4 21:51:50 1999 *************** *** 1,7 **** PROGRAM TX2BIN PARAMETER ( IMAX=128, JMAX=64 ) CHARACTER HEAD ( 64 )*16 ! REAL*8 GDATA( IMAX, JMAX ) CHARACTER HFILE*80 * READ ( 5,* ) HFILE --- 1,7 ---- PROGRAM TX2BIN PARAMETER ( IMAX=128, JMAX=64 ) CHARACTER HEAD ( 64 )*16 ! ReaL*8 GDATA( IMAX, JMAX ) CHARACTER HFILE*80 * READ ( 5,* ) HFILE diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/gtintfc/gintf.F agcm5.3-g77/src/gtintfc/gintf.F *** agcm5.3/src/gtintfc/gintf.F Sat Sep 27 03:59:13 1997 --- agcm5.3-g77/src/gtintfc/gintf.F Mon Jul 5 14:39:17 1999 *************** *** 21,28 **** CHARACTER HHEAD ( * )*(*) !" ヘッダー * #ifdef AUTO_DBL ! REAL * 8 GDATA ( IDIM, JDIM, * ) !" データ ! REAL * 8 GDATAX( * ) !" データ #else REAL GDATA ( IDIM, JDIM, * ) !" データ REAL GDATAX( * ) !" データ --- 21,28 ---- CHARACTER HHEAD ( * )*(*) !" ヘッダー * #ifdef AUTO_DBL ! ReaL * 8 GDATA ( IDIM, JDIM, * ) !" データ ! ReaL * 8 GDATAX( * ) !" データ #else REAL GDATA ( IDIM, JDIM, * ) !" データ REAL GDATAX( * ) !" データ *************** *** 46,55 **** * [INTERNAL WORK] INTEGER KDIMD PARAMETER (KDIMD=KMAX+1) ! REAL * 8 DDATA ( IMAX, JMAX, KDIMD ) ! REAL * 4 SDATA ( IMAX, JMAX, KDIMD ) ! REAL * 8 DDATAX( IMAX*JMAX*KDIMD ) ! REAL * 4 SDATAX( IMAX*JMAX*KDIMD ) EQUIVALENCE (DDATA,SDATA,DDATAX,SDATAX) CHARACTER HDFMTI *(NCC) !" データフォーマット INTEGER I, J, K, KLEVS, ISIZE --- 46,55 ---- * [INTERNAL WORK] INTEGER KDIMD PARAMETER (KDIMD=KMAX+1) ! ReaL * 8 DDATA ( IMAX, JMAX, KDIMD ) ! ReaL * 4 SDATA ( IMAX, JMAX, KDIMD ) ! ReaL * 8 DDATAX( IMAX*JMAX*KDIMD ) ! ReaL * 4 SDATAX( IMAX*JMAX*KDIMD ) EQUIVALENCE (DDATA,SDATA,DDATAX,SDATAX) CHARACTER HDFMTI *(NCC) !" データフォーマット INTEGER I, J, K, KLEVS, ISIZE *************** *** 216,222 **** CHARACTER HFMT *(*) !" フォーマット * #ifdef AUTO_DBL ! REAL*8 RVAL !" 実数(単/倍精度)データ #else REAL RVAL !" 実数(単/倍精度)データ #endif --- 216,222 ---- CHARACTER HFMT *(*) !" フォーマット * #ifdef AUTO_DBL ! ReaL*8 RVAL !" 実数(単/倍精度)データ #else REAL RVAL !" 実数(単/倍精度)データ #endif *************** *** 226,232 **** INTEGER NCHR !" 文字列長 * * [INTERNAL WORK] ! REAL * 4 SVAL * SVAL = RVAL * --- 226,232 ---- INTEGER NCHR !" 文字列長 * * [INTERNAL WORK] ! ReaL * 4 SVAL * SVAL = RVAL * *************** *** 249,261 **** * * [OUTPUT] #ifdef AUTO_DBL ! REAL * 8 RT !" 出力時刻：単位 HUNITY(単/倍) #else REAL RT !" 出力時刻：単位 HUNITY(単/倍) #endif * * [INTERNAL WORK] ! REAL * 4 ST * CALL GUCT2R !" 時間単位の変換 HUNIT → HUNITY: 単精度 I ( KT , HUNIT , HUNITY, IDELT , --- 249,261 ---- * * [OUTPUT] #ifdef AUTO_DBL ! ReaL * 8 RT !" 出力時刻：単位 HUNITY(単/倍) #else REAL RT !" 出力時刻：単位 HUNITY(単/倍) #endif * * [INTERNAL WORK] ! ReaL * 4 ST * CALL GUCT2R !" 時間単位の変換 HUNIT → HUNITY: 単精度 I ( KT , HUNIT , HUNITY, IDELT , *************** *** 307,317 **** CHARACTER HCSIGN *(*) !" データ作成者 * #ifdef AUTO_DBL ! REAL * 8 VMISS !" 欠損値の値 ! REAL * 8 DMIN !" レンジ(最小) ! REAL * 8 DMAX !" レンジ(最大) ! REAL * 8 DIVS !" 間隔(小) ! REAL * 8 DIVL !" 間隔(大) #else REAL VMISS !" 欠損値の値 REAL DMIN !" レンジ(最小) --- 307,317 ---- CHARACTER HCSIGN *(*) !" データ作成者 * #ifdef AUTO_DBL ! ReaL * 8 VMISS !" 欠損値の値 ! ReaL * 8 DMIN !" レンジ(最小) ! ReaL * 8 DMAX !" レンジ(最大) ! ReaL * 8 DIVS !" 間隔(小) ! ReaL * 8 DIVL !" 間隔(大) #else REAL VMISS !" 欠損値の値 REAL DMIN !" レンジ(最小) *************** *** 319,333 **** REAL DIVS !" 間隔(小) REAL DIVL !" 間隔(大) #endif ! * : 単/倍精度 * * [INTERNAL WORK] ! REAL SVMISS !" 欠損値の値 ! REAL SDMIN !" レンジ(最小) ! REAL SDMAX !" レンジ(最大) ! REAL SDIVS !" 間隔(小) ! REAL SDIVL !" 間隔(大) ! * : 単精度 * SVMISS = VMISS SDMIN = DMIN --- 319,334 ---- REAL DIVS !" 間隔(小) REAL DIVL !" 間隔(大) #endif ! * ! * gtool のための単精度変数 * * [INTERNAL WORK] ! ReaL*4 SVMISS !" 欠損値の値 ! ReaL*4 SDMIN !" レンジ(最小) ! ReaL*4 SDMAX !" レンジ(最大) ! ReaL*4 SDIVS !" 間隔(小) ! ReaL*4 SDIVL !" 間隔(大) ! * : 単精度に変換 * SVMISS = VMISS SDMIN = DMIN *************** *** 368,378 **** INTEGER ISTYP !" スケーリングタイプ * #ifdef AUTO_DBL ! REAL * 8 VMISS !" 欠損値の値 ! REAL * 8 DMIN !" レンジ(最小) ! REAL * 8 DMAX !" レンジ(最大) ! REAL * 8 DIVS !" 間隔(小) ! REAL * 8 DIVL !" 間隔(大) #else REAL VMISS !" 欠損値の値 REAL DMIN !" レンジ(最小) --- 369,379 ---- INTEGER ISTYP !" スケーリングタイプ * #ifdef AUTO_DBL ! ReaL * 8 VMISS !" 欠損値の値 ! ReaL * 8 DMIN !" レンジ(最小) ! ReaL * 8 DMAX !" レンジ(最大) ! ReaL * 8 DIVS !" 間隔(小) ! ReaL * 8 DIVL !" 間隔(大) #else REAL VMISS !" 欠損値の値 REAL DMIN !" レンジ(最小) *************** *** 416,429 **** * * [OUTPUT] #ifdef AUTO_DBL ! REAL * 8 RPARA !" 単/倍精度 #else REAL RPARA !" 単/倍精度 #endif * * [INTERNAL WORK] ! REAL * 4 SPARA !" 単精度 * CALL GTPGET !" GTOOL 共通パラメータ（単精度）を参照 I ( HP, --- 417,430 ---- * * [OUTPUT] #ifdef AUTO_DBL ! ReaL * 8 RPARA !" 単/倍精度 #else REAL RPARA !" 単/倍精度 #endif * * [INTERNAL WORK] ! ReaL * 4 SPARA !" 単精度 * CALL GTPGET !" GTOOL 共通パラメータ（単精度）を参照 I ( HP, diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/include/zcdim.F agcm5.3-g77/src/include/zcdim.F *** agcm5.3/src/include/zcdim.F Sat Mar 7 22:23:07 1998 --- agcm5.3-g77/src/include/zcdim.F Mon Jul 5 18:11:10 1999 *************** *** 1,6 **** ! * PACKAGE ZGT10L5 !" 格子点数，波数 * ! *" [HIS] 98/03/06(akahori) * INTEGER IMAX !" 経度格子数 INTEGER JMAX !" 緯度格子数 --- 1,7 ---- ! * PACKAGE ZGT42L16 !" 格子点数，波数 （Ｔ42Ｌ16） * ! *" [HIS] 90/08/31(numaguti) ! *" 92/08/11(takepiro) * INTEGER IMAX !" 経度格子数 INTEGER JMAX !" 緯度格子数 *************** *** 16,25 **** INTEGER JMXHF !" JMAX/2+1 INTEGER NWORK !" ワーク領域大きさ * ! PARAMETER ( IMAX = 32, JMAX = 16, & IDIM = IMAX+1, JDIM = JMAX, ! & MMAX = 10, LMAX = MMAX, NMAX = LMAX, MINT = 1, ! & KMAX = 5, KDIM = KMAX, ! & NMDIM = (MMAX+1)*(MMAX+1), & JMXHF = JMAX/2+1, ! & NWORK = IDIM*JDIM*KMAX*12 ) --- 17,27 ---- INTEGER JMXHF !" JMAX/2+1 INTEGER NWORK !" ワーク領域大きさ * ! PARAMETER ( IMAX = 128, JMAX = 64, & IDIM = IMAX+1, JDIM = JMAX, ! & MMAX = 42, LMAX = 42, NMAX = LMAX, MINT = 1, ! & KMAX = 16, KDIM = KMAX, ! & NMDIM = (MMAX/MINT+1)*(2*(NMAX+1)-MMAX) ! & - (NMAX-LMAX)/MINT*(NMAX-LMAX+1), & JMXHF = JMAX/2+1, ! & NWORK = IDIM*JDIM*KMAX*12 ) diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/io/idiag.F agcm5.3-g77/src/io/idiag.F *** agcm5.3/src/io/idiag.F Sun Apr 10 20:01:40 1994 --- agcm5.3-g77/src/io/idiag.F Sun Jul 4 21:51:49 1999 *************** *** 220,227 **** HTITLE = HTTL WRITE ( JFLOUT, 6010 ) HTITLE, VMAX, IPMAX, JPMAX, KPMAX WRITE ( JFLOUT, 6020 ) VMIN, IPMIN, JPMIN, KPMIN ! 6010 FORMAT ( ' ***[',A16,']*** MAX= ',D12.4,' AT 'I4,','I4,','I4 ) ! 6020 FORMAT ( ' ',16X,' MIN= ',D12.4,' AT 'I4,','I4,','I4 ) * RETURN END --- 220,227 ---- HTITLE = HTTL WRITE ( JFLOUT, 6010 ) HTITLE, VMAX, IPMAX, JPMAX, KPMAX WRITE ( JFLOUT, 6020 ) VMIN, IPMIN, JPMIN, KPMIN ! 6010 FORMAT ( ' ***[',A16,']*** MAX= ',D12.4,' AT ',I4,',',I4,',',I4 ) ! 6020 FORMAT ( ' ',16X,' MIN= ',D12.4,' AT ',I4,',',I4,',',I4 ) * RETURN END diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/io/ihead.F agcm5.3-g77/src/io/ihead.F *** agcm5.3/src/io/ihead.F Sun Apr 10 20:01:40 1994 --- agcm5.3-g77/src/io/ihead.F Mon Jul 5 12:42:25 1999 *************** *** 9,15 **** I DMIN , DMAX , DIVS , DIVL , ISTYP ) * * [OUTPUT] ! CHARACTER HHEAD ( * )*(*) !" ヘッダー INTEGER ISIZE !" サイズ * * [INPUT] --- 9,15 ---- I DMIN , DMAX , DIVS , DIVL , ISTYP ) * * [OUTPUT] ! * CHARACTER HHEAD ( * )*(*) !" ヘッダー INTEGER ISIZE !" サイズ * * [INPUT] *************** *** 34,39 **** --- 34,40 ---- #include "zcdim.F" !" 格子点数，波数 #include "zhdim.F" !" 文字列文字数 #endif + CHARACTER HHEAD (NDC)*(NCC) !" ヘッダー * [INTERNAL WORK] CHARACTER HAITMX *(NCC) !" 第1軸名称 CHARACTER HAITMY *(NCC) !" 第2軸名称 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/io/ihist.F agcm5.3-g77/src/io/ihist.F *** agcm5.3/src/io/ihist.F Sun Apr 10 20:01:40 1994 --- agcm5.3-g77/src/io/ihist.F Mon Jul 5 11:56:20 1999 *************** *** 134,140 **** * IF ( IMEANH( IH ) .GT. 0 ) THEN DO 1100 ID = INDEXH( IH ), INDEXH( IH )+ISIZEH( IH )-1 ! GOUT ( ID ) = GOUT ( ID ) / REAL( ITSUM ( IH ) ) 1100 CONTINUE ENDIF * --- 134,140 ---- * IF ( IMEANH( IH ) .GT. 0 ) THEN DO 1100 ID = INDEXH( IH ), INDEXH( IH )+ISIZEH( IH )-1 ! GOUT ( ID ) = GOUT ( ID ) / Real( ITSUM ( IH ) ) 1100 CONTINUE ENDIF * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/physics/p2admn.F agcm5.3-g77/src/physics/p2admn.F *** agcm5.3/src/physics/p2admn.F Sat Sep 27 03:50:01 1997 --- agcm5.3-g77/src/physics/p2admn.F Sun Jul 4 22:18:35 1999 *************** *** 337,358 **** INTEGER ISTYPL DATA ISTYPL / 1 / REAL VMISS * CALL GZDBGT( 'MISS', VMISS ) * CALL HISTRG I ( 'RAIN' , 'precipitation ' ,'W/m**2', 'ASFC', ! I 0. , VMISS , 200. , 1000. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'EVAP' , 'latent heat flux ' ,'W/m**2', 'ASFC', ! I VMISS , VMISS , 20. , 100. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG I ( 'SENS' , 'sensible heat flux ' ,'W/m**2', 'ASFC', ! I VMISS , VMISS , 20. , 100. , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG --- 337,366 ---- INTEGER ISTYPL DATA ISTYPL / 1 / REAL VMISS + REAL DMINX, DMAXX, DIVSX, DIVLX * CALL GZDBGT( 'MISS', VMISS ) * + DMINX = 0. + DIVSX = 200. + DIVLX = 1000. CALL HISTRG I ( 'RAIN' , 'precipitation ' ,'W/m**2', 'ASFC', ! I DMINX , VMISS , DIVSX , DIVLX , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * + DIVSX = 20. + DIVLX = 100. CALL HISTRG I ( 'EVAP' , 'latent heat flux ' ,'W/m**2', 'ASFC', ! I VMISS , VMISS , DIVSX , DIVLX , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * + DIVSX = 20. + DIVLX = 100. CALL HISTRG I ( 'SENS' , 'sensible heat flux ' ,'W/m**2', 'ASFC', ! I VMISS , VMISS , DIVSX , DIVLX , ISTYPL , I ' ' , ' ' , 0 , 0 , 'X' ,'(F12.3)' ) * CALL HISTRG diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/physics/p2cumk.F agcm5.3-g77/src/physics/p2cumk.F *** agcm5.3/src/physics/p2cumk.F Sun Apr 10 20:01:29 1994 --- agcm5.3-g77/src/physics/p2cumk.F Sun Jul 4 21:51:47 1999 *************** *** 10,17 **** INCLUDE (ZCDIM) !' 格子点数，波数 INCLUDE (ZCCOM) !' 標準物理定数 #else ! #include "zcdim.F" !' 格子点数，波数 ! #include "zccom.F" !' 標準物理定数 #endif * * [MODIFY] --- 10,17 ---- INCLUDE (ZCDIM) !' 格子点数，波数 INCLUDE (ZCCOM) !' 標準物理定数 #else ! #include "zcdim.F" !" 格子点数，波数 ! #include "zccom.F" !" 標準物理定数 #endif * * [MODIFY] *************** *** 67,73 **** #ifdef SYS_IBMS INCLUDE (ZQSAT) !' 飽和水蒸気圧 #else ! #include "zqsat.F" !' 飽和水蒸気圧 #endif * IF ( OFIRST ) THEN --- 67,73 ---- #ifdef SYS_IBMS INCLUDE (ZQSAT) !' 飽和水蒸気圧 #else ! #include "zqsat.F" !" 飽和水蒸気圧 #endif * IF ( OFIRST ) THEN diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/prep/xmkinit.F agcm5.3-g77/src/prep/xmkinit.F *** agcm5.3/src/prep/xmkinit.F Sun Apr 10 20:01:38 1994 --- agcm5.3-g77/src/prep/xmkinit.F Mon Jul 5 10:35:16 1999 *************** *** 75,80 **** --- 75,81 ---- WRITE ( 6,* ) ' ENTER NAMELST NMDATE?' READ ( 5,NMDATE, END=1910 ) 1910 CONTINUE + CALL INIINI * CALL GTCSET( 'DDSET' , HDSET ) * *************** *** 128,133 **** --- 129,137 ---- * [INTERNAL PARM] REAL TPRTRB DATA TPRTRB / 0.1 / + REAL TINIT, PINIT, QINIT + DATA TINIT, PINIT, QINIT / 250.0, 1000.0, 1.0E-10 / + NAMELIST /NMINIP/ TPRTRB, TINIT, PINIT, QINIT * * [INTERNAL WORK] REAL RANDOM( IDIM, JDIM, KMAX ) *************** *** 147,159 **** DO 1000 I = 1, IDIM GAU ( I,J,K ) = 0. GAV ( I,J,K ) = 0. ! GAT ( I,J,K ) = 250. ! GAQ ( I,J,K ) = 1.E-10 1000 CONTINUE * DO 2000 J = 1, JDIM DO 2000 I = 1, IDIM ! GAPS ( I,J ) = 1000. 2000 CONTINUE * *" << RANDOM : 乱数の付加 >> --- 151,163 ---- DO 1000 I = 1, IDIM GAU ( I,J,K ) = 0. GAV ( I,J,K ) = 0. ! GAT ( I,J,K ) = TINIT ! GAQ ( I,J,K ) = QINIT 1000 CONTINUE * DO 2000 J = 1, JDIM DO 2000 I = 1, IDIM ! GAPS ( I,J ) = PINIT 2000 CONTINUE * *" << RANDOM : 乱数の付加 >> *************** *** 163,170 **** DO 3000 J = 1, JDIM DO 3000 I = 1, IDIM RANDOM( I,J,K ) ! & = SIN( REAL(I*J*(KMAX-K)) ! & / REAL(IMAX*JMAX*KMAX)*10. ) 3000 CONTINUE * DO 3100 K = 1, KMAX --- 167,174 ---- DO 3000 J = 1, JDIM DO 3000 I = 1, IDIM RANDOM( I,J,K ) ! & = SIN( Real(I*J*(KMAX-K)) ! & / Real(IMAX*JMAX*KMAX)*10. ) 3000 CONTINUE * DO 3100 K = 1, KMAX *************** *** 174,178 **** --- 178,192 ---- & + RANDOM( I,J,K ) * TPRTRB 3100 CONTINUE * + RETURN + + * 初期値に含まれる定数の設定 + ENTRY INIINI + + WRITE(6, NMINIP) + WRITE(6, *) ' ENTER NAMELST NMINIP:' + READ(5, NMINIP, END=4000) + WRITE(6, NMINIP) + 4000 CONTINUE RETURN END diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/Makefile agcm5.3-g77/src/sysdep/Makefile *** agcm5.3/src/sysdep/Makefile Thu Mar 5 18:58:22 1998 --- agcm5.3-g77/src/sysdep/Makefile Wed Feb 10 21:57:58 1999 *************** *** 19,32 **** SRCS = y$(SYSTEM).F ! all: lib ! lib: $(OBJS) ! $(AR) $(LIBRARY) $(OBJS) $(RANLIB) $(LIBRARY) libsm: $(OBJS) ! $(AR) $(LIBRARYSM) $(OBJS) $(RANLIB) $(LIBRARYSM) clean: --- 19,38 ---- SRCS = y$(SYSTEM).F ! # hack by toyoda: service routine emulations written in C ! # ! # --- if you need not it, just let $(CSERV) empty --- ! # CSERV = ! CSERV = g77linux.o ! all: lib $(CSERV) ! ! lib: $(OBJS) $(CSERV) ! $(AR) $(LIBRARY) $(OBJS) $(CSERV) $(RANLIB) $(LIBRARY) libsm: $(OBJS) ! $(AR) $(LIBRARYSM) $(OBJS) $(CSERV) $(RANLIB) $(LIBRARYSM) clean: *************** *** 55,57 **** --- 61,64 ---- $(TOUCH) $(PACKFILE) ### + diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/g77linux.c agcm5.3-g77/src/sysdep/g77linux.c *** agcm5.3/src/sysdep/g77linux.c Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/sysdep/g77linux.c Tue Jul 6 17:45:41 1999 *************** *** 0 **** --- 1,46 ---- + /* Copyright by TOYODA Eizi, 1999. All rights reserved. + * g77.c: Fortran service routine emulation for g77 + UNIX environment. + */ + + #include + #include + #include + + double etime_(float *tarray) { + struct tms etm; + + (void)times(&etm); + tarray[0] = (double)etm.tms_utime / CLOCKS_PER_SEC; + tarray[1] = (double)etm.tms_stime / CLOCKS_PER_SEC; + return tarray[0] + tarray[1]; + } + + int idate_(int *ans) { + struct tm *tmp; + time_t now; + time(&now); + tmp = localtime(&now); + ans[0] = tmp->tm_mday; + ans[1] = tmp->tm_mon; + ans[2] = tmp->tm_year + 1900; + return 0; + } + + int itime_(int *ans) { + struct tm *tmp; + time_t now; + time(&now); + tmp = localtime(&now); + ans[0] = tmp->tm_hour; + ans[1] = tmp->tm_min; + ans[2] = tmp->tm_sec; + return 0; + } + + int getpid_(void) { + return getpid(); + } + + int kill_(int pid, int signal) { + return kill(pid, signal); + } diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/yCONVEX.F agcm5.3-g77/src/sysdep/yCONVEX.F *** agcm5.3/src/sysdep/yCONVEX.F Wed Sep 24 05:31:51 1997 --- agcm5.3-g77/src/sysdep/yCONVEX.F Sun Jul 4 21:51:49 1999 *************** *** 60,69 **** SAVE * * [INTERNAL WORK] ! REAL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! REAL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 --- 60,69 ---- SAVE * * [INTERNAL WORK] ! ReaL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! ReaL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/yCRAY90.F agcm5.3-g77/src/sysdep/yCRAY90.F *** agcm5.3/src/sysdep/yCRAY90.F Wed Sep 24 05:29:11 1997 --- agcm5.3-g77/src/sysdep/yCRAY90.F Sun Jul 4 21:51:49 1999 *************** *** 58,67 **** SAVE * * [INTERNAL WORK] ! REAL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! REAL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = SECOND() CPUTIM = TICKS - TICK0 --- 58,67 ---- SAVE * * [INTERNAL WORK] ! ReaL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! ReaL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = SECOND() CPUTIM = TICKS - TICK0 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/yDECA90.F agcm5.3-g77/src/sysdep/yDECA90.F *** agcm5.3/src/sysdep/yDECA90.F Wed Sep 24 05:31:51 1997 --- agcm5.3-g77/src/sysdep/yDECA90.F Sun Jul 4 21:51:49 1999 *************** *** 60,69 **** SAVE * * [INTERNAL WORK] ! REAL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! REAL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 --- 60,69 ---- SAVE * * [INTERNAL WORK] ! ReaL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! ReaL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/yLinux.F agcm5.3-g77/src/sysdep/yLinux.F *** agcm5.3/src/sysdep/yLinux.F Thu Jan 1 09:00:00 1970 --- agcm5.3-g77/src/sysdep/yLinux.F Sun Jul 4 21:51:49 1999 *************** *** 0 **** --- 1,223 ---- + * PACKAGE YSUN !" 機種依存ルーチン for Sun Fortran + * + *" [HIS] 90/05/19(numaguti) + *" 92/06/05(takepiro) + *" 97/09/23(horinout) NFILN=38->100 + * + ********************************************************************* + SUBROUTINE ERRTRA !" エラートレースバック + * + * [INTERNAL WORK] + INTEGER IPID, IRET + * + * [INTRINSIC FUNC] + INTEGER GETPID + INTEGER KILL + * + IPID=GETPID() + IRET=KILL(IPID,5) + * + RETURN + END + ********************************************************************* + SUBROUTINE YCLOCP !" CPU時間を出力 + I ( HREM ) + * + * [INPUT] + CHARACTER HREM*(*) !" 表題 + * + * [INTERNAL SAVE] + REAL CPUOLD !" 前の呼出の経過時間 + REAL VPUOLD !" 前の呼出のユーザ時間 + DATA CPUOLD, VPUOLD / 0. , 0. / + SAVE + * + * [INTERNAL WORK] + REAL CPUTIM !" 経過時間 + REAL VPUTIM !" ユーザ時間 + * + CALL YCLOCK( CPUTIM, VPUTIM ) + WRITE ( 6,9000 ) HREM, CPUTIM - CPUOLD, VPUTIM - VPUOLD, + & CPUTIM , VPUTIM + 9000 FORMAT( ' ',A8,': CPU/VPU TIME = ', 2E10.5, ' : ', 2E10.5 ) + CPUOLD = CPUTIM + VPUOLD = VPUTIM + * + RETURN + END + ******************************************************************** + SUBROUTINE YCLOCK !" CPU時間を与える + O ( CPUTIM, VPUTIM ) + * + * [OUTPUT] + REAL CPUTIM !" 経過時間 + REAL VPUTIM !" ユーザ時間 + * + * [INTERNAL SAVE] + REAL TICKS + REAL TICK0, TUSR0 !" スタート + DATA TICK0, TUSR0 / 0. , 0. / + SAVE + * + * [INTERNAL WORK] + ReaL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 + * + * [INTRINSIC FUNC] + ReaL*4 ETIME !" *4 に変更 by takepiro 92/06/15 + * + TICKS = ETIME( TARRAY ) + CPUTIM = TICKS - TICK0 + VPUTIM = TARRAY( 1 ) - TUSR0 + * + RETURN + *====================================================================== + ENTRY YCLOCL !" CPU時間クリアー + * + TICK0 = ETIME( TARRAY ) + TUSR0 = TARRAY( 1 ) + * + RETURN + END + *********************************************************************** + SUBROUTINE MKFILN !" ファイル名 HCH を HREP に置き換え + M ( HFILE , + I HCH , HREP ) + * + * [MODIFY] + CHARACTER HFILE *(*) + * + * [INPUT] + CHARACTER HCH *1 + CHARACTER HREP *(*) + * + * [INTERNAL WORK] + INTEGER NFILN + PARAMETER (NFILN=100) + CHARACTER HFILX *(NFILN) + INTEGER NFILE, NREP, I, II + * + * [EXTERNAL FUNC] + INTEGER LENC + * + HFILX = HFILE + HFILE = ' ' + NREP = LENC( HREP ) + NFILE = LEN ( HFILE ) + II = 1 + * + DO 1100 I = 1, MIN( NFILE, NFILN ) + IF ( II .GT. NFILE ) GOTO 1200 + IF ( HFILX(I:I) .EQ. HCH ) THEN + HFILE(II:II+NREP-1) = HREP(1:NREP) + II = II + NREP + ELSE + HFILE(II:II) = HFILX(I:I) + II = II + 1 + ENDIF + 1100 CONTINUE + 1200 CONTINUE + * + RETURN + END + *********************************************************************** + SUBROUTINE REWNML !" NAMELISTファイル, 入力巻き戻し + O ( IFILE, JFILE ) + * + * [OUTPUT] + INTEGER IFILE + INTEGER JFILE + * + * [INTERNAL SAVE] + INTEGER IFILEZ, JFILEZ + DATA IFILEZ / 5 / + DATA JFILEZ / 6 / + SAVE + * + REWIND ( IFILEZ, ERR = 1999 ) + IFILE = IFILEZ + JFILE = JFILEZ + RETURN + * + 1999 IF ( IFILEZ .EQ. 5 ) THEN + CALL MSGDMP( 'W','REWNML','UNIT 5 MAY BE A STANDARD INPUT' ) + IFILE = IFILEZ + JFILE = JFILEZ + ELSE + CALL MSGDMP( 'E','REWNML','ERROR IN REWINDING' ) + * STOP + ENDIF + * + RETURN + *====================================================================== + ENTRY SETNML !" NAMELIST入出力ファイルセット + I ( IFILE, JFILE ) + * + IFILEZ = IFILE + JFILEZ = JFILE + * + RETURN + END + *********************************************************************** + SUBROUTINE YPREP !" システム前処理 + * + RETURN + END + *********************************************************************** + SUBROUTINE YFINE !" システム後処理 + * + RETURN + END + *********************************************************************** + INTEGER FUNCTION IOSLEV !" 入出力エラーレベル + I ( IOS ) + * + * [INPUT] + INTEGER IOS + * + IF ( IOS.EQ.0 ) THEN + IOSLEV = 0 + ELSE + IOSLEV = 2 + ENDIF + * + RETURN + END + *********************************************************************** + SUBROUTINE YDATE !" 現在日付(yyyy mm dd)の取得 + O ( HDATE ) + * + * [OUTPUT] + CHARACTER HDATE *(*) !" 日付(yyyy/mm/dd) + * + * [INTERNAL WORK] + INTEGER IDATE1 ( 3 ) !" 日付(yyyy, mm, dd) + * + CALL IDATE ( IDATE1 ) + IF ( IDATE1(1) .LT. 100 ) IDATE1(1) = IDATE1(1) + 1900 + + HDATE = '****/0*/0*' + WRITE ( HDATE(1:4) , '(I4)' ) IDATE1(1) + WRITE ( HDATE(6:7) , '(I2.2)' ) IDATE1(2) + WRITE ( HDATE(9:10), '(I2.2)' ) IDATE1(3) + * + RETURN + END + *********************************************************************** + SUBROUTINE YTIME !" 現在時刻(hh mm ss)の取得 + O ( HTIME ) + * + * [OUTPUT] + CHARACTER HTIME *(*) !" 時刻(hh:mm:ss) + * + * [INTERNAL WORK] + INTEGER ITIME1 ( 3 ) !" 時刻(hh, mm, ss) + * + CALL ITIME ( ITIME1 ) + + HTIME = '0*:0*:0*' + WRITE ( HTIME(1:2), '(I2.2)' ) ITIME1(1) + WRITE ( HTIME(4:5), '(I2.2)' ) ITIME1(2) + WRITE ( HTIME(7:8), '(I2.2)' ) ITIME1(3) + * + RETURN + END diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/ySX3.F agcm5.3-g77/src/sysdep/ySX3.F *** agcm5.3/src/sysdep/ySX3.F Wed Sep 24 05:32:02 1997 --- agcm5.3-g77/src/sysdep/ySX3.F Sun Jul 4 21:51:49 1999 *************** *** 44,50 **** REAL VPUTIM !" ベクトル時間 * * [INTERNAL SAVE] ! REAL * 8 CPUT0, CPUT1 DATA CPUT0 / 0. / SAVE CPUT0 * --- 44,50 ---- REAL VPUTIM !" ベクトル時間 * * [INTERNAL SAVE] ! ReaL * 8 CPUT0, CPUT1 DATA CPUT0 / 0. / SAVE CPUT0 * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/sysdep/ySun.F agcm5.3-g77/src/sysdep/ySun.F *** agcm5.3/src/sysdep/ySun.F Wed Sep 24 05:31:51 1997 --- agcm5.3-g77/src/sysdep/ySun.F Sun Jul 4 21:51:49 1999 *************** *** 60,69 **** SAVE * * [INTERNAL WORK] ! REAL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! REAL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 --- 60,69 ---- SAVE * * [INTERNAL WORK] ! ReaL*4 TARRAY( 2 ) !" *4 に変更 by takepiro 92/06/15 * * [INTRINSIC FUNC] ! ReaL*4 ETIME !" *4 に変更 by takepiro 92/06/15 * TICKS = ETIME( TARRAY ) CPUTIM = TICKS - TICK0 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/usr/horinout/prep/ltinit.F agcm5.3-g77/src/usr/horinout/prep/ltinit.F *** agcm5.3/src/usr/horinout/prep/ltinit.F Mon Feb 17 19:28:06 1997 --- agcm5.3-g77/src/usr/horinout/prep/ltinit.F Sun Jul 4 21:51:50 1999 *************** *** 146,153 **** * * [INTERNAL WORK] * ! REAL * 8 QSINLA( JDIM ) !" sin（緯度）：倍精度 ! REAL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 * *" < 2. 球面調和関数展開パラメーター > * --- 146,153 ---- * * [INTERNAL WORK] * ! ReaL * 8 QSINLA( JDIM ) !" sin（緯度）：倍精度 ! ReaL * 8 QGW ( JDIM ) !" ガウス荷重 ：倍精度 * *" < 2. 球面調和関数展開パラメーター > * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/util/ufftr.F agcm5.3-g77/src/util/ufftr.F *** agcm5.3/src/util/ufftr.F Sun Apr 10 20:01:37 1994 --- agcm5.3-g77/src/util/ufftr.F Sun Jul 4 22:09:44 1999 *************** *** 142,148 **** MM = 0 1100 CONTINUE IF ( L .LT. N ) THEN ! THETA = 2. * PI / REAL( 2*L ) * DO 1110 M = 0, L-1 MM = MM + 1 --- 142,148 ---- MM = 0 1100 CONTINUE IF ( L .LT. N ) THEN ! THETA = 2. * PI / Real( 2*L ) * DO 1110 M = 0, L-1 MM = MM + 1 diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/util/usphe.F agcm5.3-g77/src/util/usphe.F *** agcm5.3/src/util/usphe.F Sun Apr 10 20:01:37 1994 --- agcm5.3-g77/src/util/usphe.F Sun Jul 4 22:19:14 1999 *************** *** 627,634 **** MI = 2*MM + 2 DO 2110 K = 1, KMAX DO 2110 J = 1, JDIM ! ZDATA( MR, J, K ) = - REAL( M ) * ZDW ( MI, J, K ) ! ZDATA( MI, J, K ) = REAL( M ) * ZDW ( MR, J, K ) 2110 CONTINUE 2100 CONTINUE * --- 627,634 ---- MI = 2*MM + 2 DO 2110 K = 1, KMAX DO 2110 J = 1, JDIM ! ZDATA( MR, J, K ) = - Real( M ) * ZDW ( MI, J, K ) ! ZDATA( MI, J, K ) = Real( M ) * ZDW ( MR, J, K ) 2110 CONTINUE 2100 CONTINUE * diff --exclude=tags --exclude=TAGS* --exclude=*.o --exclude=*.f -rcN agcm5.3/src/util/uspst.F agcm5.3-g77/src/util/uspst.F *** agcm5.3/src/util/uspst.F Sun Apr 10 20:01:37 1994 --- agcm5.3-g77/src/util/uspst.F Sun Jul 4 21:51:49 1999 *************** *** 25,35 **** INTEGER NMO ( 2, 0:MMAX, 0:LMAX ) !" スペクトルの添字順番 * * [INPUT] ! REAL * 8 QSINLA( JMAX ) !" sin（緯度）：倍精度 * * [WORK] ! REAL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! REAL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * * [INTERNAL WORK] INTEGER J, L, M, LEND, MEND, NMH --- 25,35 ---- INTEGER NMO ( 2, 0:MMAX, 0:LMAX ) !" スペクトルの添字順番 * * [INPUT] ! ReaL * 8 QSINLA( JMAX ) !" sin（緯度）：倍精度 * * [WORK] ! ReaL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! ReaL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * * [INTERNAL WORK] INTEGER J, L, M, LEND, MEND, NMH *************** *** 88,102 **** INTEGER MMAX * * [INPUT] ! REAL * 8 CTHETA !" μ=sin(緯度) * * [OUTPUT] ! REAL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! REAL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * * [INTERNAL WORK] INTEGER M, N ! REAL * 8 STHETA, EPS, EPSM, EPSP * * *" < 1. ( １−μ＾2 )^1/2 > --- 88,102 ---- INTEGER MMAX * * [INPUT] ! ReaL * 8 CTHETA !" μ=sin(緯度) * * [OUTPUT] ! ReaL * 8 QPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm ルジャンドル ! ReaL * 8 QDPNM ( 0:NMAX+1, 0:MMAX ) !" Ｐnm μ微分 * * [INTERNAL WORK] INTEGER M, N ! ReaL * 8 STHETA, EPS, EPSM, EPSP * * *" < 1. ( １−μ＾2 )^1/2 > *************** *** 166,178 **** INTEGER JMAX * * [OUTPUT] ! REAL * 8 CTHETA( JMAX ) !" sin(緯度)＝μ ! REAL * 8 GW ( JMAX ) !" ガウス荷重 * * [INTERNAL PARM] ! REAL * 8 DELTP ! REAL * 8 DELTPD ! REAL * 8 DELTPS DATA DELTPD / 1.0 D-28 / !' ４倍精度の収束条件 DATA DELTPS / 1.0 D-14 / !' 倍精度の収束条件 INTEGER ITRMAX --- 166,178 ---- INTEGER JMAX * * [OUTPUT] ! ReaL * 8 CTHETA( JMAX ) !" sin(緯度)＝μ ! ReaL * 8 GW ( JMAX ) !" ガウス荷重 * * [INTERNAL PARM] ! ReaL * 8 DELTP ! ReaL * 8 DELTPD ! ReaL * 8 DELTPS DATA DELTPD / 1.0 D-28 / !' ４倍精度の収束条件 DATA DELTPS / 1.0 D-14 / !' 倍精度の収束条件 INTEGER ITRMAX *************** *** 181,190 **** * [INTERNAL WORK] INTEGER JMAXD PARAMETER (JMAXD=200) ! REAL * 8 QPN( 0:JMAXD+1 ) !" Ｐn (倍精度) INTEGER N, J, ITER ! REAL * 8 PI, X0, DELTX, R, RP ! REAL * 8 EPS, EPSM, EPSP, QDPN * IF ( JMAX .GT. JMAXD) THEN CALL MSGDMP( 'E','GAUSS','INSUFFICIENT WORK STRAGE.') --- 181,190 ---- * [INTERNAL WORK] INTEGER JMAXD PARAMETER (JMAXD=200) ! ReaL * 8 QPN( 0:JMAXD+1 ) !" Ｐn (倍精度) INTEGER N, J, ITER ! ReaL * 8 PI, X0, DELTX, R, RP ! ReaL * 8 EPS, EPSM, EPSP, QDPN * IF ( JMAX .GT. JMAXD) THEN CALL MSGDMP( 'E','GAUSS','INSUFFICIENT WORK STRAGE.')