Changeset 1771 for trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90

Timestamp:

Aug 27, 2017, 8:22:48 AM (7 years ago)

Author:

aslmd

Message:

LMDZ.MARS dimradmars 1. add the possibility to deallocate common modules with subroutine (as in previous commit) 2. moved totcloudfrac allocation to dimradmars as suggested by AP it has nothing to do in physiq_mod

File:

• trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90

@@ -98 +98 @@
   REAL,SAVE,ALLOCATABLE :: albedo(:,:) ! Surface albedo in each solar band
   REAL,SAVE,ALLOCATABLE :: tauscaling(:)   ! Convertion factor for qdust and Ndust
+  REAL,SAVE,ALLOCATABLE :: totcloudfrac(:) ! total cloud fraction over the column
 ! aerosol (dust or ice) extinction optical depth  at reference wavelength 
 ! "longrefvis" set in dimradmars_mod , for one of the "naerkind"  kind of
@@ -169 +170 @@
    allocate(tauscaling(ngrid))
    allocate(nueffdust(ngrid,nlayer))
+   allocate(totcloudfrac(ngrid))
 
   end subroutine ini_dimradmars_mod
+
+  subroutine end_dimradmars_mod
+
+  implicit none
+
+   if (allocated(albedo))      deallocate(albedo)
+   if (allocated(dtrad))       deallocate(dtrad)
+   if (allocated(fluxrad_sky)) deallocate(fluxrad_sky)
+   if (allocated(fluxrad))     deallocate(fluxrad)
+   if (allocated(tauscaling))  deallocate(tauscaling)
+   if (allocated(nueffdust))   deallocate(nueffdust)
+   if (allocated(totcloudfrac))   deallocate(totcloudfrac)
+
+  end subroutine end_dimradmars_mod
+
  
   subroutine ini_scatterers(ngrid,nlayer)
@@ -179 +196 @@
   integer,intent(in) :: nlayer ! number of atmospheric layers
 
-   allocate(name_iaer(naerkind))
+   !! domain-dependent 
+   !! -- only used in physiq_mod & intent(out) in callradite
+   if (allocated(aerosol)) deallocate(aerosol)
    allocate(aerosol(ngrid,nlayer,naerkind))
-   allocate(longrefir(naerkind))
-   allocate(longrefvis(naerkind))
-   allocate(iaerdust(naerkind))
-   allocate(file_id(naerkind,2))
-
-   allocate(QVISsQREF(nsun,naerkind,nsizemax))
-   allocate(omegavis(nsun,naerkind,nsizemax))
-   allocate(gvis(nsun,naerkind,nsizemax))
-   allocate(QIRsQREF(nir,naerkind,nsizemax))
-   allocate(omegaIR(nir,naerkind,nsizemax))
-   allocate(gIR(nir,naerkind,nsizemax))
-   allocate(nsize(naerkind,2))
-   allocate(radiustab(naerkind,2,nsizemax))
-   allocate(QREFvis(naerkind,nsizemax))
-   allocate(QREFir(naerkind,nsizemax))
-   allocate(omegaREFvis(naerkind,nsizemax))
-   allocate(omegaREFir(naerkind,nsizemax))
+
+   !! not domain-dependent
+   if (.not.allocated(name_iaer)) allocate(name_iaer(naerkind))
+   if (.not.allocated(longrefir)) allocate(longrefir(naerkind))
+   if (.not.allocated(longrefvis)) allocate(longrefvis(naerkind))
+   if (.not.allocated(iaerdust)) allocate(iaerdust(naerkind))
+   if (.not.allocated(file_id)) allocate(file_id(naerkind,2))
+   if (.not.allocated(QVISsQREF)) allocate(QVISsQREF(nsun,naerkind,nsizemax))
+   if (.not.allocated(omegavis)) allocate(omegavis(nsun,naerkind,nsizemax))
+   if (.not.allocated(gvis)) allocate(gvis(nsun,naerkind,nsizemax))
+   if (.not.allocated(QIRsQREF)) allocate(QIRsQREF(nir,naerkind,nsizemax))
+   if (.not.allocated(omegaIR)) allocate(omegaIR(nir,naerkind,nsizemax))
+   if (.not.allocated(gIR)) allocate(gIR(nir,naerkind,nsizemax))
+   if (.not.allocated(nsize)) allocate(nsize(naerkind,2))
+   if (.not.allocated(radiustab)) allocate(radiustab(naerkind,2,nsizemax))
+   if (.not.allocated(QREFvis)) allocate(QREFvis(naerkind,nsizemax))
+   if (.not.allocated(QREFir)) allocate(QREFir(naerkind,nsizemax))
+   if (.not.allocated(omegaREFvis)) allocate(omegaREFvis(naerkind,nsizemax))
+   if (.not.allocated(omegaREFir)) allocate(omegaREFir(naerkind,nsizemax))
 
   end subroutine ini_scatterers