Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/guide_loc_mod.F90

Timestamp:

Jul 23, 2024, 3:29:36 PM (4 months ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/guide_loc_mod.F90 (modified) (10 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/guide_loc_mod.F90

@@ -92 +92 @@
     CALL ini_getparam("nudging_parameters_out.txt")
 ! Variables guidees
-    CALL getpar('guide_u',.true.,guide_u,'guidage de u')
-    CALL getpar('guide_v',.true.,guide_v,'guidage de v')
-    CALL getpar('guide_T',.true.,guide_T,'guidage de T')
-    CALL getpar('guide_P',.true.,guide_P,'guidage de P')
-    CALL getpar('guide_Q',.true.,guide_Q,'guidage de Q')
-    CALL getpar('guide_hr',.true.,guide_hr,'guidage de Q par H.R')
-    CALL getpar('guide_teta',.false.,guide_teta,'guidage de T par Teta')
-
-    CALL getpar('guide_add',.false.,guide_add,'foréage constant?')
-    CALL getpar('guide_zon',.false.,guide_zon,'guidage moy zonale')
+    CALL getpar('guide_u',.TRUE.,guide_u,'guidage de u')
+    CALL getpar('guide_v',.TRUE.,guide_v,'guidage de v')
+    CALL getpar('guide_T',.TRUE.,guide_T,'guidage de T')
+    CALL getpar('guide_P',.TRUE.,guide_P,'guidage de P')
+    CALL getpar('guide_Q',.TRUE.,guide_Q,'guidage de Q')
+    CALL getpar('guide_hr',.TRUE.,guide_hr,'guidage de Q par H.R')
+    CALL getpar('guide_teta',.FALSE.,guide_teta,'guidage de T par Teta')
+
+    CALL getpar('guide_add',.FALSE.,guide_add,'foréage constant?')
+    CALL getpar('guide_zon',.FALSE.,guide_zon,'guidage moy zonale')
     if (guide_zon .and. abs(grossismx - 1.) > 0.01) &
          CALL abort_gcm("guide_init", &
@@ -117 +117 @@
     CALL getpar('tau_min_P',0.02,tau_min_P,'Cste de rappel min, P')
     CALL getpar('tau_max_P', 10.,tau_max_P,'Cste de rappel max, P')
-    CALL getpar('gamma4',.false.,gamma4,'Zone sans rappel elargie')
-    CALL getpar('guide_BL',.true.,guide_BL,'guidage dans C.Lim')
+    CALL getpar('gamma4',.FALSE.,gamma4,'Zone sans rappel elargie')
+    CALL getpar('guide_BL',.TRUE.,guide_BL,'guidage dans C.Lim')
     CALL getpar('plim_guide_BL',85000.,plim_guide_BL,'BL top presnivs value')
 
 ! Sauvegarde du forçage
-    CALL getpar('guide_sav',.false.,guide_sav,'sauvegarde guidage')
+    CALL getpar('guide_sav',.FALSE.,guide_sav,'sauvegarde guidage')
     CALL getpar('iguide_sav',4,iguide_sav,'freq. sauvegarde guidage')
     ! frequences f>0: fx/jour; f<0: tous les f jours; f=0: 1 seule fois.
@@ -134 +134 @@
 
 ! Guidage regional seulement (sinon constant ou suivant le zoom)
-    CALL getpar('guide_reg',.false.,guide_reg,'guidage regional')
+    CALL getpar('guide_reg',.FALSE.,guide_reg,'guidage regional')
     CALL getpar('lat_min_g',-90.,lat_min_g,'Latitude mini guidage ')
     CALL getpar('lat_max_g', 90.,lat_max_g,'Latitude maxi guidage ')
@@ -154 +154 @@
     CALL getpar('guide_plevs',0,guide_plevs,'niveaux pression fichiers guidage')
     ! Pour compatibilite avec ancienne version avec guide_modele
-    CALL getpar('guide_modele',.false.,guide_modele,'niveaux pression ap+bp*psol')
+    CALL getpar('guide_modele',.FALSE.,guide_modele,'niveaux pression ap+bp*psol')
     IF (guide_modele) THEN
         guide_plevs=1
     ENDIF
 !FC
-    CALL getpar('convert_Pa',.true.,convert_Pa,'Convert Pressure levels in Pa')
+    CALL getpar('convert_Pa',.TRUE.,convert_Pa,'Convert Pressure levels in Pa')
     ! Fin raccord
-    CALL getpar('ini_anal',.false.,ini_anal,'Etat initial = analyse')
-    CALL getpar('guide_invertp',.true.,invert_p,'niveaux p inverses')
-    CALL getpar('guide_inverty',.true.,invert_y,'inversion N-S')
-    CALL getpar('guide_2D',.false.,guide_2D,'fichier guidage lat-P')
+    CALL getpar('ini_anal',.FALSE.,ini_anal,'Etat initial = analyse')
+    CALL getpar('guide_invertp',.TRUE.,invert_p,'niveaux p inverses')
+    CALL getpar('guide_inverty',.TRUE.,invert_y,'inversion N-S')
+    CALL getpar('guide_2D',.FALSE.,guide_2D,'fichier guidage lat-P')
 
     CALL fin_getparam
@@ -2366 +2366 @@
   
 !===========================================================================
-  subroutine correctbid(iim,nl,x)
+  SUBROUTINE correctbid(iim,nl,x)
     integer iim,nl
     real x(iim+1,nl)
@@ -2376 +2376 @@
             if(abs(x(i,l))>1.e10) then
                zz=0.5*(x(i-1,l)+x(i+1,l))
-              print*,'correction ',i,l,x(i,l),zz
+              PRINT*,'correction ',i,l,x(i,l),zz
                x(i,l)=zz
             endif
@@ -2382 +2382 @@
      enddo
      return
-  end subroutine correctbid
+  END SUBROUTINE  correctbid
 
 
@@ -2389 +2389 @@
 !====================================================================
 
-subroutine dump2du(var,varname)
+SUBROUTINE dump2du(var,varname)
 use parallel_lmdz
 use mod_hallo
@@ -2416 +2416 @@
 
     return
-    end subroutine dump2du
+    END SUBROUTINE  dump2du
 
 !====================================================================
 ! Ascii debug output. Could be reactivated
 !====================================================================
-subroutine dumpall
+SUBROUTINE dumpall
      implicit none
      include "dimensions.h"
@@ -2431 +2431 @@
      CALL dump2du(ugui1(ijb_u:ije_u,1)*sqrt(unscu2(ijb_u:ije_u)),'  ugui1 couche 1')
      return
-end subroutine dumpall
+END SUBROUTINE  dumpall
 
 !===========================================================================