Changeset 2994

Timestamp:

Jul 13, 2023, 3:37:43 PM (17 months ago)

Author:

llange

Message:

MARS PCM
Introduce paleoclimate modul which will contains all the stuff used for
the paleoclimates studies. For now, just the lag layer thicknesses (for
preliminary tests with the 1D model).
LL

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• conf_phys.F (modified) (2 diffs)
• paleoclimate_mod.F90 (added)
• phyetat0_mod.F90 (modified) (2 diffs)
• phyredem.F90 (modified) (3 diffs)
• phys_state_var_init_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -54 +54 @@
       use aeropacity_mod, only: iddist, topdustref
       USE mod_phys_lmdz_transfert_para, ONLY: bcast
+      USE paleoclimate_mod,ONLY: paleoclimate
       IMPLICIT NONE
       include "callkeys.h"
@@ -346 +347 @@
          call getin_p("rayleigh",rayleigh)
          write(*,*)" rayleigh = ",rayleigh
+
+! PALEOCLIMATE
+
+         write(*,*)"Is it paleoclimate run?"
+         paleoclimate=.false. ! default value
+         call getin_p("paleoclimate",paleoclimate)
+         write(*,*)" paleoclimate = ",paleoclimate
 
 

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

@@ -27 +27 @@
   use comsoil_h, only: flux_geo
   USE comslope_mod, ONLY: nslope, major_slope
+  USE paleoclimate_mod, ONLY: paleoclimate, lag_h2o_ice, lag_co2_ice
   implicit none
   
@@ -752 +753 @@
 endif !startphy_file
 
-!
+
+if (paleoclimate) then
+  if (startphy_file) then
+   call get_field("lag_h2o_ice",lag_h2o_ice,found,indextime)
+   write(*,*) 'paleo found start?',found
+   
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <lag_h2o_ice> : ", &
+                          "<lag_h2o_ice> is set as -1 (no subsurface ice)"
+     lag_h2o_ice(:,:) = -1.
+   endif
+
+   call get_field("lag_co2_ice",lag_co2_ice,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <lag_co2_ice> : ", &
+                          "<lag_co2_ice> is set as -1 (no subsurface ice)"
+     lag_co2_ice(:,:) = -1.
+   endif
+  else
+     lag_h2o_ice(:,:) = -1.
+     lag_co2_ice(:,:) = -1.
+  endif !startphy_file
+else
+   write(*,*) 'paleo found? nostart',found
+   
+  lag_h2o_ice(:,:) = -1.
+  lag_co2_ice(:,:) = -1.
+endif !paleoclimate
+
 ! close file:
 !

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

@@ -25 +25 @@
   use time_phylmdz_mod, only: daysec
   use comslope_mod, ONLY: nslope
+  USE paleoclimate_mod, ONLY: paleoclimate, lag_h2o_ice, lag_co2_ice  
   implicit none
  
@@ -46 +47 @@
   integer :: ig
   real :: watercaptag_tmp(ngrid)
-  
+
   ! Create physics start file
   call open_restartphy(filename)
@@ -153 +154 @@
   call put_field("subslope_dist","under mesh slope distribution",subslope_dist)
 
+  ! Paleoclimate outputs
+  if(paleoclimate) then
+    call put_field("lag_h2o_ice","Depth of the H2O lags",lag_h2o_ice)
+    call put_field("lag_co2_ice","Depth of the CO2 lags",lag_co2_ice)
+  endif
   ! Close file
   call close_restartphy

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

@@ -62 +62 @@
                                      end_dust_rad_adjust_mod
       use comslope_mod, ONLY: nslope,end_comslope_h,ini_comslope_h
+      use paleoclimate_mod, ONLY: end_paleoclimate_h,ini_paleoclimate_h
       use netcdf
       USE mod_phys_lmdz_para, ONLY: is_master,bcast
@@ -186 +187 @@
       call ini_comslope_h(ngrid,nslope)
 
+      ! allocate arrays in "paleoclimate_mod"
+      call end_paleoclimate_h
+      call ini_paleoclimate_h(ngrid,nslope)
+
       END SUBROUTINE phys_state_var_init