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

Timestamp:

Jul 19, 2013, 10:25:38 AM (12 years ago)

Author:

emillour

Message:

Generic GCM:

• update the sponge layer: trun it into a module and (more important) compute the sponge quenching analytically rather than via Forward Euler approximation.

EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/callsedim.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/callsedim.F

@@ -75 +75 @@
       IF (firstcall) THEN
         firstcall=.false.
+        ! add some tests on presence of required tracers/aerosols:
+        if (water) then
+          if (igcm_h2o_ice.eq.0) then
+            write(*,*) "callsedim error: water=.true.",
+     &                 " but igcm_h2o_ice=0"
+          stop
+          endif
+          if (iaero_h2o.eq.0) then
+            write(*,*) "callsedim error: water=.true.",
+     &                 " but iaero_ho2=0"
+          stop
+          endif
+        endif
       ENDIF ! of IF (firstcall)