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

Timestamp:

Jan 7, 2014, 10:40:56 AM (11 years ago)

Author:

milmd

Message:

LMDZ.GENERIC. added a condition to avoid a small bug in mucorr when running in 1D. but 'correction de rotondite' seems not generic at all in mucorr although the impact is small because hardcoded correction is pretty negligible (AS).

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• mucorr.F (modified) (1 diff)
• physiq.F90 (modified) (1 diff)

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

@@ -92 +92 @@
 c   ------------------------
 
+      ! condition added to avoid errors when rad is not set (e.g. 1D runs)
+      IF (prad.ne.0) THEN
+  
       alph=phaut/prad
       DO 30 j=1,npts
 c !!!!!!
+ !!!!!!! AS: how generic is this???
          pmu(j)=sqrt(1224.*pmu(j)*pmu(j)+1.)/35.
 c    $          (sqrt(alph*alph*pmu(j)*pmu(j)+2.*alph+1.)-alph*pmu(j))
 30    CONTINUE
 
+      ENDIF
+
       RETURN
       END

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

@@ -700 +700 @@
 
                call mucorr(ngrid,declin,lati,mu0,fract,10000.,rad)
-               ! WARNING: this function appears not to work in 1D
 
            endif