Changeset 5084 for LMDZ6/trunk/libf/misc

Timestamp:

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

Author:

Laurent Fairhead

Message:

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

Location:

LMDZ6/trunk/libf/misc

Files:

• lmdz_netcdf.F90 (deleted)
• lmdz_xios.F90 (modified) (1 diff)
• write_field.F90 (modified) (6 diffs)
• wxios.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/misc/lmdz_xios.F90

@@ -12 +12 @@
 MODULE lmdz_xios
   !!!! Wrapper XIOS
-  !! => must be replaced later by official xios wrapper when available
+  !! => must be replaced latter by official xios wrapper when available
 
   LOGICAL,PARAMETER :: using_xios = .FALSE.

LMDZ6/trunk/libf/misc/write_field.F90

@@ +1 @@
+!
+! $Id$
+!
 module write_field
-  USE lmdz_netcdf, ONLY: nf_sync, nf90_put_var, nf_enddef, nf_def_dim, nf_unlimited, &
-      nf_clobber, nf90_format, nf_create, nf_def_var
-
-  implicit none
+implicit none
 
   integer, parameter :: MaxWriteField = 100
@@ -73 +73 @@
     subroutine WriteField_gen(name,Field,dimx,dimy,dimz)
     implicit none
+    include 'netcdf.inc'
       character(len=*) :: name
       integer :: dimx,dimy,dimz
@@ -101 +102 @@
       count(4)=1
 
-      status = nf90_put_var(FieldId(Index),FieldVarId(Index),Field,start,count)
+      status = NF_PUT_VARA_DOUBLE(FieldId(Index),FieldVarId(Index),start,count,Field)
       status = NF_SYNC(FieldId(Index))
       
@@ -108 +109 @@
     subroutine CreateNewField(name,dimx,dimy,dimz)
     implicit none
+    include 'netcdf.inc'  
       character(len=*) :: name
       integer :: dimx,dimy,dimz
@@ -124 +126 @@
       status = NF_DEF_DIM(FieldId(NbField),'Z',dimz,TabDim(3))
       status = NF_DEF_DIM(FieldId(NbField),'iter',NF_UNLIMITED,TabDim(4))
-      status = NF_DEF_VAR(FieldId(NbField),FieldName(NbField),NF90_FORMAT,4,TabDim,FieldVarId(NbField))
+      status = NF_DEF_VAR(FieldId(NbField),FieldName(NbField),NF_DOUBLE,4,TabDim,FieldVarId(NbField))
       status = NF_ENDDEF(FieldId(NbField))
 
     end subroutine CreateNewField
+    
+    
     
   subroutine write_field1D(name,Field)
@@ -281 +285 @@
                       //trim(int2str(pos+offset))      &    
                       //'," ---> ",g22.16," | ")'
-! d�pent de l'impl�mention, sur compaq, c'est necessaire
+! dépent de l'implémention, sur compaq, c'est necessaire
 !            Pos=Pos+ColumnSize
           endif

LMDZ6/trunk/libf/misc/wxios.F90

@@ -70 +70 @@
         reformaop = "average"
         
-        IF (op=="inst(X)") THEN
+        IF (op.EQ."inst(X)") THEN
             reformaop = "instant"
         END IF
         
-        IF (op=="once") THEN
+        IF (op.EQ."once") THEN
             reformaop = "once"
         END IF
         
-        IF (op=="t_max(X)") THEN
+        IF (op.EQ."t_max(X)") THEN
             reformaop = "maximum"
         END IF
         
-        IF (op=="t_min(X)") THEN
+        IF (op.EQ."t_min(X)") THEN
             reformaop = "minimum"
         END IF
@@ -604 +604 @@
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     SUBROUTINE wxios_add_field(fieldname, fieldgroup, fieldlongname, fieldunit)
-        USE lmdz_netcdf, only: nf90_fill_real
+        USE netcdf, only: nf90_fill_real
 
         IMPLICIT NONE
@@ -621 +621 @@
         def = nf90_fill_real
         
-        IF (fieldunit == " ") THEN
+        IF (fieldunit .EQ. " ") THEN
             newunit = "-"
         ELSE
@@ -666 +666 @@
         
         ! Ajout Abd pour NMC:
-        IF (fid<=6) THEN
+        IF (fid.LE.6) THEN
           axis_id="presnivs"
         ELSE
@@ -682 +682 @@
         
         !On selectionne le bon groupe de champs:
-        IF (fdim==2) THEN
+        IF (fdim.EQ.2) THEN
           CALL xios_get_handle("fields_2D", fieldgroup)
         ELSE
@@ -726 +726 @@
             CALL xios_set_attr(field, level=field_level, enabled=.TRUE.)
             
-            IF (fdim==2) THEN
+            IF (fdim.EQ.2) THEN
                 !Si c'est un champ 2D:
                 IF (prt_level >= 10) THEN