Changeset 3132 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Nov 21, 2023, 8:25:02 AM (20 months ago)

Author:

emillour

Message:

Mars PCM:
Bug fix in phyetat0, igcm_co2_tmp was not initialed when paleoclimate==.false.
but used nonetheless. Added an if(paleoclimate) around computation of
perennial_co2ice() as it should only be computed in that context; also added
an if(paleoclimate) around writting perennial_co2ice() to restartfi.nc in
physdem1 for the same reason.
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• phyetat0_mod.F90 (modified) (3 diffs)
• phyredem.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -781 +781 @@
 ! Perennial CO2 ice
 perennial_co2ice(:,:) = 0.
-if (startphy_file) then
+if (paleoclimate) then
+  if (startphy_file) then
     call get_field("perennial_co2ice",perennial_co2ice,found,indextime)
     if (.not. found) then
@@ -793 +794 @@
         endif
     endif ! not found
-else
+  else
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
         do islope = 1,nslope
@@ -800 +801 @@
         enddo
     endif
-endif !startphy_file
+  endif !startphy_file
+endif ! of if (paleoclimate)
 
 ! close file:

trunk/LMDZ.MARS/libf/phymars/phyredem.F90

@@ -227 +227 @@
  
   ! Perennial CO2 ice layer
-  call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
+  if(paleoclimate) then
+    call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
+  endif
 
   ! Surface temperature