Changeset 2798 for LMDZ5/trunk/libf

Timestamp:

Feb 24, 2017, 3:40:57 PM (8 years ago)

Author:

jghattas

Message:

Add option no_ter_antartique that will remove all land fractions souther than 60deg south. This option is only used by ce0l.

File:

• LMDZ5/trunk/libf/dynphy_lonlat/phylmd/init_ssrf_m.F90 (modified) (3 diffs)

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

@@ -9 +9 @@
   USE grid_atob_m,        ONLY: grille_m
   USE ioipsl,             ONLY: flininfo, flinopen, flinget, flinclo
+  USE ioipsl_getin_p_mod, ONLY: getin_p
   USE comconst_mod, ONLY: im, pi
 
@@ -42 +43 @@
   REAL, ALLOCATABLE :: dlat_lic(:), lat_lic(:,:), flic_tmp(:,:), vtmp(:,:)
   REAL              :: date, lev(1), dt, deg2rad
+  LOGICAL           :: no_ter_antartique   ! If true, no land points are allowed at Antartic
 !-------------------------------------------------------------------------------
   deg2rad= pi/180.0
@@ -97 +99 @@
   END DO 
 
+
+  !--- Option no_ter_antartique removes all land fractions souther than 60S.
+  !--- Land ice is set instead of the land fractions on these latitudes.
+  !--- The ocean and sea-ice fractions are not changed.
+  no_ter_antartique=.FALSE.
+  CALL getin_p('no_ter_antartique',no_ter_antartique)
+  WRITE(lunout,*)"no_ter_antartique=",no_ter_antartique
+  IF (no_ter_antartique) THEN
+     ! Remove all land fractions souther than 60S and set land-ice instead
+     WRITE(lunout,*) "Remove land fractions souther than 60deg south by increasing"
+     WRITE(lunout,*) "the continental ice fractions. No land can now be found at Antartic."
+     DO ji=1, klon
+        IF (latitude_deg(ji)<-60.0) THEN
+           pctsrf(ji,is_lic) = pctsrf(ji,is_lic) + pctsrf(ji,is_ter)
+           pctsrf(ji,is_ter) = 0
+        END IF
+     END DO
+  END IF
+
+
 !--- Sub-surface ocean and sea ice (sea ice set to zero for start).
   pctsrf(:,is_oce)=(1.-zmasq(:))