- Timestamp:
- 09/25/19 16:40:08 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dcmip/dcmip_initial_conditions_test_1_2_3_v5.f90
r960 r976 1026 1026 1027 1027 SUBROUTINE test2_schaer_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,shear,u,v,w,t,phis,ps,rho,q) 1028 1028 USE getin_mod, ONLY : getin 1029 1029 IMPLICIT NONE 1030 1030 !----------------------------------------------------------------------- … … 1074 1074 as = a/X, & ! New Radius of small Earth 1075 1075 ueq = 20.d0, & ! Reference Velocity 1076 Teq = 300.d0, & ! Temperature at Equator1077 1076 Peq = 100000.d0, & ! Reference PS at Equator 1078 1077 ztop = 30000.d0, & ! Model Top … … 1090 1089 real(rstd) :: c ! Shear 1091 1090 1091 REAL(rstd) :: Teq = 300.d0 ! default Temperature at Equator, can be modified in run.def 1092 1093 CALL getin("dcmip2_Teq",Teq) 1092 1094 !----------------------------------------------------------------------- 1093 1095 ! PHIS (surface geopotential) … … 1197 1199 !========== 1198 1200 SUBROUTINE test3_gravity_wave (X,lon,lat,p,z,zcoords,u,v,w,t,phis,ps,rho,q) 1199 1201 USE getin_mod, ONLY : getin 1200 1202 IMPLICIT NONE 1201 1203 !----------------------------------------------------------------------- … … 1232 1234 u0 = 20.d0, & ! Reference Velocity 1233 1235 ! u0 = 0.d0, & ! FIXME : no zonal wind for NH tests 1234 Teq = 300.d0, & ! Temperature at Equator1235 1236 Peq = 100000.d0, & ! Reference PS at Equator 1236 1237 ztop = 10000.d0, & ! Model Top … … 1253 1254 real(rstd) :: theta_pert ! Pot-temp perturbation 1254 1255 1256 REAL(rstd) :: Teq = 300.d0 ! default Temperature at Equator, can be modified in run.def 1257 1258 CALL getin("dcmip2_Teq",Teq) 1259 1255 1260 as = a/X 1256 1261
Note: See TracChangeset
for help on using the changeset viewer.