Changeset 544 for codes/icosagcm/trunk/src/etat0_heldsz.f90
- Timestamp:
- 08/10/17 16:02:21 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_heldsz.f90
r494 r544 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.