Changeset 1673 for LMDZ5/trunk/libf/dyn3dmem/etat0_netcdf.F90
- Timestamp:
- Oct 27, 2012, 4:23:07 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3dmem/etat0_netcdf.F90
r1658 r1673 1 1 ! 2 ! $Id : etat0_netcdf.F90 1425 2010-09-02 13:45:23Z lguez$2 ! $Id$ 3 3 ! 4 4 !------------------------------------------------------------------------------- 5 5 ! 6 SUBROUTINE etat0_netcdf(ib, masque, letat0)6 SUBROUTINE etat0_netcdf(ib, masque, phis, letat0) 7 7 ! 8 8 !------------------------------------------------------------------------------- … … 37 37 LOGICAL, INTENT(IN) :: ib ! barycentric interpolat. 38 38 REAL, DIMENSION(iip1,jjp1), INTENT(INOUT) :: masque ! land mask 39 REAL, DIMENSION(iip1,jjp1), INTENT(OUT) :: phis ! geopotentiel au sol 39 40 LOGICAL, INTENT(IN) :: letat0 ! F: masque only required 40 41 #ifndef CPP_EARTH … … 51 52 REAL, DIMENSION(klon) :: tsol, qsol 52 53 REAL, DIMENSION(klon) :: sn, rugmer, run_off_lic_0 53 REAL, DIMENSION(iip1,jjp1) :: orog, rugo, psol , phis54 REAL, DIMENSION(iip1,jjp1) :: orog, rugo, psol 54 55 REAL, DIMENSION(iip1,jjp1,llm+1) :: p3d 55 56 REAL, DIMENSION(iip1,jjp1,llm) :: uvent, t3d, tpot, qsat, qd … … 98 99 REAL :: dummy 99 100 LOGICAL :: ok_newmicro, ok_journe, ok_mensuel, ok_instan, ok_hf 100 LOGICAL :: ok_LES, ok_ade, ok_aie, aerosol_couple, new_aod 101 LOGICAL :: ok_LES, ok_ade, ok_aie, aerosol_couple, new_aod, callstats 101 102 INTEGER :: iflag_radia, flag_aerosol 102 103 REAL :: bl95_b0, bl95_b1, fact_cldcon, facttemps, ratqsbas, ratqshaut … … 130 131 !--- CONSTRUCT A GRID 131 132 CALL conf_phys( ok_journe, ok_mensuel, ok_instan, ok_hf, ok_LES, & 133 callstats, & 132 134 solarlong0,seuil_inversion, & 133 135 fact_cldcon, facttemps,ok_newmicro,iflag_radia, & … … 137 139 flag_aerosol, new_aod, & 138 140 bl95_b0, bl95_b1, & 139 iflag_thermals,nsplit_thermals,tau_thermals, & 140 iflag_thermals_ed,iflag_thermals_optflux, & 141 iflag_coupl,iflag_clos,iflag_wake, read_climoz, & 141 read_climoz, & 142 142 alp_offset) 143 143 … … 251 251 !******************************************************************************* 252 252 CALL pression(ip1jmp1, ap, bp, psol, p3d) 253 CALL exner_hyb(ip1jmp1, psol, p3d, alpha, beta, pks, pk, y) 253 if (pressure_exner) then 254 CALL exner_hyb(ip1jmp1, psol, p3d, alpha, beta, pks, pk, y) 255 else 256 CALL exner_milieu(ip1jmp1,psol,p3d,beta,pks,pk,y) 257 endif 254 258 pls(:,:,:)=preff*(pk(:,:,:)/cpp)**(1./kappa) 255 259 ! WRITE(lunout,*) 'P3D :', p3d(10,20,:)
Note: See TracChangeset
for help on using the changeset viewer.