- Timestamp:
- 12/30/19 12:37:59 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dcmip/dcmip_initial_conditions_test_1_2_3_v5.f90
r976 r995 874 874 ! Test 2-0: Steady-State Atmosphere at Rest in the Presence of Orography 875 875 !========================================================================= 876 SUBROUTINE test2_steady_state_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,u,v,w,t,phis,ps,rho,q )876 SUBROUTINE test2_steady_state_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,u,v,w,t,phis,ps,rho,q,h0_test2) 877 877 878 878 IMPLICIT NONE … … 922 922 lambdam = 3.d0*pi/2.d0, & ! mountain longitude center point (radians) 923 923 phim = 0.d0, & ! mountain latitude center point (radians) 924 h0 = 2000.d0, & ! peak height of the mountain range (m)925 924 Rm = 3.d0*pi/4.d0, & ! mountain radius (radians) 926 925 zetam = pi/16.d0, & ! mountain oscillation half-width (radians) … … 930 929 real(rstd) :: r ! Great circle distance (radians) 931 930 real(rstd) :: zs ! Surface elevation (m) 932 real(rstd) :: exponent ! exponent: g/(Rd * gamma) 933 real(rstd) :: exponent_rev ! reversed exponent 934 931 real(rstd) :: exponent ! exponent: g/(Rd * gamma) 932 real(rstd) :: exponent_rev ! reversed exponent 933 real(rstd) :: h0 ! peak height of the mountain range (m) 934 real(rstd), optional :: h0_test2 935 936 if (present(h0_test2)) then 937 h0 = h0_test2 ! peak height of the mountain range (m), can be modified in run.def 938 else 939 h0 = 250.d0 ! Default value for call without optional variable h0_test2 940 endif 935 941 936 942 !----------------------------------------------------------------------- … … 1025 1031 !===================================================================================== 1026 1032 1027 SUBROUTINE test2_schaer_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,shear,u,v,w,t,phis,ps,rho,q )1028 USE getin_mod, ONLY : getin 1033 SUBROUTINE test2_schaer_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,shear,u,v,w,t,phis,ps,rho,q,h0_test2,Teq_test2) 1034 1029 1035 IMPLICIT NONE 1030 1036 !----------------------------------------------------------------------- … … 1078 1084 lambdac = pi/4.d0, & ! Lon of Schar Mountain Center 1079 1085 phic = 0.d0, & ! Lat of Schar Mountain Center 1080 h0 = 250.d0, & ! Height of Mountain1081 1086 d = 5000.d0, & ! Mountain Half-Width 1082 1087 xi = 4000.d0, & ! Mountain Wavelength … … 1088 1093 real(rstd) :: zs ! Surface height 1089 1094 real(rstd) :: c ! Shear 1090 1091 REAL(rstd) :: Teq = 300.d0 ! default Temperature at Equator, can be modified in run.def 1092 1093 CALL getin("dcmip2_Teq",Teq) 1095 real(rstd) :: h0 ! peak height of the mountain range (m) 1096 real(rstd) :: Teq ! Temperature at Equator 1097 real(rstd), optional :: h0_test2,Teq_test2 ! FIX default value in MODULE etat0_dcmip2_mod 1098 1099 if (present(h0_test2)) then 1100 h0 = h0_test2 ! peak height of the mountain range (m), can be modified in run.def 1101 else 1102 h0 = 250.d0 ! Default value for call without optional variable h0_test2 1103 endif 1104 if (present(Teq_test2)) then 1105 Teq = Teq_test2 ! Temperature at Equator, can be modified in run.def 1106 else 1107 Teq = 300.d0 ! Default value for call without optional variable Teq_test2 1108 endif 1094 1109 !----------------------------------------------------------------------- 1095 1110 ! PHIS (surface geopotential) … … 1199 1214 !========== 1200 1215 SUBROUTINE test3_gravity_wave (X,lon,lat,p,z,zcoords,u,v,w,t,phis,ps,rho,q) 1201 USE getin_mod, ONLY : getin 1216 1202 1217 IMPLICIT NONE 1203 1218 !----------------------------------------------------------------------- … … 1234 1249 u0 = 20.d0, & ! Reference Velocity 1235 1250 ! u0 = 0.d0, & ! FIXME : no zonal wind for NH tests 1251 Teq = 300.d0, & ! Temperature at Equator 1236 1252 Peq = 100000.d0, & ! Reference PS at Equator 1237 1253 ztop = 10000.d0, & ! Model Top … … 1254 1270 real(rstd) :: theta_pert ! Pot-temp perturbation 1255 1271 1256 REAL(rstd) :: Teq = 300.d0 ! default Temperature at Equator, can be modified in run.def1257 1258 CALL getin("dcmip2_Teq",Teq)1259 1260 1272 as = a/X 1261 1273
Note: See TracChangeset
for help on using the changeset viewer.