Changeset 549
- Timestamp:
- 09/18/17 13:55:05 (7 years ago)
- Location:
- codes/icosagcm/devel/src/initial
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/initial/etat0.f90
r531 r549 16 16 USE etat0_database_mod, ONLY: init_etat0_database => init_etat0 17 17 USE etat0_start_file_mod, ONLY: init_etat0_start_file => init_etat0 18 USE etat0_heldsz_mod, ONLY: init_etat0_held_suarez => init_etat0 18 19 IMPLICIT NONE 19 20 … … 32 33 CASE ('academic') 33 34 CASE ('held_suarez') 35 CALL init_etat0_held_suarez 34 36 CASE ('venus') 35 37 CASE ('dcmip1') -
codes/icosagcm/devel/src/initial/etat0_heldsz.f90
r531 r549 3 3 IMPLICIT NONE 4 4 PRIVATE 5 SAVE 5 6 6 7 TYPE(t_field),POINTER :: f_theta_eq(:) … … 12 13 !$OMP THREADPRIVATE(done) 13 14 14 REAL(rstd),SAVE :: teta0,ttp,delt_y,delt_z,eps15 !$OMP THREADPRIVATE( teta0,ttp,delt_y,delt_z,eps)15 REAL(rstd),SAVE :: p0,teta0,ttp,delt_y,delt_z,eps 16 !$OMP THREADPRIVATE(p0,teta0,ttp,delt_y,delt_z,eps) 16 17 17 18 REAL(rstd),SAVE :: knewt_g, k_f,k_c_a,k_c_s 18 19 !$OMP THREADPRIVATE(knewt_g, k_f,k_c_a,k_c_s) 19 20 20 PUBLIC :: etat0, held_suarez21 PUBLIC :: init_etat0, etat0, held_suarez 21 22 22 23 CONTAINS … … 47 48 END SUBROUTINE test_etat0_heldsz 48 49 50 SUBROUTINE init_etat0 51 p0=1e5 ! default value of initial surface pressure as in H&S paper 52 ! p0=101080 for CMIP5 aquaplanets, cf LMDZ5 ini_aqua 53 CALL getin('heldsz_p0',p0) 54 END SUBROUTINE init_etat0 55 49 56 SUBROUTINE etat0(f_ps,f_phis,f_theta_rhodz,f_u, f_q) 50 57 USE theta2theta_rhodz_mod … … 78 85 phis=f_phis(ind) 79 86 u=f_u(ind) 80 ps(:)= 1e587 ps(:)=p0 81 88 phis(:)=0. 82 89 u(:,:)=0
Note: See TracChangeset
for help on using the changeset viewer.