Changeset 3025 for trunk

Timestamp:

Aug 7, 2023, 6:02:51 AM (16 months ago)

Author:

emillour

Message:

Mars PCM:
Fixed a bug in newstart call to lect_start_archive; missing perenial_co2ice
argument. Turned lect_start_archive into a module so this cant't happen again.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/dynphy_lonlat/phymars/lect_start_archive.F (modified) (3 diffs)
• libf/dynphy_lonlat/phymars/newstart.F (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4162 +4162 @@
 Complete rewriting (much simpler) of script "run_pem1d_1" in deftank/pem/ to launch chained simulations of 1D PCM/PEM. To be adapted to 3D in the future.
 Addition of a bash script "modify_startfi.sh" to modify orbital parameters in "startfi.nc". It is needed to initialize the PCM/PEM simulation with right orbital parameters according to the year in Laskar data.
+
+== 07/08/2023 == EM
+Fixed a bug in newstart call to lect_start_archive; missing perenial_co2ice
+argument. Turned lect_start_archive into a module so this cant't happen again.

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F

@@ +1 @@
+      MODULE lect_start_archive_mod
+      
+      IMPLICIT NONE
+      
+      CONTAINS
+      
       SUBROUTINE lect_start_archive(ngrid,nlayer,nqtot,
      &     date,tsurf,tsoil,inertiesoil,albedo,emis,q2,
@@ -1277 +1283 @@
      &       watercap(:,1))
 
-c Watercap
+c Perenial CO2 ice cap
       call interp_horiz (peren_co2iceold(:,:,1),peren_co2iceS(:,:,1),
      &       imold,jmold,iim,jjm,1,rlonuold,rlatvold,rlonu,rlatv)
@@ -1686 +1692 @@
 
 !      write(*,*)'lect_start_archive: END'
-      return
-      end
+
+      END SUBROUTINE lect_start_archive
+
+      END MODULE lect_start_archive_mod

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -37 +37 @@
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
       use geometry_mod, only: longitude,latitude,cell_area
+      use lect_start_archive_mod, only: lect_start_archive
       use phyetat0_mod, only: phyetat0
       use phyredem, only: physdem0, physdem1
@@ -444 +445 @@
      &   date,tsurf,tsoil,inertiesoil,albedo,emis,q2,
      &   t,ucov,vcov,ps,teta,phisold_newgrid,q,qsurf,
-     &   tauscaling,totcloudfrac,surfith,nid,watercap)
+     &   tauscaling,totcloudfrac,surfith,nid,watercap,perenial_co2ice)
         write(*,*) "OK, read start_archive file"
         ! copy soil thermal inertia