- Timestamp:
- 07/16/19 01:45:04 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dcmip/dcmip_initial_conditions_test_1_2_3_v5.f90
r531 r957 1230 1230 real(rstd), parameter :: & 1231 1231 Om = 0.d0, & ! Rotation Rate of Earth 1232 u0 = 20.d0, & ! Reference Velocity1233 !u0 = 0.d0, & ! FIXME : no zonal wind for NH tests1232 ! u0 = 20.d0, & ! Reference Velocity 1233 u0 = 0.d0, & ! FIXME : no zonal wind for NH tests 1234 1234 Teq = 300.d0, & ! Temperature at Equator 1235 1235 Peq = 100000.d0, & ! Reference PS at Equator 1236 1236 ztop = 10000.d0, & ! Model Top 1237 1237 lambdac = 2.d0*pi/3.d0, & ! Lon of Pert Center 1238 d = 5000.d0, & ! Width for Pert1238 d = 625000.d0, & ! Width for Pert for X=1 1239 1239 phic = 0.d0, & ! Lat of Pert Center 1240 1240 delta_theta = 1.d0, & ! Max Amplitude of Pert … … 1333 1333 r = as * ACOS (sin_tmp + cos_tmp*cos(lon-lambdac)) 1334 1334 1335 s = ( d**2)/(d**2 + r**2)1335 s = ((d/X)**2)/((d/X)**2 + r**2) 1336 1336 1337 1337 theta_pert = delta_theta*s*sin(2.d0*pi*height/Lz)
Note: See TracChangeset
for help on using the changeset viewer.