Changeset 1983 for trunk/LMDZ.MARS/libf/phymars/lwmain_mod.F

Timestamp:

Aug 1, 2018, 3:35:24 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:
Cosmetic/practical changes:

• swmain and lwmain become modules swmain_mod, lwmain_mod
• Addition of the intent in/out characteristics of variables in swmain_mod and lwmain_mod subroutines

Correction:

• in callsedim_mod, declaration of variable tau(ngrid,nlay) corrected to tau(ngrid,naerkind)

MV

File:

• trunk/LMDZ.MARS/libf/phymars/lwmain_mod.F (moved) (moved from trunk/LMDZ.MARS/libf/phymars/lwmain.F) (5 diffs)

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

@@ +1 @@
+       MODULE lwmain_mod
+
+       IMPLICIT NONE
+
+       CONTAINS
+       
        subroutine lwmain (ig0,icount,kdlon,kflev
      .                   ,dp,dt0,emis
@@ -16 +22 @@
       use yomlw_h, only: nlaylte, xi
       implicit none
- 
-#include "callkeys.h"
-#include "comg1d.h"
+      
+c     declarations
+c     -------------   
+      include "callkeys.h"
+      include "comg1d.h"
 
 c----------------------------------------------------------------------
 c         0.1   arguments
 c               ---------
-c                                                            inputs:
-c                                                            -------
-      integer ig0
-      integer icount
-      integer kdlon            ! part of ngrid
-      integer kflev            ! part of nlayer
-
-      real dp (ndlo2,kflev)         ! layer pressure thickness (Pa)
-      real dt0 (ndlo2)              ! surface temperature discontinuity (K)
-      real emis (ndlo2)             ! surface emissivity
-      real plev (ndlo2,kflev+1)     ! level pressure (Pa)
-      real tlev (ndlo2,kflev+1)     ! level temperature (K)
-      real tlay (ndlo2,kflev)       ! layer temperature (K)
-      real aerosol(ndlo2,kflev,naerkind)      !  aerosol extinction optical
+c     inputs/outputs:
+c     -------
+      integer, intent(in) :: ig0
+      integer, intent(in) :: icount
+      integer, intent(in) :: kdlon            ! part of ngrid
+      integer, intent(in) :: kflev            ! part of nlayer
+
+      real, intent(in) :: dp (ndlo2,kflev)         ! layer pressure thickness (Pa)
+      real, intent(in) :: dt0 (ndlo2)              ! surface temperature discontinuity (K)
+      real, intent(in) :: emis (ndlo2)             ! surface emissivity
+      real, intent(in) :: plev (ndlo2,kflev+1)     ! level pressure (Pa)
+      real, intent(in) :: tlev (ndlo2,kflev+1)     ! level temperature (K)
+      real, intent(in) :: tlay (ndlo2,kflev)       ! layer temperature (K)
+      real, intent(in) :: aerosol(ndlo2,kflev,naerkind)      !  aerosol extinction optical
 c                         depth at reference wavelength "longrefvis" set
 c                         in dimradmars_mod , in each layer, for one of
@@ -44 +52 @@
 c                                                            outputs:
 c                                                            --------
-      real coolrate(ndlo2,kflev)      ! cooling rate (K/s)
-      real fluxground(ndlo2)          ! downward ground flux (W/m2)
-      real fluxtop(ndlo2)             ! outgoing upward flux (W/m2) ("OLR")
-      real netrad (ndlo2,kflev)       ! radiative budget (W/m2)
+      real, intent(out) :: coolrate(ndlo2,kflev)      ! cooling rate (K/s)
+      real, intent(out) :: fluxground(ndlo2)          ! downward ground flux (W/m2)
+      real, intent(out) :: fluxtop(ndlo2)             ! outgoing upward flux (W/m2) ("OLR")
+      real, intent(out) :: netrad (ndlo2,kflev)       ! radiative budget (W/m2)
 c     Aerosol optical properties
-      REAL :: QIRsQREF3d(ndlo2,kflev,nir,naerkind)
-      REAL :: omegaIR3d(ndlo2,kflev,nir,naerkind)
-      REAL :: gIR3d(ndlo2,kflev,nir,naerkind)
-
+      real, intent(in) :: QIRsQREF3d(ndlo2,kflev,nir,naerkind)
+      real, intent(in) :: omegaIR3d(ndlo2,kflev,nir,naerkind)
+      real, intent(in) :: gIR3d(ndlo2,kflev,nir,naerkind)
+      real, intent(in) :: co2ice(ndlo2)           ! co2 ice surface layer (kg.m-2)
 c----------------------------------------------------------------------
 c         0.2   local arrays
@@ -75 +83 @@
       real newcoolrate(ndlon,nflev) ! cooling rate (K/s) / with implicite scheme
 
-      REAL co2ice(ndlo2)           ! co2 ice surface layer (kg.m-2)
-      REAL emis_gaz(ndlo2)         ! emissivity for gaz computations
+      real emis_gaz(ndlo2)         ! emissivity for gaz computations
 
       integer jk,jkk,ja,jl
@@ -193 +200 @@
 c----------------------------------------------------------------------
 
-      return
-      end
+      END SUBROUTINE lwmain
+
+      END MODULE lwmain_mod 
+