Changeset 3910 for trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90

Timestamp:

Sep 3, 2025, 11:19:36 PM (4 months ago)

Author:

tbertrand

Message:

Pluto PCM:

• adding option to use N2 ice fractional maps (n2frac) read in the startfi.nc
• adding option in newstart to correct (tsurf and tsoil) for too warm or too cold N2-free (correct_t_non2) or N2-rich (correct_t_n2) patches

TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90

@@ -23 +23 @@
   use tabfi_mod, only: tabfi,ini_tab_controle_dyn_xios,tab_cntrl_mod
   USE tracer_h, ONLY: noms
-  USE surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam, zthe
+  USE surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam, zthe, n2frac
   use iostart, only: nid_start, open_startphy, close_startphy, &
                      get_field, get_var, inquire_field, &
@@ -91 +91 @@
 !
 IF (.not. ALLOCATED(albedodat)) ALLOCATE(albedodat(ngrid))
+IF (.not. ALLOCATED(n2frac)) ALLOCATE(n2frac(ngrid))
 IF (.not. ALLOCATED(phisfi)) ALLOCATE(phisfi(ngrid))
 IF (.not. ALLOCATED(zmea)) ALLOCATE(zmea(ngrid))
@@ -140 +141 @@
 if (is_master) write(*,*) "phyetat0: Bare ground albedo <albedodat> range:", &
              minval(albedodat), maxval(albedodat)
+
+!n2frac
+
+if (startphy_file) then
+  ! Load bare ground albedo:
+  call get_field(nid_start,"n2frac",n2frac,found)
+  if (.not.found) then
+    call abort_physic(modname,"Failed loading <n2frac>",1)
+  endif
+else
+  ! If no startfi file, use fixed value to test
+  n2frac(:)=1
+endif ! of if (startphy_file)
+if (is_master) write(*,*) "phyetat0: n2 fraction of ice  <n2frac> range:", &
+             minval(n2frac), maxval(n2frac)
 
 ! ZMEA