Changeset 2665 for LMDZ5/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90

Timestamp:

Oct 12, 2016, 2:53:20 PM (8 years ago)

Author:

dcugnet

Message:

• A (re)startphy.nc file (standard name: "startphy0.nc") can be read by ce0l to get land mask, so mask can be defined (in decreasing priority order) from: 1) "o2a.nc file" if this file is found 2) "startphy0.nc" if this file is found 3) "Relief.nc" otherwise
• Sub-cell scales parameters for orographic gravity waves can be read from file "oro_params.nc" if the configuration key "read_orop" is TRUE. The effect is to bypass the "grid_noro" routine in ce0l, so that any pre-defined mask (from o2a.nc or startphy0.nc) is then overwritten.
• The gcm stops if the "limit.nc" records number differs from the current year number of days. A warning is issued in case the gcm calendar does not match the time axis attribute "calendar" (if available) from the "limit.nc" file. This attribute is now added to the "limit.nc" time axis.
• Few simplifications in grid_noro
• Few parameters changes in acama_gwd and flott_gwd.
• Variable d_u can be saved in the outputs.

File:

• LMDZ5/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90

@@ -12 +12 @@
 !     * "masque" can be:
 !       - read from file "o2a.nc"          (for coupled runs).
+!       - read from file "startphy0.nc"    (from a previous run).
 !       - created in etat0phys or etat0dyn (for forced  runs).
 !     It is then passed to limit_netcdf to ensure consistancy.
@@ -20 +21 @@
   USE etat0phys,      ONLY: etat0phys_netcdf
   USE limit,          ONLY: limit_netcdf
-  USE netcdf,         ONLY: NF90_OPEN, NF90_NOWRITE, NF90_CLOSE, NF90_NOERR
+  USE netcdf,         ONLY: NF90_OPEN, NF90_NOWRITE, NF90_CLOSE, NF90_NOERR,    &
+         NF90_INQUIRE_DIMENSION, NF90_INQ_DIMID, NF90_INQ_VARID, NF90_GET_VAR
   USE infotrac,       ONLY: type_trac, infotrac_init
   USE dimphy,         ONLY: klon
@@ -60 +62 @@
   REAL, ALLOCATABLE  :: lat_omask(:,:), dlat_omask(:), ocetmp (:,:)
   REAL               :: date, lev(1)
+
+!--- Local variables for land mask from startphy0 file reading
+  INTEGER            :: nid_sta, nid_nph, nid_msk, nphys
+  REAL, ALLOCATABLE  :: masktmp(:)
+
 #ifndef CPP_PARA
 ! for iniphysiq in serial mode
@@ -133 +140 @@
   ENDIF
 
-!--- LAND MASK. TWO CASES:
+!--- LAND MASK. THREE CASES:
 !   1) read from ocean model    file "o2a.nc"    (coupled runs)
-!   2) computed from topography file "Relief.nc" (masque(:,:)=-99999.)
-! Coupled simulations (case 1) use the ocean model mask to compute the
+!   2) read from previous run   file="startphy0.nc"
+!   3) computed from topography file "Relief.nc" (masque(:,:)=-99999.)
+! In the first case, the mask from the ocean model is used compute the
 ! weights to ensure ocean fractions are the same for atmosphere and ocean.
 !*******************************************************************************
-  IF(NF90_OPEN("o2a.nc", NF90_NOWRITE, nid_o2a)/=NF90_NOERR) THEN
-    WRITE(lunout,*)'BEWARE !! No ocean mask "o2a.nc" file found'
-    WRITE(lunout,*)'Forced run.'
-    masque(:,:)=-99999.
-  ELSE
+  IF(NF90_OPEN("o2a.nc", NF90_NOWRITE, nid_o2a)==NF90_NOERR) THEN
     iret=NF90_CLOSE(nid_o2a)
     WRITE(lunout,*)'BEWARE !! Ocean mask "o2a.nc" file found'
@@ -175 +179 @@
     masque(iip1 ,:)=masque(1,:)
     DEALLOCATE(ocemask)
+  ELSE IF(NF90_OPEN("startphy0.nc", NF90_NOWRITE, nid_sta)==NF90_NOERR) THEN
+    WRITE(lunout,*)'BEWARE !! File "startphy0.nc" found.'
+    WRITE(lunout,*)'Getting the land mask from a previous run.'
+    iret=NF90_INQ_DIMID(nid_sta,'points_physiques',nid_nph)
+    iret=NF90_INQUIRE_DIMENSION(nid_sta,nid_nph,len=nphys)
+    IF(nphys/=klon) THEN
+      WRITE(lunout,*)'Mismatching dimensions for land mask'
+      WRITE(lunout,*)'nphys  = ',nphys ,' klon = ',klon
+      iret=NF90_CLOSE(nid_sta)
+      CALL abort_gcm(modname,'',1)
+    END IF
+    ALLOCATE(masktmp(klon))
+    iret=NF90_INQ_VARID(nid_sta,'masque',nid_msk)
+    iret=NF90_GET_VAR(nid_sta,nid_msk,masktmp)
+    iret=NF90_CLOSE(nid_sta)
+    CALL gr_fi_dyn(1,klon,iip1,jjp1,masktmp,masque)
+    IF(prt_level>=1) THEN
+      WRITE(fmt,"(i4,'i1)')")iip1 ; fmt='('//ADJUSTL(fmt)
+      WRITE(lunout,*)'LAND MASK :'
+      WRITE(lunout,fmt) NINT(masque)
+    END IF
+    DEALLOCATE(masktmp)
+  ELSE
+    WRITE(lunout,*)'BEWARE !! No ocean mask "o2a.nc" file or "startphy0.nc" file found'
+    WRITE(lunout,*)'Land mask will be built from the topography file.'
+    masque(:,:)=-99999.
   END IF
   phis(:,:)=-99999.