Changeset 551 for trunk/LMDZ.MARS/libf/phymars/inifis.F

Timestamp:

Mar 1, 2012, 6:17:09 PM (13 years ago)

Author:

acolaitis

Message:

Added a fail-safe in inifis to avoid wrong configurations with callrichsl and calladj without calltherm.

File:

• trunk/LMDZ.MARS/libf/phymars/inifis.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/inifis.F

@@ -249 +249 @@
      $              "1-> new correction",
      $              "(matters only if callnirco2=T)"
+         nircorr=0
          call getin("nircorr",nircorr)
          write(*,*) " nircorr = ",nircorr
@@ -287 +288 @@
           print*,'you use the new surface layer scheme '
           print*,'by setting callrichsl=T '
+         endif
+
+         if (calladj .and. callrichsl) then
+          print*,'You should not be calling the convective adjustment
+     & scheme with the Richardson surface-layer. This approach is not
+     & physically consistent and can lead to unrealistic friction
+     & values.'
+          print*,'If you want to use the Ri. surface-layer, either
+     & activate thermals OR de-activate the convective adjustment.'
+          stop
          endif