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

Timestamp:

May 5, 2023, 3:15:04 PM (19 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

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F (modified) (7 diffs)

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)