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

Timestamp:

May 5, 2023, 3:15:04 PM (20 months ago)

Author:

romain.vande

Message:

Mars PCM :
Correct start2archive to write watercaptag correctly.
Watercaptag will be set to false and correctly handle by the PCM in the case where we change resolution.
+ Correct inertiesoil writting in start2archive
RV

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• lect_start_archive.F (modified) (7 diffs)
• start2archive.F (modified) (10 diffs)
• start2archive_SSO.F (modified) (6 diffs)
• write_archive.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F

@@ -26 +26 @@
      &             subslope_dist,end_comslope_h,ini_comslope_h
       use netcdf
+      use surfdat_h, ONLY: watercaptag
       implicit none
 
@@ -139 +140 @@
       real totcloudfracS(iip1,jjp1)
       real watercapS(iip1,jjp1,nslope)
+      real watercaptagS(iip1,jjp1)
       real albedoS(iip1,jjp1,nslope)
 
@@ -177 +179 @@
       real, dimension(:,:,:), allocatable :: watercapold
       real, dimension(:,:), allocatable :: watercapold_noslope
+      real, dimension(:,:), allocatable :: watercaptagold
+      real, dimension(:), allocatable :: watercaptag_tmp
       real, dimension(:,:,:), allocatable :: albedoold
       real, dimension(:,:), allocatable :: albedoold_noslope
@@ -426 +430 @@
       allocate(totcloudfracold(imold+1,jmold+1))
       allocate(watercapold(imold+1,jmold+1,nslope))
+      allocate(watercaptagold(imold+1,jmold+1))
+      allocate(watercaptag_tmp(ngrid))
       allocate(albedoold(imold+1,jmold+1,nslope))
       
@@ -908 +914 @@
         IF (ierr .NE. NF_NOERR) THEN
            PRINT*, "lect_start_archive: Failed loading <watercap>"
+           PRINT*, NF_STRERROR(ierr)
+           CALL abort
+        ENDIF
+      ENDIF
+
+c
+      ierr = NF_INQ_VARID (nid, "watercaptag", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: <watercaptag> not in file"
+         PRINT*, "watercaptag set to false, will be adapted in 
+     &            surfini of PCM"
+         watercaptagold(:,:) = 0
+      ELSE
+        ierr = nf90_get_var(nid, nvarid,watercaptagold)
+        IF (ierr .NE. NF_NOERR) THEN
+           PRINT*, "lect_start_archive: Failed loading <watercaptag>"
            PRINT*, NF_STRERROR(ierr)
            CALL abort
@@ -1230 +1252 @@
      &       watercap(:,1))
 
+c Watercaptag
+      if(imold.eq.iim .and. jmold.eq.jjm) then
+      else
+       print *, "We are doing an horizontal interpolation, 
+     &    watercaptag will be set to false and redefined proprely 
+     &    in the PCM(surfini)"
+       watercaptagold(:,:)=0
+      endif
+        call interp_horiz (watercaptagold(:,:),watercaptags(:,:),
+     &       imold,jmold,iim,jjm,1,rlonuold,rlatvold,rlonu,rlatv)
+        call gr_dyn_fi (1,iim+1,jjm+1,ngrid,watercaptags(:,:),
+     &       watercaptag_tmp(:))
+
+        do i=1,ngrid
+          if(watercaptag_tmp(i).gt. 0.5) then
+            watercaptag(i)=.true.
+          else
+            watercaptag(i)=.false.
+          endif
+        enddo
+
+
 c Surface albedo
       call interp_horiz (albedoold(:,:,1),albedoS(:,:,1),
@@ -1605 +1649 @@
       deallocate(totcloudfracold)
       deallocate(watercapold)
+      deallocate(watercaptagold)
       deallocate(albedoold)
       deallocate(var,varp1)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -20 +20 @@
 
       use infotrac, only: infotrac_init, nqtot, tname
-      use comsoil_h, only: nsoilmx, inertiedat
-      use surfdat_h, only: ini_surfdat_h, qsurf
+      use comsoil_h, only: nsoilmx, inertiedat, inertiesoil
+      use surfdat_h, only: ini_surfdat_h, qsurf,watercaptag
       use comsoil_h, only: ini_comsoil_h
 !      use comgeomphy, only: initcomgeomphy
@@ -66 +66 @@
       REAL,ALLOCATABLE :: tsurf(:,:)        ! Surface temperature
       REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
-      REAL,ALLOCATABLE :: inertiesoil(:,:,:) ! Soil thermal inertia (!= inertiedat which is for present day climate)
       REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
       REAL :: tauscaling(ngridmx) ! dust conversion factor
@@ -87 +86 @@
       REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate)
       REAL,ALLOCATABLE :: watercapS(:,:)
+      REAL,ALLOCATABLE :: watercaptag_tmp(:)
+      REAL,ALLOCATABLE :: watercaptagS(:)
       REAL :: tauscalingS(ip1jmp1)
       REAL :: totcloudfracS(ip1jmp1)
@@ -167 +168 @@
       allocate(tsurf(ngridmx,nslope))
       allocate(tsoil(ngridmx,nsoilmx,nslope))
-      allocate(inertiesoil(ngridmx,nsoilmx,nslope))
       allocate(watercap(ngridmx,nslope))
       allocate(emis(ngridmx,nslope))
@@ -177 +177 @@
       allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
       allocate(watercapS(ip1jmp1,nslope))
+      allocate(watercaptagS(ip1jmp1))
       allocate(emisS(ip1jmp1,nslope))
       allocate(albedoS(ip1jmp1,nslope))
@@ -207 +208 @@
        ENDIF
 
+         CALL surfini(ngridmx,qsurf)
+
       ierr = NF_CLOSE(nid1)
 
@@ -287 +290 @@
      &    tsoilS(:,:,islope))
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiesoil(:,:,islope),
-     &    inertiesoil(:,:,islope))
+     &    inertiesoilS(:,:,islope))
       ! Note: thermal inertia "inertiedat" is in comsoil.h
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis(:,islope),
@@ -302 +305 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,tauscaling,tauscalingS)
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,totcloudfrac,totcloudfracS)
+
+      allocate(watercaptag_tmp(ngridmx))
+      do ij=1,ngridmx
+        if(watercaptag(ij)) then
+        watercaptag_tmp(ij)=1
+        else
+        watercaptag_tmp(ij)=0
+        endif
+      enddo
+
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercaptag_tmp(:),
+     &    watercaptagS(:))
 
 c=======================================================================
@@ -402 +417 @@
       call write_archive(nid,ntime,'watercap','couche de glace h2o',
      &  'kg/m2',2,watercapS)
+      call write_archive(nid,ntime,'watercaptag','couche de glace h2o',
+     &  'kg/m2',2,watercaptagS)
       call write_archive(nid,ntime,'tauscaling',
      &  'dust conversion factor',' ',2,tauscalingS)
@@ -462 +479 @@
      &     'J.s-1/2.m-2.K-1',-3,inertiesoilS(:,:,:))
 ! Write soil thermal inertia for current climate
-      call write_archive(nid,ntime,'inertiedat - present day TI',
-     &     'Soil thermal inertia',
+      call write_archive(nid,ntime,'inertiedat',
+     &     'Soil thermal inertia (present day TI)',
      &     'J.s-1/2.m-2.K-1',-3,ithS)
 

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive_SSO.F

@@ -20 +20 @@
 
       use infotrac, only: infotrac_init, nqtot, tname
-      use comsoil_h, only: nsoilmx, inertiedat
+      use comsoil_h, only: nsoilmx, inertiedat,inertiesoil
       use surfdat_h, only: ini_surfdat_h, qsurf
-      use comsoil_h, only: ini_comsoil_h
 !      use comgeomphy, only: initcomgeomphy
       use filtreg_mod, only: inifilr
@@ -69 +68 @@
       REAL,ALLOCATABLE :: tsurf(:,:)        ! Surface temperature
       REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
-      REAL,ALLOCATABLE :: inertiesoil(:,:,:) ! Soil thermal inertia (!= inertiedat which is for present day climate)
       REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
       REAL tauscaling(ngridmx) ! dust conversion factor
@@ -185 +183 @@
       allocate(tsurf(ngridmx,nslope))
       allocate(tsoil(ngridmx,nsoilmx,nslope))
-      allocate(inertiesoil(ngridmx,nsoilmx,nslope))
       allocate(watercap(ngridmx,nslope))
       allocate(emis(ngridmx,nslope))
@@ -302 +299 @@
      &    tsoilS(:,:,islope))
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiesoil(:,:,islope),
-     &    inertiesoil(:,:,islope))
+     &    inertiesoilS(:,:,islope))
       ! Note: thermal inertia "inertiedat" is in comsoil.h
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis(:,islope),
@@ -313 +310 @@
      &    subslope_distS(:,islope))
       enddo
+
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiedat,ithS)
       call gr_fi_dyn(llm+1,ngridmx,iip1,jjp1,q2,q2S)
@@ -497 +495 @@
       call write_archive(nid,ntime,'tsoil','Soil temperature',
      &     'K',-3,tsoilS(:,:,:))
+
 ! Write soil thermal inertia
       call write_archive(nid,ntime,'inertiesoil','Soil TI',

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F

@@ -226 +226 @@
      &        nom.eq."ZGAM" .or. nom.eq."albedodat" .or.
      &        nom.eq."z0" .or. nom.eq."summit" .or. nom.eq."hmons" 
-     &        .or. nom.eq."base") then
+     &        .or. nom.eq."base".or. nom.eq."watercaptag") then
 
            if (ierr /= NF_NOERR) then