Changeset 752

Timestamp:

10/06/09 15:20:26 (15 years ago)

Author:

bellier

Message:

Added the possibility to choose the kind of NETCDF variables (R4/R8) by calling histdef

File:

• IOIPSL/trunk/src/histcom.f90 (modified) (6 diffs)

IOIPSL/trunk/src/histcom.f90

@@ -76 +76 @@
  &                     nb_hax_max=5,nb_zax_max=10,nbopp_max=10
   REAL,PARAMETER :: missing_val=nf90_fill_real
+  INTEGER,PARAMETER,PUBLIC :: &
+ &  hist_r4=nf90_real4, hist_r8=nf90_real8
 !-
   INTEGER :: bufftmp_max(nb_files_max) = 1
@@ -116 +118 @@
 ! Informations on each variable
 !-
-  INTEGER,DIMENSION(nb_files_max,nb_var_max),SAVE :: &
- &  nbopp
+  INTEGER,DIMENSION(nb_files_max,nb_var_max),SAVE :: v_typ
   CHARACTER(LEN=20),DIMENSION(nb_files_max,nb_var_max),SAVE :: &
  &  name,unit_name
@@ -142 +143 @@
  &  freq_opp,freq_wrt
   INTEGER,DIMENSION(nb_files_max,nb_var_max),SAVE :: &
- &  last_opp,last_wrt,last_opp_chk,last_wrt_chk,nb_opp,nb_wrt,point
+ &  nbopp,last_opp,last_wrt,last_opp_chk,last_wrt_chk,nb_opp,nb_wrt,point
 !-
 ! Book keeping for the buffers
@@ -1104 +1105 @@
 !- pzid     : ID of the vertical axis to use. It has to have
 !-            the size of the zoom.
-!- pnbbyt   : Number of bytes on which to store in netCDF (Not opp.)
+!- pnbbyt   : External netCDF type (hist_r4/hist_r8)
 !- popp     : Operation to be performed. The following options
 !-            exist today :
@@ -1197 +1198 @@
     sopps(pfileid,iv,nbopp(pfileid,iv)) = 'ident'
     scal(pfileid,iv,nbopp(pfileid,iv)) = missing_val
+  ENDIF
+!-
+! 1.3 External type of the variable
+!-
+  IF (pnbbyt == hist_r8) THEN
+    v_typ(pfileid,iv) = hist_r8
+  ELSE
+    v_typ(pfileid,iv) = hist_r4
   ENDIF
 !-
@@ -1613 +1622 @@
       ENDIF
 !-
-      iret = NF90_DEF_VAR (ncid,TRIM(name(pfileid,iv)),NF90_FLOAT, &
- &                         dims(1:ABS(ndim)),ncvarid)
+      iret = NF90_DEF_VAR (ncid,TRIM(name(pfileid,iv)), &
+ &                         v_typ(pfileid,iv),dims(1:ABS(ndim)),ncvarid)
 !-
       ncvar_ids(pfileid,iv) = ncvarid