Changeset 4013 for LMDZ6/branches/Ocean_skin/libf/dyn3dmem/iniacademic_loc.F90

Timestamp:

Nov 19, 2021, 4:58:59 PM (3 years ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/dyn3dmem/iniacademic_loc.F90 (modified) (5 diffs)

LMDZ6/branches/Ocean_skin/libf/dyn3dmem/iniacademic_loc.F90

@@ -73 +73 @@
   LOGICAL ok_geost             ! Initialisation vent geost. ou nul
   LOGICAL ok_pv                ! Polar Vortex
-  REAL phi_pv,dphi_pv,gam_pv   ! Constantes pour polar vortex 
+  REAL phi_pv,dphi_pv,gam_pv,tetanoise   ! Constantes pour polar vortex 
 
   real zz,ran1
@@ -122 +122 @@
   CALL inigeom
   CALL inifilr
+
+  ! Initialize pressure and mass field if read_start=.false.
+  IF (.NOT. read_start) THEN
+    ! allocate global fields:
+!    allocate(vcov_glo(ip1jm,llm))
+    allocate(ucov_glo(ip1jmp1,llm))
+    allocate(teta_glo(ip1jmp1,llm))
+    allocate(ps_glo(ip1jmp1))
+    allocate(masse_glo(ip1jmp1,llm))
+    allocate(phis_glo(ip1jmp1))
+
+     ! surface pressure
+     if (iflag_phys>2) then
+        ! specific value for CMIP5 aqua/terra planets
+        ! "Specify the initial dry mass to be equivalent to
+        !  a global mean surface pressure (101325 minus 245) Pa."
+        ps_glo(:)=101080.  
+     else
+        ! use reference surface pressure
+        ps_glo(:)=preff
+     endif
+     
+     ! ground geopotential
+     phis_glo(:)=0.
+
+     CALL pression ( ip1jmp1, ap, bp, ps_glo, p       )
+     if (pressure_exner) then
+       CALL exner_hyb( ip1jmp1, ps_glo, p, pks, pk )
+     else
+       call exner_milieu(ip1jmp1,ps_glo,p,pks,pk)
+     endif
+     CALL massdair(p,masse_glo)
+  ENDIF
+
 
   if (llm == 1) then
@@ -172 +206 @@
      gam_pv=4.              ! -dT/dz vortex (in K/km)
      CALL getin('gam_pv',gam_pv)
+     tetanoise=0.005
+     CALL getin('tetanoise',tetanoise)
 
      ! 2. Initialize fields towards which to relax
@@ -224 +260 @@
      ! 3. Initialize fields (if necessary)
      IF (.NOT. read_start) THEN
-       ! allocate global fields:
-!       allocate(vcov_glo(ip1jm,llm))
-       allocate(ucov_glo(ip1jmp1,llm))
-       allocate(teta_glo(ip1jmp1,llm))
-       allocate(ps_glo(ip1jmp1))
-       allocate(masse_glo(ip1jmp1,llm))
-       allocate(phis_glo(ip1jmp1))
-
-        ! surface pressure
-        if (iflag_phys>2) then
-           ! specific value for CMIP5 aqua/terra planets
-           ! "Specify the initial dry mass to be equivalent to
-           !  a global mean surface pressure (101325 minus 245) Pa."
-           ps_glo(:)=101080.  
-        else
-           ! use reference surface pressure
-           ps_glo(:)=preff
-        endif
-        
-        ! ground geopotential
-        phis_glo(:)=0.
-
-        CALL pression ( ip1jmp1, ap, bp, ps_glo, p       )
-        if (pressure_exner) then
-          CALL exner_hyb( ip1jmp1, ps_glo, p, pks, pk )
-        else
-          call exner_milieu(ip1jmp1,ps_glo,p,pks,pk)
-        endif
-        CALL massdair(p,masse_glo)
-
         ! bulk initialization of temperature
-        teta_glo(:,:)=tetarappel(:,:)
-
+        IF (iflag_phys>10000) THEN
+        ! Particular case to impose a constant temperature T0=0.01*iflag_phys
+           teta_glo(:,:)= 0.01*iflag_phys/(pk(:,:)/cpp)
+        ELSE
+           teta_glo(:,:)=tetarappel(:,:)
+        ENDIF
         ! geopotential
         CALL geopot(ip1jmp1,teta_glo,pk,pks,phis_glo,phi)
@@ -306 +316 @@
         do l=1,llm
            do ij=iip2,ip1jm
-              teta_glo(ij,l)=teta_glo(ij,l)*(1.+0.005*ran1(idum))
+              teta_glo(ij,l)=teta_glo(ij,l)*(1.+tetanoise*ran1(idum))
            enddo
         enddo