Changeset 3025 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Aug 7, 2023, 6:02:51 AM (17 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/libf/dynphy_lonlat/phymars

Files:

• lect_start_archive.F (modified) (3 diffs)
• newstart.F (modified) (2 diffs)

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