Changeset 2603 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jan 7, 2022, 8:27:58 AM (3 years ago)

Author:

emillour

Message:

Generic GCM:
Minor bug fix for the (rather specific) case when there is no surface and
inertiedat() is not defined.
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -48 +48 @@
       use time_phylmdz_mod, only: daysec
       use callkeys_mod, only: albedo_spectral_mode, calladj, calldifv, &
-                              calllott_nonoro, callrad, callsoil, callstats, &
+                              calllott_nonoro, callrad, callsoil, nosurf, &
+                              callstats, &
                               calltherm, CLFvarying, co2cond, corrk, diagdtau, &
                               diurnal, enertest, fat1au, flatten, j2, &
@@ -633 +634 @@
          
            rnat(:)=1.
-           do ig=1,ngrid
+           if (.not.nosurf) then ! inertiedat only defined if there is a surface
+            do ig=1,ngrid
               if(inertiedat(ig,1).gt.1.E4)then
                  rnat(ig)=0
               endif
-           enddo
-
-           print*,'WARNING! Surface type currently decided by surface inertia'
-           print*,'This should be improved e.g. in newstart.F'
-           
+            enddo
+
+            print*,'WARNING! Surface type currently decided by surface inertia'
+            print*,'This should be improved e.g. in newstart.F'
+           endif
          endif ! end of 'ok_slab_ocean'.