Changeset 3468 for trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

Timestamp:

Oct 24, 2024, 9:18:07 AM (4 weeks ago)

Author:

emillour

Message:

Mars PCM:
Remove obsolete/depreciated lwrite flag (which would trigger some very specific
extra text outputs), in code and in reference callphys.def files.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (5 diffs)

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

@@ -6 +6 @@
       
       SUBROUTINE vdifc(ngrid,nlay,nsoil,nq,nqsoil,ppopsk,
-     $                ptimestep,pcapcal,lecrit,
+     $                ptimestep,pcapcal,
      $                pplay,pplev,pzlay,pzlev,pz0,
      $                pu,pv,ph,pq,ptsrf,ptsoil,pemis,pqsurf,qsoil,
@@ -85 +85 @@
 c    Argument added for condensation:
       REAL,INTENT(IN) :: ppopsk(ngrid,nlay)
-      logical,INTENT(IN) :: lecrit
       REAL,INTENT(IN) :: pcondicea_co2microp(ngrid,nlay)! tendency due to CO2 condensation (kg/kg.s-1)
       
@@ -341 +340 @@
       ENDDO
 
-c    ** diagnostique pour l'initialisation
-c       ----------------------------------
-
-      IF(lecrit) THEN
-         ig=ngrid/2+1
-         PRINT*,'Pression (mbar) ,altitude (km),u,v,theta, rho dz'
-         DO ilay=1,nlay
-            WRITE(*,'(6f11.5)')
-     s      .01*pplay(ig,ilay),.001*pzlay(ig,ilay),
-     s      pu(ig,ilay),pv(ig,ilay),ph(ig,ilay),za(ig,ilay)
-         ENDDO
-         PRINT*,'Pression (mbar) ,altitude (km),zb'
-         DO ilev=1,nlay
-            WRITE(*,'(3f15.7)')
-     s      .01*pplev(ig,ilev),.001*pzlev(ig,ilev),
-     s      zb0(ig,ilev)
-         ENDDO
-      ENDIF
-
 c     -----------------------------------
 c     Potential Condensation temperature:
@@ -503 +483 @@
         end do
       end if
-
-c    ** diagnostique pour le schema de turbulence
-c       -----------------------------------------
-
-      IF(lecrit) THEN
-         PRINT*
-         PRINT*,'Diagnostic for the vertical turbulent mixing'
-         PRINT*,'Cd for momentum and potential temperature'
-
-         PRINT*,zcdv_tmp(ngrid/2+1),zcdh_tmp(ngrid/2+1)
-         PRINT*,'Mixing coefficient for momentum and pot.temp.'
-         DO ilev=1,nlay
-            PRINT*,zkv(ngrid/2+1,ilev),zkh(ngrid/2+1,ilev)
-         ENDDO
-      ENDIF
 
 c-----------------------------------------------------------------------
@@ -1465 +1430 @@
      &                                    *ptimestep))/ptimestep
 
-c    ** diagnostique final 
-c       ------------------
-
-      IF(lecrit) THEN
-         PRINT*,'In vdif'
-         PRINT*,'Ts (t) and Ts (t+st)'
-         WRITE(*,'(a10,3a15)')
-     s   'theta(t)','theta(t+dt)','u(t)','u(t+dt)'
-         PRINT*,ptsrf(ngrid/2+1,:),ztsrf2(ngrid/2+1)
-         DO ilev=1,nlay
-            WRITE(*,'(4f15.7)')
-     s      ph(ngrid/2+1,ilev),zhs(ngrid/2+1,ilev),
-     s      pu(ngrid/2+1,ilev),zu(ngrid/2+1,ilev)
-
-         ENDDO
-      ENDIF
-
       END SUBROUTINE vdifc