Changeset 5084 for LMDZ6/trunk/libf/phylmd/interfoce_lim.F90

Timestamp:

Jul 19, 2024, 6:40:44 PM (4 months ago)

Author:

Laurent Fairhead

Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

File:

• LMDZ6/trunk/libf/phylmd/interfoce_lim.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/interfoce_lim.F90

@@ -10 +10 @@
   USE mod_phys_lmdz_para
   USE indice_sol_mod
-  USE lmdz_netcdf, ONLY: nf90_get_var,nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite
   
   IMPLICIT NONE
+  
+  INCLUDE "netcdf.inc"
 
 ! Cette routine sert d'interface entre le modele atmospherique et un fichier
@@ -115 +116 @@
         fich = TRIM(fich)
         ierr = NF_OPEN (fich, NF_NOWRITE,nid)
-        IF (ierr/=NF_NOERR) THEN
+        IF (ierr.NE.NF_NOERR) THEN
            abort_message = 'Pb d''ouverture du fichier de conditions aux limites'
            CALL abort_physic(modname,abort_message,1)
@@ -136 +137 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           ierr = nf90_get_var(nid,nvarid,pct_tmp(:,is_oce),start,epais)
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais,pct_tmp(1,is_oce))
+#else
+           ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais,pct_tmp(1,is_oce))
+#endif
            IF (ierr /= NF_NOERR) THEN
               abort_message = 'Lecture echouee pour <FOCE>'
@@ -149 +154 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           ierr = nf90_get_var(nid,nvarid,pct_tmp(:,is_sic),start,epais)
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais,pct_tmp(1,is_sic))
+#else
+           ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais,pct_tmp(1,is_sic))
+#endif
            IF (ierr /= NF_NOERR) THEN
               abort_message = 'Lecture echouee pour <FSIC>'
@@ -162 +171 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           ierr = nf90_get_var(nid,nvarid,pct_tmp(:,is_ter),start,epais)
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais,pct_tmp(1,is_ter))
+#else
+           ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais,pct_tmp(1,is_ter))
+#endif
            IF (ierr /= NF_NOERR) THEN
               abort_message = 'Lecture echouee pour <FTER>'
@@ -175 +188 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           ierr = nf90_get_var(nid,nvarid,pct_tmp(:,is_lic),start,epais)
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais,pct_tmp(1,is_lic))
+#else
+           ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais,pct_tmp(1,is_lic))
+#endif
            IF (ierr /= NF_NOERR) THEN
               abort_message = 'Lecture echouee pour <FLIC>'
@@ -188 +205 @@
               CALL abort_physic(modname,abort_message,1)
            ENDIF
-           ierr = nf90_get_var(nid,nvarid,nat_lu,start,epais)
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais, nat_lu)
+#else
+           ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais, nat_lu)
+#endif
            IF (ierr /= NF_NOERR) THEN
               abort_message = 'Lecture echouee pour <NAT>'
@@ -218 +239 @@
            CALL abort_physic(modname,abort_message,1)
         ENDIF
-        ierr = nf90_get_var(nid,nvarid,sst_lu,start,epais)
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,epais, sst_lu)
+#else
+        ierr = NF_GET_VARA_REAL(nid,nvarid,start,epais, sst_lu)
+#endif
         IF (ierr /= NF_NOERR) THEN
            abort_message = 'Lecture echouee pour <SST>'