Changeset 217 for codes/icosagcm
- Timestamp:
- 07/16/14 10:08:16 (10 years ago)
- Location:
- codes/icosagcm/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/param_sets/make_rundefs.sh
r215 r217 55 55 function itau_adv { print 'Advection called every itau_adv time steps' integer 2 itau_adv $1 56 56 } 57 function nqtot { print 'Number of tracers' integer 1 nqtot $1 57 function nqtot { print 'Number of tracers' integer 1 nqtot $1 58 58 } 59 59 function caldyn { print 'Equations solved' '[gcm|adv]' gcm caldyn $1 … … 84 84 function etat0 { 85 85 OPTS=$(printf '\n# [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear|dcmip3]') 86 print 'Initial state' "$OPTS" jablonowsky06 etat0 $1 86 print 'Initial state' "$OPTS" jablonowsky06 etat0 $1 87 87 } 88 88 function dissip … … 98 98 } 99 99 100 function physics { print 'Physics package' '[ automatic|dcmip]' automaticphysics $1100 function physics { print 'Physics package' '[none|held_suarez|dcmip]' none physics $1 101 101 } 102 102 … … 160 160 group 'Time' time_style lmd ndays 1200 write_period 86400 161 161 group 'Planet' radius 6.371e6 g 9.8 omega 7.292e-5 kappa 0.2857143 cpp 1004 preff 1e5 162 group 'Physical parameters' etat0 held_suarez nqtot 1162 group 'Physical parameters' etat0 held_suarez physics held_suarez 163 163 dissip 50000 2 50000 2 50000 2 164 164 } … … 186 186 dissip 1800 1 1800 2 1800 2 187 187 group 'Time' run_length 3.6e6 write_period 50000 188 group 'Physical parameters' scale_factor 500 guided_type dcmip1 nqtot 1etat0 dcmip2_schaer_noshear188 group 'Physical parameters' scale_factor 500 guided_type dcmip1 etat0 dcmip2_schaer_noshear 189 189 } 190 190 … … 196 196 dissip 1800 1 1800 2 1800 2 197 197 group 'Time' run_length 4.5e5 write_period 12500 198 group 'Physical parameters' omega 0 scale_factor 125 nqtot 1etat0 dcmip3198 group 'Physical parameters' omega 0 scale_factor 125 etat0 dcmip3 199 199 } 200 200 -
codes/icosagcm/trunk/src/physics.f90
r215 r217 11 11 TYPE(t_field),POINTER :: f_dulon(:), f_dulat(:) 12 12 13 CHARACTER(LEN=255) :: physics_type ="automatic"13 CHARACTER(LEN=255) :: physics_type 14 14 !$OMP THREADPRIVATE(physics_type) 15 15 … … 26 26 IMPLICIT NONE 27 27 28 CALL allocate_field(f_dulon,field_t,type_real,llm, name='dulon') 29 CALL allocate_field(f_dulat,field_t,type_real,llm, name='dulat') 30 CALL init_pack_before ! Compute physics_inout%ngrid and offsets used by pack/unpack 31 32 physics_type='automatic' 28 physics_inout%dt_phys = dt*itau_physics 29 physics_type='none' 33 30 CALL getin("physics",physics_type) 34 35 31 SELECT CASE(TRIM(physics_type)) 36 CASE ('automatic') 37 etat0_type='jablonowsky06' 38 CALL getin("etat0",etat0_type) 39 SELECT CASE(TRIM(etat0_type)) 40 CASE('held_suarez') 41 phys_type = phys_HS94 42 CASE DEFAULT 43 IF(is_mpi_root) PRINT*,"NO PHYSICAL PACKAGE USED" 44 phys_type = phys_none 45 END SELECT 46 32 CASE ('none') 33 IF(is_mpi_root) PRINT*,"NO PHYSICAL PACKAGE USED" 34 phys_type = phys_none 35 CASE ('held_suarez') 36 phys_type = phys_HS94 47 37 CASE ('dcmip') 38 CALL allocate_field(f_dulon,field_t,type_real,llm, name='dulon') 39 CALL allocate_field(f_dulat,field_t,type_real,llm, name='dulat') 40 CALL init_pack_before ! Compute physics_inout%ngrid and offsets used by pack/unpack 48 41 CALL init_physics_dcmip 42 CALL init_pack_after ! Defines Ai, lon, lat in physics_inout 49 43 phys_type = phys_DCMIP 50 44 CASE DEFAULT 51 45 IF(is_mpi_root) PRINT*, 'init_physics : Bad selector for variable physics <',& 52 TRIM(physics_type), '> options are < automatic>, <dcmip>, <dry>'46 TRIM(physics_type), '> options are <none>, <held_suarez>, <dcmip>' 53 47 STOP 54 48 END SELECT 55 49 56 50 IF(is_mpi_root) PRINT *, 'phys_type = ',phys_type 57 physics_inout%dt_phys = dt*itau_physics58 CALL init_pack_after ! Defines Ai, lon, lat in physics_inout59 51 END SUBROUTINE init_physics 60 52
Note: See TracChangeset
for help on using the changeset viewer.