Changeset 4595 for LMDZ6/trunk/libf/phylmd/create_etat0_unstruct.F90

Timestamp:

Jun 30, 2023, 12:26:07 PM (11 months ago)

Author:

jghattas

Message:

Option no_ter_antartique is only read and used for the case landice_opt < 2. This was already the case on the regular grid.

File:

• LMDZ6/trunk/libf/phylmd/create_etat0_unstruct.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/create_etat0_unstruct.F90

@@ -147 +147 @@
   !--- 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
+  !--- This option is only available if landice_opt<2.   
+  IF (landice_opt .LT. 2) THEN
+     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
   END IF