Changeset 3114 for trunk/LMDZ.COMMON/libf
- Timestamp:
- Nov 3, 2023, 5:47:39 PM (14 months ago)
- Location:
- trunk/LMDZ.COMMON/libf/evolution
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/changelog.txt
r3108 r3114 127 127 == 27/10/2023 == JBC 128 128 Addition of the possibility to ask for a number of Mars years to be simulated in total (previously it was only in Earth years). 129 130 == 03/11/2023 == JBC 131 Following r3113, addition of 'nqsoil' and 'qsoil' in the arguments of the subroutines 'phyetat0' and 'physdem1' to be able to compile. -
trunk/LMDZ.COMMON/libf/evolution/pem.F90
r3096 r3114 73 73 74 74 #ifndef CPP_STD 75 use comsoil_h, only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil 75 use comsoil_h, only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil, nqsoil, qsoil 76 76 use surfdat_h, only: tsurf, emis, qsurf, watercap, ini_surfdat_h, & 77 77 albedodat, zmea, zstd, zsig, zgam, zthe, & … … 367 367 ! First we read the initial state (starfi.nc) 368 368 #ifndef CPP_STD 369 call phyetat0(FILE_NAME,0,0,nsoilmx,ngrid,nlayer,nq, day_ini,time_phys,tsurf, &370 tsoil,albedo,emis,q2,qsurf, tauscaling,totcloudfrac,wstar,&369 call phyetat0(FILE_NAME,0,0,nsoilmx,ngrid,nlayer,nq,nqsoil,day_ini,time_phys,tsurf, & 370 tsoil,albedo,emis,q2,qsurf,qsoil,tauscaling,totcloudfrac,wstar, & 371 371 watercap,perenial_co2ice,def_slope,def_slope_mean,subslope_dist) 372 372 … … 383 383 allocate(qsurf_read_generic(ngrid,nq)) 384 384 allocate(tsoil_read_generic(ngrid,nsoilmx)) 385 allocate(qsoil_read_generic(ngrid,nsoilmx,nqsoil,nslope)) 385 386 allocate(emis_read_generic(ngrid)) 386 387 allocate(tsurf(ngrid,1)) … … 393 394 allocate(albedo(ngrid,2,1)) 394 395 allocate(inertiesoil(ngrid,nsoilmx,1)) 395 call phyetat0(.true.,ngrid,nlayer,FILE_NAME,0,0,nsoilmx,nq, day_ini,time_phys, &396 tsurf_read_generic,tsoil_read_generic,emis_read_generic,q2, &397 qsurf_read_generic, cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic,&398 tslab,tsea_ice,sea_ice)396 call phyetat0(.true.,ngrid,nlayer,FILE_NAME,0,0,nsoilmx,nq,nqsoil,day_ini,time_phys, & 397 tsurf_read_generic,tsoil_read_generic,emis_read_generic,q2, & 398 qsurf_read_generic,qsoil_read_generic,cloudfrac,totcloudfrac,hice, & 399 rnat,pctsrf_sic,tslab,tsea_ice,sea_ice) 399 400 call surfini(ngrid,nq,qsurf_read_generic,albedo_read_generic,albedo_bareground,albedo_snow_SPECTV,albedo_co2_ice_SPECTV) 400 401 … … 1112 1113 inertiedat,def_slope,subslope_dist) 1113 1114 1114 call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq, 1115 ptimestep,ztime_fin,tsurf,tsoil,inertiesoil, &1116 albedo,emis,q2,qsurf, tauscaling,totcloudfrac, &1115 call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,nqsoil, & 1116 ptimestep,ztime_fin,tsurf,tsoil,inertiesoil, & 1117 albedo,emis,q2,qsurf,qsoil,tauscaling,totcloudfrac, & 1117 1118 wstar,watercap,perenial_co2ice) 1118 1119 #else … … 1121 1122 albedo_bareground,inertiedat,zmea,zstd,zsig,zgam,zthe) 1122 1123 1123 call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq, 1124 ptimestep,ztime_fin,tsurf,tsoil,emis,q2,qsurf, &1125 cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic, &1126 ts lab,tsea_ice,sea_ice)1124 call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,nqsoil, & 1125 ptimestep,ztime_fin,tsurf,tsoil,emis,q2,qsurf,qsoil, & 1126 cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic,tslab, & 1127 tsea_ice,sea_ice) 1127 1128 #endif 1128 1129 write(*,*) "restartfi_evol.nc has been written"
Note: See TracChangeset
for help on using the changeset viewer.