Changeset 2534 for trunk/LMDZ.VENUS/libf/phyvenus/turb_mod.F90

Timestamp:

Jun 9, 2021, 9:38:56 PM (3 years ago)

Author:

emillour

Message:

Venus GCM:
Some cosmetics in routines, along with adding the possibility to read/write q2 from/to startphy.nc/restartsphy.nc (not really used, but in preparation of future commits).
VB+EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/turb_mod.F90 (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/turb_mod.F90

@@ -4 +4 @@
   REAL,SAVE,ALLOCATABLE :: q2(:,:)    ! Turbulent Kinetic Energy
   REAL,SAVE,ALLOCATABLE :: sens(:)
-  REAL,allocatable,SAVE :: l0(:)
+  REAL,allocatable,SAVE :: l0(:)      ! Characteristic mixing length
   REAL,SAVE,ALLOCATABLE :: yustar(:)
   REAL,SAVE,ALLOCATABLE :: wstar(:)
@@ -15 +15 @@
       ! mostly for LES use. default is FALSE (for GCM and mesoscale)
 
-!contains
+contains
 
-!  subroutine ini_turb_mod(ngrid,nlayer)
+  subroutine ini_turb_mod(ngrid,nlayer)
 
-!  implicit none
-!  integer,intent(in) :: ngrid ! number of atmospheric columns
-!  integer,intent(in) :: nlayer ! number of atmospheric layers
+  implicit none
+  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: nlayer ! number of atmospheric layers
 
-!    allocate(q2(ngrid,nlayer+1))
-!    allocate(l0(ngrid))
-!    allocate(wstar(ngrid))
-!   allocate(ustar(ngrid))
-!   allocate(tstar(ngrid))
-!   allocate(hfmax_th(ngrid))
-!   allocate(zmax_th(ngrid))
+    allocate(q2(ngrid,nlayer+1))
+    allocate(sens(ngrid))
+    allocate(l0(ngrid))
+    allocate(wstar(ngrid))
+    allocate(yustar(ngrid))
+    allocate(tstar(ngrid))
+    allocate(hfmax_th(ngrid))
+    allocate(zmax_th(ngrid))
 !   allocate(sensibFlux(ngrid))
 
-!  end subroutine ini_turb_mod
+  end subroutine ini_turb_mod
+
+  subroutine end_turb_mod
+
+  implicit none
+
+    if (allocated(q2)) deallocate(q2)
+    if (allocated(sens)) deallocate(sens)
+    if (allocated(l0)) deallocate(l0)
+    if (allocated(wstar)) deallocate(wstar)
+    if (allocated(yustar)) deallocate(yustar)
+    if (allocated(tstar)) deallocate(tstar)
+    if (allocated(hfmax_th)) deallocate(hfmax_th)
+    if (allocated(zmax_th)) deallocate(zmax_th)
+!    if (allocated(sensibFlux)) deallocate(sensibFlux)
+
+  end subroutine end_turb_mod
+
 
 end module turb_mod