Changeset 3976 for LMDZ6/trunk/libf/dyn3d

Timestamp:

Aug 6, 2021, 4:30:26 PM (3 years ago)

Author:

fhourdin

Message:

Options supplementaires pour iniacademic

File:

• LMDZ6/trunk/libf/dyn3d/iniacademic.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/dyn3d/iniacademic.F90

@@ -67 +67 @@
   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
@@ -117 +117 @@
   CALL inigeom
   CALL inifilr
+
+  ! Initialize pressure and mass field if read_start=.false.
+  IF (.NOT. read_start) THEN
+     ! 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(:)=101080.  
+     else
+        ! use reference surface pressure
+        ps(:)=preff
+     endif
+     ! ground geopotential
+     phis(:)=0.
+     CALL pression ( ip1jmp1, ap, bp, ps, p       )
+     if (pressure_exner) then
+       CALL exner_hyb( ip1jmp1, ps, p, pks, pk)
+     else
+       call exner_milieu(ip1jmp1,ps,p,pks,pk)
+     endif
+     CALL massdair(p,masse)
+  ENDIF
 
   if (llm == 1) then
@@ -167 +190 @@
      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
@@ -219 +245 @@
      ! 3. Initialize fields (if necessary)
      IF (.NOT. read_start) THEN
-        ! 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(:)=101080.  
-        else
-           ! use reference surface pressure
-           ps(:)=preff
-        endif
-        
-        ! ground geopotential
-        phis(:)=0.
-
-        CALL pression ( ip1jmp1, ap, bp, ps, p       )
-        if (pressure_exner) then
-          CALL exner_hyb( ip1jmp1, ps, p, pks, pk)
-        else
-          call exner_milieu(ip1jmp1,ps,p,pks,pk)
-        endif
-        CALL massdair(p,masse)
-
         ! bulk initialization of temperature
-        teta(:,:)=tetarappel(:,:)
+
+        IF (iflag_phys>10000) THEN
+        ! Particular case to impose a constant temperature T0=0.01*iflag_physx
+           teta(:,:)= 0.01*iflag_phys/(pk(:,:)/cpp)
+        ELSE
+           teta(:,:)=tetarappel(:,:)
+        ENDIF
 
         ! geopotential
         CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
+
+        DO l=1,llm
+          print*,'presnivs,play,l',presnivs(l),(pk(1,l)/cpp)**(1./kappa)*preff
+         !pks(ij) = (cpp/preff) * ps(ij)
+         !pk(ij,1) = .5*pks(ij)
+         ! pk = cpp * (p/preff)^kappa
+        ENDDO
 
         ! winds
@@ -292 +309 @@
         do l=1,llm
            do ij=iip2,ip1jm
-              teta(ij,l)=teta(ij,l)*(1.+0.005*ran1(idum))
+              teta(ij,l)=teta(ij,l)*(1.+tetanoise*ran1(idum))
            enddo
         enddo