Changeset 3203 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Feb 7, 2024, 12:14:38 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:

• Addition of the possibility to use subslopes parametization in 1D. To do so, put 'nslope=x' in "run.def" where x (1, 3, 5 or 7) is the number of slopes you want to, or modify it in an appropriate "startfi.nc". Then, default subslopes definition and distribution are used by the model. The already existing case of using one slope whose inclination and orientation are defined in "run.def" is still possible.
• Some cleanings throughout the code, in particular related to unused variables.

JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• orbit_param_criterion_mod.F90 (modified) (1 diff)
• pem.F90 (modified) (1 diff)
• pemetat0.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/orbit_param_criterion_mod.F90

@@ -77 +77 @@
 do ilask = 1,size(yearlask,1)
     read(73,*) yearlask(ilask), obllask(ilask), ecclask(ilask), lsplask(ilask)
-    yearlask(ilask) = yearlask(ilask)*10./convert_years ! Conversion from Laskar's kilo Earth years to PEM Martian years
+    yearlask(ilask) = yearlask(ilask)*1000./convert_years ! Conversion from Laskar's kilo Earth years to PEM Martian years
     if (yearlask(ilask) > Year) last_ilask = ilask + 1
 enddo

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -320 +320 @@
     endif
 
-    call init_testphys1d('start1D_evol.txt','startfi_evol.nc',.true.,therestart1D,therestartfi,ngrid,nlayer,610., &
-                         nq,q,time_0,ps(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,  &
+    call init_testphys1d('start1D_evol.txt','startfi_evol.nc',therestart1D,therestartfi,ngrid,nlayer,610.,nq,q, &
+                         time_0,ps(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,     &
                          play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,CO2cond_ps)
     ps(2) = ps(1)

trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

@@ -253 +253 @@
                 tsoil_tmp_yr2(:,:,islope) = tsoil_tmp_yr1(:,:,islope)
                 call compute_tsoil_pem(ngrid,nsoil_PEM,.false.,TI_PEM(:,:,islope),timestep,tsurf_avg_yr2(:,islope),tsoil_tmp_yr2(:,:,islope))
-                
+
                 do iloop = nsoil_PCM+1,nsoil_PEM
                     tsoil_PEM(:,iloop,islope) = tsoil_tmp_yr2(:,iloop,islope)