Changeset 32 for codes/icosagcm/trunk
- Timestamp:
- 07/27/12 16:04:22 (12 years ago)
- Location:
- codes/icosagcm/trunk/src
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/caldyn_gcm.f90
r26 r32 24 24 CALL allocate_caldyn 25 25 CALL getin('write_period',write_period) 26 write_period=write_period/scale_factor 26 27 27 28 itau_out=INT(write_period/dt) -
codes/icosagcm/trunk/src/dissip_gcm.f90
r29 r32 79 79 tau_graddiv(:)=5000 80 80 CALL getin("tau_graddiv",tau) 81 tau_graddiv(:)=tau 81 tau_graddiv(:)=tau/scale_factor 82 82 83 83 CALL getin("nitergdiv",nitergdiv) … … 85 85 tau_gradrot(:)=5000 86 86 CALL getin("tau_gradrot",tau_gradrot) 87 tau_gradrot(:)=tau 87 tau_gradrot(:)=tau/scale_factor 88 88 89 89 CALL getin("nitergrot",nitergrot) … … 92 92 tau_divgrad(:)=5000 93 93 CALL getin("tau_divgrad",tau) 94 tau_divgrad(:)=tau 94 tau_divgrad(:)=tau/scale_factor 95 95 96 96 CALL getin("niterdivgrad",niterdivgrad) -
codes/icosagcm/trunk/src/earth_const.f90
r12 r32 3 3 USE math_const 4 4 5 REAL(rstd),PARAMETER :: radius=6.37122E6 6 ! REAL(rstd),PARAMETER :: radius=1. 7 REAL(rstd),PARAMETER :: g=9.80616 5 REAL(rstd),SAVE :: radius=6.37122E6 6 REAL(rstd),SAVE :: g=9.80616 8 7 REAL(rstd),PARAMETER :: daysec=86400 9 ! REAL(rstd),PARAMETER :: omega=2*Pi/daysec 10 REAL(rstd),PARAMETER :: omega=7.292E-5 11 ! REAL(rstd),PARAMETER :: omega=0 12 REAL(rstd),PARAMETER :: kappa=0.2857143 13 REAL(rstd),PARAMETER :: cpp=1004.70885 14 REAL(rstd),PARAMETER :: preff=101325. 15 REAL(rstd),PARAMETER :: pa=50000. 8 REAL(rstd),SAVE :: omega=7.292E-5 9 REAL(rstd),SAVE :: kappa=0.2857143 10 REAL(rstd),SAVE :: cpp=1004.70885 11 REAL(rstd),SAVE :: preff=101325. 12 REAL(rstd),SAVE :: pa=50000. 13 REAL(rstd),SAVE :: scale_factor=1. 16 14 15 CONTAINS 16 17 SUBROUTINE init_earth_const 18 USE ioipsl 19 IMPLICIT NONE 20 REAL(rstd) :: X=1 21 22 CALL getin("radius",radius) 23 CALL getin("g",g) 24 CALL getin("X",scale_factor) 25 CALL getin("omega",omega) 26 CALL getin("kappa",kappa) 27 CALL getin("cpp",cpp) 28 CALL getin("preff",preff) 29 30 radius=radius/scale_factor 31 omega=omega/scale_factor 32 PRINT *,"radius = ",radius 33 34 END SUBROUTINE init_earth_const 35 36 17 37 END MODULE earth_const 18 38 -
codes/icosagcm/trunk/src/icosa_gcm.f90
r26 r32 18 18 19 19 CALL init_mpipara 20 20 CALL init_earth_const 21 21 CALL init_grid_param 22 22 CALL compute_metric -
codes/icosagcm/trunk/src/timeloop_gcm.f90
r25 r32 44 44 dt=90. 45 45 CALL getin('dt',dt) 46 46 dt=dt/scale_factor 47 47 48 itaumax=100 48 49 CALL getin('itaumax',itaumax) … … 50 51 write_period=0 51 52 CALL getin('write_period',write_period) 53 write_period=write_period/scale_factor 54 52 55 itau_out=INT(write_period/dt) 53 56 … … 93 96 94 97 CALL etat0(f_ps,f_phis,f_theta_rhodz,f_u, f_q) 98 99 CALL writefield('ps',f_ps) 100 CALL writefield('theta_rhodz',f_theta_rhodz) 95 101 96 102 DO it=0,itaumax
Note: See TracChangeset
for help on using the changeset viewer.