Changeset 2884

Timestamp:

Jan 31, 2023, 9:15:06 PM (22 months ago)

Author:

romain.vande

Message:

Mars PCM:
Watercaptag is now outputed in the starfi.nc.
There is a new mode : icelocationmode=5 to read watercaptag from startfi and to do as iceloctaionmode=4 if the variable is not present.
RV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/phyetat0_mod.F90 (modified) (3 diffs)
• libf/phymars/phyredem.F90 (modified) (3 diffs)
• libf/phymars/surfini.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3849 +3849 @@
 == 31/01/2023 == RV
 Small fix : move an endif misplaced for watercap
+
+== 31/01/2023 == RV
+Watercaptag is now outputed in the starfi.nc.
+There is a new mode : icelocationmode=5 to read watercaptag from startfi and to do as iceloctaionmode=4 if the variable is not present.

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

@@ -14 +14 @@
   use tracer_mod, only: noms ! tracer names
   use surfdat_h, only: phisfi, albedodat, z0, z0_default,&
-                       zmea, zstd, zsig, zgam, zthe, hmons, summit, base
+                       zmea, zstd, zsig, zgam, zthe, hmons, summit, base,&
+                       watercaptag
   use iostart, only: nid_start, open_startphy, close_startphy, &
                      get_field, get_var, inquire_field, &
@@ -102 +103 @@
       REAL :: surfemis  ! constant emissivity when no startfi
       REAL :: surfalbedo  ! constant albedo when no startfi
+
+      REAL :: watercaptag_tmp(ngrid)
+
       CHARACTER(len=5) :: modname="phyetat0"
 
@@ -641 +645 @@
   call soil_settings(nid_start,ngrid,nsoil,tsurf,tsoil,indextime)
 endif ! of if (startphy_file)
+
+if (startphy_file) then
+   call get_field("watercaptag",watercaptag_tmp,found,indextime)
+   if (.not.found) then
+     write(*,*) "phyetat0: Failed loading <watercaptag> : ", &
+                          "<watercaptag> is set as defined by icelocationmode in surfini.F"
+     watercaptag(:)=.false.
+   else
+     do ig=1,ngrid
+       if(watercaptag_tmp(ig).lt.0.5) then
+          watercaptag(ig)=.false.
+       else
+          watercaptag(ig)=.true.
+       endif
+     enddo
+   endif
+endif !startphy_file
+
 !
 ! close file:

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

@@ -14 +14 @@
                        z0_default, albedice, emisice, emissiv, &
                        iceradius, dtemisice, phisfi, z0,   &
-                       hmons,summit,base
+                       hmons,summit,base,watercaptag
   use dimradmars_mod, only: tauvis
   use iostart, only : open_restartphy, close_restartphy, & 
@@ -40 +40 @@
 
   real :: tab_cntrl(length) ! nb "length=100" defined in iostart module
+  integer :: ig
+  real :: watercaptag_tmp(ngrid)
   
   ! Create physics start file
@@ -133 +135 @@
   call put_field("z0","Surface roughness length",z0)
   
+  do ig=1,ngrid
+    if(watercaptag(ig)) then
+       watercaptag_tmp(ig)=1
+    else
+       watercaptag_tmp(ig)=0
+    endif
+  enddo
+
+  call put_field("watercaptag","Infinite water reservoir",watercaptag_tmp)
+
   ! Close file
   call close_restartphy

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

@@ -43 +43 @@
 ! resolution, and easily adaptable for dynamico 
 ! For visualisation : > /u/tnalmd/bin/watercaps gcm_txt_output_file
-      INTEGER,SAVE :: icelocationmode = 4
+      INTEGER,SAVE :: icelocationmode = 5
       
 !$OMP THREADPRIVATE(icelocationmode)
@@ -118 +118 @@
 c=======================================================================
 ! Initialize watercaptag (default is false)
-      watercaptag_glo(:)=.false.
+!      watercaptag_glo(:)=.false. !Already done in phyetat0 if needed
 
 c     water ice outliers
@@ -142 +142 @@
          
          if (ngrid .ne. 1) then
-           watercaptag(:) = .false.
+!           watercaptag(:) = .false.
            longwatercaptag(:) = .false.
          endif
@@ -159 +159 @@
       call gather(latitude,lati_glo)
       call gather(longitude,long_glo)
+      call gather(watercaptag,watercaptag_glo)
 
       if (is_master) then
@@ -447 +448 @@
 !         print*,'watercaptag_glo(:), ',watercaptag_glo(:)
 
+        ELSE IF (icelocationmode .eq. 5) THEN
+      
+         print*,'icelocationmode = 5'
+         print*,'Surfini: ice caps defined using startfi.nc data'
+         do ig=1,ngrid
+           if(any(watercaptag_glo)) then
+           else
+              call locate_watercaptag(klon_glo,lati_glo,
+     &            long_glo,watercaptag_glo)
+           endif
+         enddo
+
+!         print*,'watercaptag_glo(:), ',watercaptag_glo(:)
+
         ELSE
       
          print*, 'In surfini.F, icelocationmode is ', icelocationmode
-         print*, 'It should be 1, 2, 3 or 4 (default is 4)'
+         print*, 'It should be 1, 2, 3, 4 or 5 (default is 5)'
          call abort_physic("surfini","wrong icelocationmode",1)