Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/interfoce_lim.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/interfoce_lim.F90 (modified) (17 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/interfoce_lim.F90

@@ -1 @@
-!
+
 ! $Header$
-!
+
 SUBROUTINE interfoce_lim(itime, dtime, jour, &
      knon, knindex, &
@@ -16 +16 @@
 ! Cette routine sert d'interface entre le modele atmospherique et un fichier
 ! de conditions aux limites
-!
+
 ! L. Fairhead 02/2000
-!
+
 ! input:
 !   itime        numero du pas de temps courant
@@ -28 +28 @@
 !   klon         taille de la grille
 !   debut        logical: 1er appel a la physique (initialisation)
-!
+
 ! output:
 !   lmt_sst_p      SST lues dans le fichier de CL
 !   pctsrf_new-p   sous-maille fractionnelle
-!
-
 
 ! Parametres d'entree
@@ -83 +81 @@
   REAL, DIMENSION(klon_glo)                 :: sst_lu
   REAL, DIMENSION(klon_glo)                 :: nat_lu
-!
+
 ! Fin declaration
 !****************************************************************************************
@@ -89 +87 @@
 !****************************************************************************************
 ! Start calculation
-!
+
 !****************************************************************************************
   IF (debut .AND. .NOT. ALLOCATED(sst_lu_p)) THEN
@@ -105 +103 @@
 !****************************************************************************************
 ! Ouverture et lecture du fichier pour le master process si c'est le bon moment
-!
+
 !****************************************************************************************
 ! Tester d'abord si c'est le moment de lire le fichier
@@ -128 +126 @@
 
         IF (newlmt) THEN
-           !
+
            ! Fraction "ocean" 
-           !
+
            ierr = nf90_inq_varid(nid, 'FOCE', nvarid)
            IF (ierr /= nf90_noerr) THEN
@@ -141 +139 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           !
+
            ! Fraction "glace de mer" 
-           !
+
            ierr = nf90_inq_varid(nid, 'FSIC', nvarid)
            IF (ierr /= nf90_noerr) THEN
@@ -154 +152 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           !
+
            ! Fraction "terre" 
-           !
+
            ierr = nf90_inq_varid(nid, 'FTER', nvarid)
            IF (ierr /= nf90_noerr) THEN
@@ -167 +165 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           !
+
            ! Fraction "glacier terre" 
-           !
+
            ierr = nf90_inq_varid(nid, 'FLIC', nvarid)
            IF (ierr /= nf90_noerr) THEN
@@ -180 +178 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           !
+
         ELSE  ! on en est toujours a rnatur
-           ! 
+
            ierr = nf90_inq_varid(nid, 'NAT', nvarid)
            IF (ierr /= nf90_noerr) THEN
@@ -193 +191 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-!
+
 ! Remplissage des fractions de surface
 ! nat = 0, 1, 2, 3 pour ocean, terre, glacier, seaice
-! 
+
            pct_tmp = 0.0
            DO ii = 1, klon_glo
@@ -202 +200 @@
            ENDDO
 
-!
 !  On se retrouve avec ocean en 1 et terre en 2 alors qu'on veut le contraire
-!
+
            pctsrf_new = pct_tmp
            pctsrf_new (:,2)= pct_tmp (:,1)
@@ -210 +207 @@
            pct_tmp = pctsrf_new 
         ENDIF ! fin test sur newlmt
-!
+
 ! Lecture SST
-!
+
         ierr = nf90_inq_varid(nid, 'SST', nvarid)
         IF (ierr /= nf90_noerr) THEN
@@ -226 +223 @@
 !****************************************************************************************
 ! Fin de lecture, fermeture de fichier
-!
+
 !****************************************************************************************
         ierr = nf90_close(nid)
@@ -237 +234 @@
 !****************************************************************************************
 ! Distribue les variables sur tous les processus
-!
+
 !****************************************************************************************
      CALL Scatter(sst_lu,sst_lu_p)
@@ -248 +245 @@
 !****************************************************************************************
 ! Recopie des variables dans les champs de sortie
-!
+
 !****************************************************************************************
   lmt_sst_p = 999999999.