Changeset 520 for trunk/LMDZ.MARS/libf/phymars/aeropacity.F

Timestamp:

Feb 10, 2012, 12:33:14 PM (13 years ago)

Author:

tnavarro

Message:

10/02/12 == TN

Major update on watercycle: a smaller integration timestep is now used

in watercloud.F, sedimentation of clouds is done in watercloud instead of
callsedim.F

Temperature-dependant contact parameter in nuclea.F
No dust lifting if CO2 ice in vdif.c
Ice integrated column opacity is written in diagfi from physiq.F, instead

of aeropacity.F. Mandatory if iradia is not 1.

New definition of permanent ice in surfini.F and possibility to have an ice

cap in it in 1d.

Update in deftank: callphys.def.outliers,run.def.1d; added traceur.def.scavenging

File:

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

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

@@ -1 +1 @@
       SUBROUTINE aeropacity(ngrid,nlayer,nq,zday,pplay,pplev,ls,
-     &    pq,tauscaling,tauref,tau,aerosol,reffrad,nueffrad,
+     &    pq,tauscaling,tauref,tau,taucloudtes,aerosol,reffrad,nueffrad,
      &    QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d)
                                                    
@@ -386 +386 @@
           ENDDO
         ENDDO
-c       3. Outputs
-        IF (ngrid.NE.1) THEN
-          CALL WRITEDIAGFI(ngridmx,'tauVIS','tauext VIS refwvl',
-     &      ' ',2,taucloudvis)
-          CALL WRITEDIAGFI(ngridmx,'tauTES','tauabs IR refwvl',
-     &      ' ',2,taucloudtes)
-          IF (callstats) THEN
-            CALL wstats(ngridmx,'tauVIS','tauext VIS refwvl',
-     &        ' ',2,taucloudvis)
-            CALL wstats(ngridmx,'tauTES','tauabs IR refwvl',
-     &        ' ',2,taucloudtes)
-          ENDIF
-        ELSE
-c         CALL writeg1d(ngrid,1,taucloudtes,'tautes','NU')
-        ENDIF
+c       3. Outputs -- Now done in physiq.F
+!        IF (ngrid.NE.1) THEN
+!          CALL WRITEDIAGFI(ngridmx,'tauVIS','tauext VIS refwvl',
+!     &      ' ',2,taucloudvis)
+!          CALL WRITEDIAGFI(ngridmx,'tauTES','tauabs IR refwvl',
+!     &      ' ',2,taucloudtes)
+!          IF (callstats) THEN
+!            CALL wstats(ngridmx,'tauVIS','tauext VIS refwvl',
+!     &        ' ',2,taucloudvis)
+!            CALL wstats(ngridmx,'tauTES','tauabs IR refwvl',
+!     &        ' ',2,taucloudtes)
+!          ENDIF
+!        ELSE
+!         CALL writeg1d(ngrid,1,taucloudtes,'tautes','NU')
+!        ENDIF
 c==================================================================
         END SELECT aerkind