Changeset 4223 for dynamico_lmdz/simple_physics/phyparam/physics

Timestamp:

Jan 10, 2020, 12:09:43 AM (6 years ago)

Author:

dubos

Message:

simple_physics : cleanup phyparam and iniphyparam

Location:

dynamico_lmdz/simple_physics/phyparam/physics

Files:

• callkeys.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/callkeys.F90)
• iniphyparam_mod.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/iniphyparam.F90)
• phyparam_mod.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/phyparam.F90) (7 diffs)
• radiative_lw.F90 (modified) (1 diff)
• read_param_mod.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/physics/read_param.F90)
• writefield_mod.F90 (added)

dynamico_lmdz/simple_physics/phyparam/physics/phyparam_mod.F90

@@ -27 +27 @@
 
 CONTAINS
-  
-  SUBROUTINE phyparam(ngrid,nlayer,nq, &
-       &            firstcall,lastcall, &
+
+  SUBROUTINE phyparam(ngrid,nlayer,             &
+       &            firstcall,lastcall,         &
        &            rjourvrai,gmtime,ptimestep, &
-       &            pplev,pplay,pphi,pphis,presnivs, &
-       &            pu,pv,pt,pq, &
-       &            pw, &
-       &            pdu,pdv,pdt,pdq,pdpsrf)
+       &            pplev,pplay,pphi,           &
+       &            pu,pv,pt,                   &
+       &            pdu,pdv,pdt,pdpsrf)
     USE phys_const, ONLY : g, rcp, r, unjours
     USE surface,    ONLY : soil
     USE turbulence, ONLY : vdif
     USE convection, ONLY : convadj
+    USE writefield_mod, ONLY : writefield
 
     !=======================================================================
@@ -53 +53 @@
     INTEGER, INTENT(IN) ::      &
          ngrid,                 & ! Size of the horizontal grid.
-         nlayer,                & ! Number of vertical layers.
-         nq                       ! Number of advected fields (tracers)
+         nlayer                   ! Number of vertical layers.
     LOGICAL, INTENT(IN) ::      &
          firstcall,             & ! True at the first call 
@@ -65 +64 @@
          pplay(ngrid,nlayer),   & ! Pressure at the middle of the layers (Pa)
          pphi(ngrid,nlayer),    & ! Geopotential at the middle of the layers (m2s-2)
-         pphis(ngrid),          & ! surface geopotential (unused)
-         presnivs(nlayer),      & 
          pu(ngrid,nlayer),      & ! u component of the wind (ms-1)
          pv(ngrid,nlayer),      & ! v component of the wind (ms-1)
-         pw(ngrid,nlayer),      & ! vertical velocity (unused)
-         pt(ngrid,nlayer),      & ! Temperature (K)
-         pq(ngrid,nlayer,nq)      ! Advected fields (unused)
+         pt(ngrid,nlayer)         ! Temperature (K)
     REAL, INTENT(OUT)   ::      & ! output : physical tendencies
          pdu(ngrid,nlayer),     &
          pdv(ngrid,nlayer),     &
          pdt(ngrid,nlayer),     &
-         pdq(ngrid,nlayer,nq),  &
          pdpsrf(ngrid)
     
     !    Local variables :
     REAL, DIMENSION(ngrid) :: mu0,fract
-    INTEGER :: j,l,ig,ierr,aslun,nlevel,igout,it1,it2,isoil,iq
-    INTEGER*4 day_ini
+    INTEGER :: j,l,ig,nlevel,igout
     !
     REAL :: zday, zdtime
@@ -96 +89 @@
     REAL zdtlw(ngrid,nlayer),zdtsw(ngrid,nlayer)
     REAL zfluxsw(ngrid),zfluxlw(ngrid)
-    REAL factq(nq),tauq(nq)
-    character*3 nomq
-    
-    !   Local saved variables:
-    !   ----------------------
     
     print*,'OK DANS PHYPARAM'
-    print*,'nq ',nq
     print*,'latitude0',ngrid,lati(1:2),lati(ngrid-1:ngrid)
     print*,'nlayer',nlayer
-    print*,'size pdq ',ngrid*nlayer*4,ngrid*nlayer*nq, &
-         &      size(pdq),size(lati),size(pq),size(factq)
     
     IF (ngrid.NE.ngridmax) THEN
@@ -142 +127 @@
           ENDDO
        ENDIF
-       
-       PRINT*,'FIRSTCALL B '
-       print*,'INIIO avant iophys_ini '
-       call iophys_ini('phys.nc    ',presnivs)       
     ENDIF
 
@@ -154 +135 @@
     icount=icount+1
     
-    pdq(:,:,:) = 0. ! we do not use tracers in this physics package    
     pdv(:,:)  = 0.
     pdu(:,:)  = 0.
@@ -315 +295 @@
        enddo
        
-       call iophys_ecrit('u',nlayer,'Vent zonal moy','m/s',pu)
-       call iophys_ecrit('v',nlayer,'Vent meridien moy','m/s',pv)
-       call iophys_ecrit('temp',nlayer,'Temperature','K',pt)
-       call iophys_ecrit('geop',nlayer,'Geopotential','m2/s2',pphi)
-       call iophys_ecrit('plev',nlayer,'plev','Pa',pplev(:,1:nlayer))
-       
-       call iophys_ecrit('du',nlayer,'du',' ',pdu)
-       call iophys_ecrit('dv',nlayer,'du',' ',pdv)
-       call iophys_ecrit('dt',nlayer,'du',' ',pdt)
-       call iophys_ecrit('dtsw',nlayer,'dtsw',' ',zdtsw)
-       call iophys_ecrit('dtlw',nlayer,'dtlw',' ',zdtlw)
-       
-       do iq=1,nq
-          nomq="tr."
-          write(nomq(2:3),'(i1.1)') iq
-          call iophys_ecrit(nomq,nlayer,nomq,' ',pq(:,:,iq))
-       enddo
-       
-       call iophys_ecrit('ts',1,'Surface temper','K',tsurf)
-       call iophys_ecrit('coslon',1,'coslon',' ',coslon)
-       call iophys_ecrit('sinlon',1,'sinlon',' ',sinlon)
-       call iophys_ecrit('coslat',1,'coslat',' ',coslat)
-       call iophys_ecrit('sinlat',1,'sinlat',' ',sinlat)
-       call iophys_ecrit('mu0',1,'mu0',' ',mu0)
-       call iophys_ecrit('alb',1,'alb',' ',albedo)
-       call iophys_ecrit('fract',1,'fract',' ',fract)
-       call iophys_ecrit('ps',1,'Surface pressure','Pa',pplev)
-       call iophys_ecrit('slp',1,'Sea level pressure','Pa',zpmer)
-       call iophys_ecrit('swsurf',1,'SW surf','Pa',zfluxsw)
-       call iophys_ecrit('lwsurf',1,'LW surf','Pa',zfluxlw)
+       call writefield('u','Vent zonal moy','m/s',pu)
+       call writefield('v','Vent meridien moy','m/s',pv)
+       call writefield('temp','Temperature','K',pt)
+       call writefield('geop','Geopotential','m2/s2',pphi)
+       call writefield('plev','plev','Pa',pplev(:,1:nlayer))
+       
+       call writefield('du','du',' ',pdu)
+       call writefield('dv','du',' ',pdv)
+       call writefield('dt','du',' ',pdt)
+       call writefield('dtsw','dtsw',' ',zdtsw)
+       call writefield('dtlw','dtlw',' ',zdtlw)
+       
+       call writefield('ts','Surface temper','K',tsurf)
+       call writefield('coslon','coslon',' ',coslon)
+       call writefield('sinlon','sinlon',' ',sinlon)
+       call writefield('coslat','coslat',' ',coslat)
+       call writefield('sinlat','sinlat',' ',sinlat)
+       call writefield('mu0','mu0',' ',mu0)
+       call writefield('alb','alb',' ',albedo)
+       call writefield('fract','fract',' ',fract)
+       call writefield('ps','Surface pressure','Pa',pplev(:,1))
+       call writefield('slp','Sea level pressure','Pa',zpmer)
+       call writefield('swsurf','SW surf','Pa',zfluxsw)
+       call writefield('lwsurf','LW surf','Pa',zfluxlw)
        
     endif
-    
-    !-----------------------------------------------------------------------
-    IF(lastcall) THEN
-       call iotd_fin
-       PRINT*,'Ecriture du fichier de reinitialiastion de la physique'
-       write(75)  tsurf,tsoil
-    ENDIF    
     
   END SUBROUTINE phyparam

dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90

@@ -68 +68 @@
     REAL zup(ngrid,nlayer+1),zdup(ngrid)
 
-    CHARACTER(:), PARAMETER :: tag='rad/lw'
+    CHARACTER(6), PARAMETER :: tag='rad/lw'
     !-----------------------------------------------------------------------
     !   initialisations: