Changeset 1992 for LMDZ5/trunk/libf/phylmd/ini_wake.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (11 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/ini_wake.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/ini_wake.F) (1 diff)

LMDZ5/trunk/libf/phylmd/ini_wake.F90

@@ -1 @@
-!
+
 ! $Id$
-!
-      SUBROUTINE INI_WAKE(wape,fip,it_wape_prescr,
-     :     wape_prescr, fip_prescr, alp_bl_prescr, ale_bl_prescr)
-***************************************************************
-*                                                             *
-*        INI_WAKE : variables d'initialisation de la poche    *
-*                   froide, necessaires au declenchement      *
-*                   de la convection.                         *
-*                                                             *
-*                                                             *
-***************************************************************
-c Arguments
-c =========
-c Input
-c -----
-c   wape           : valeur de l'energie potentielle de la poche (WAPE)
-c                    dans l'etat initial
-c   fip            : valeur de la puissance incidente sur le front (FIP)
-c                    dans l'etat initial
-c Output
-c ------
-c   it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
-c            doit etre imposee.
-c   wape_prescr    : valeur prescrite de la WAPE.
-c   fip_prescr     : valeur prescrite de la FIP.
-c   ale_bl_prescr  : valeur prescrite de la Ale de PBL.
-c   alp_bl_prescr  : valeur prescrite de la Alp de PBL.
-c
-c Variables internes
-c ==================
-c   it = nbre de pas de temps lu
-c   w  = WAPE lue
-c   f  = FIP lue
-c   alebl  = Ale de PBL lue
-c   alpbl  = Alp de PBL lue
-c
-      include 'iniprint.h'
-cdeclarations
-      real ale_bl_prescr
-      real alp_bl_prescr
-      real it
 
-! FH A mettre si besoin dans physiq.def
-! FH : voir avec JYG
-      it=0.
-      w=4.
-      f=0.1
-      alebl=4.
-      alpbl=0.1
-c
-cCR: on rajoute ale et alp de la PBL precrits
-      open (99,file='ini_wake_param.data',form='formatted',
-     s      status='old',err=902)
-      read (99,*) it
-      read (99,*) w
-      read (99,*) f
-      read (99,*,end=901) alebl
-      read (99,*,end=901) alpbl
-901   close (99)
-902   continue
-c
-      write(lunout,*)' it,wape ',it,wape
-      it_wape_prescr = it
-      if (w .lt. 0) then
-         wape_prescr = wape
-         fip_prescr = fip
-      else
-         wape_prescr = w
-         fip_prescr = f
-      endif
-c
-      write(lunout,*)' alebl, alpbl ',alebl,alpbl
-      ale_bl_prescr=alebl
-      alp_bl_prescr=alpbl
-      print *,'Initialisation de la poche : WAPE, FIP imposees ='
-     $               ,wape_prescr, fip_prescr
-      print *, '                   pendant ',it_wape_prescr,' steps'
-c
-      print *,'Initialisation de la BL: ALP, ALE imposees ='
-     $               ,alp_bl_prescr, ale_bl_prescr
-      return
-      end
+SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, &
+    alp_bl_prescr, ale_bl_prescr)
+  ! **************************************************************
+  ! *
+  ! INI_WAKE : variables d'initialisation de la poche    *
+  ! froide, necessaires au declenchement      *
+  ! de la convection.                         *
+  ! *
+  ! *
+  ! **************************************************************
+  ! Arguments
+  ! =========
+  ! Input
+  ! -----
+  ! wape           : valeur de l'energie potentielle de la poche (WAPE)
+  ! dans l'etat initial
+  ! fip            : valeur de la puissance incidente sur le front (FIP)
+  ! dans l'etat initial
+  ! Output
+  ! ------
+  ! it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE
+  ! doit etre imposee.
+  ! wape_prescr    : valeur prescrite de la WAPE.
+  ! fip_prescr     : valeur prescrite de la FIP.
+  ! ale_bl_prescr  : valeur prescrite de la Ale de PBL.
+  ! alp_bl_prescr  : valeur prescrite de la Alp de PBL.
+
+  ! Variables internes
+  ! ==================
+  ! it = nbre de pas de temps lu
+  ! w  = WAPE lue
+  ! f  = FIP lue
+  ! alebl  = Ale de PBL lue
+  ! alpbl  = Alp de PBL lue
+
+  include 'iniprint.h'
+  ! declarations
+  REAL ale_bl_prescr
+  REAL alp_bl_prescr
+  REAL it
+
+  ! FH A mettre si besoin dans physiq.def
+  ! FH : voir avec JYG
+  it = 0.
+  w = 4.
+  f = 0.1
+  alebl = 4.
+  alpbl = 0.1
+
+  ! CR: on rajoute ale et alp de la PBL precrits
+  OPEN (99, FILE='ini_wake_param.data', FORM='formatted', STATUS='old', &
+    ERR=902)
+  READ (99, *) it
+  READ (99, *) w
+  READ (99, *) f
+  READ (99, *, END=901) alebl
+  READ (99, *, END=901) alpbl
+901 CLOSE (99)
+902 CONTINUE
+
+  WRITE (lunout, *) ' it,wape ', it, wape
+  it_wape_prescr = it
+  IF (w<0) THEN
+    wape_prescr = wape
+    fip_prescr = fip
+  ELSE
+    wape_prescr = w
+    fip_prescr = f
+  END IF
+
+  WRITE (lunout, *) ' alebl, alpbl ', alebl, alpbl
+  ale_bl_prescr = alebl
+  alp_bl_prescr = alpbl
+  PRINT *, 'Initialisation de la poche : WAPE, FIP imposees =', wape_prescr, &
+    fip_prescr
+  PRINT *, '                   pendant ', it_wape_prescr, ' steps'
+
+  PRINT *, 'Initialisation de la BL: ALP, ALE imposees =', alp_bl_prescr, &
+    ale_bl_prescr
+  RETURN
+END SUBROUTINE ini_wake