Changeset 614

Timestamp:

Mar 30, 2005, 4:56:22 PM (20 years ago)

Author:

Laurent Fairhead

Message:

Rajout d'un cnat 'mixte'
Rajout d'un axe y bidon pour le lecture avec grads
LF

Location:

BOL/trunk/Class_Reg

Files:

• config.def (modified) (1 diff)
• geo2reg.F90 (modified) (7 diffs)

BOL/trunk/Class_Reg/config.def

@@ -1 +1 @@
 # type de classification (w500 pour l'instant)
 bintype=w500
-# sur quelle surface? (glob, ocea, land)
-cnat=ocea
+# sur quelle surface? (glob, ocea, mixt, land)
+cnat=land
 # limites des latitudes a traiter -lat0 < donnees < lat0
 lat0=30.0

BOL/trunk/Class_Reg/geo2reg.F90

@@ -62 +62 @@
 ! netcdf stuff
   real         :: time
-  integer, dimension(:), allocatable :: var_dim
+  real, dimension(:), allocatable :: var_dim
   INTEGER      :: ierr, nvdims, lunout, in_file_id, ierrs, ndims, natts
   integer      :: lonid, latid, in_time_id, varid, out_file_id, ngatts, levid
@@ -140 +140 @@
 !---------------------------------------------------------------------
 
-  if (cnat.ne.'ocea'.and.cnat.ne.'land'.and.cnat.ne.'glob') then
+  if (cnat.ne.'ocea'.and.cnat.ne.'land'.and.cnat.ne.'glob'.and.cnat.ne.'mixt')&
+& then
     write(*,*) 'erreur cnat ',cnat
     stop
@@ -273 +274 @@
     do j = 1, jm
       ok_msk(i,j) = .FALSE.
-      if (cnat.eq.'ocea' .and. msk(i,j).le.0.5) ok_msk(i,j) = .TRUE.
-      if (cnat.eq.'land' .and. msk(i,j).gt.0.5) ok_msk(i,j) = .TRUE.
+      if (cnat.eq.'ocea' .and. msk(i,j).le.0.01) ok_msk(i,j) = .TRUE.
+      if (cnat.eq.'mixt' .and. msk(i,j).gt.0.01.and. msk(i,j).le.0.99) &
+  &                                             ok_msk(i,j) = .TRUE.
+      if (cnat.eq.'land' .and. msk(i,j).gt.0.99) ok_msk(i,j) = .TRUE.
       if (cnat.eq.'glob') ok_msk(i,j) = .TRUE.
     enddo
@@ -315 +318 @@
   nvdims = 1
   vdims(1) = nbin_id
-  ierr = NF_DEF_VAR(out_file_id, 'nbin', NF_INT, nvdims, vdims, var_bin_id)
-  ierrs = ierrs + ierr
-  ierr = NF_PUT_ATT_TEXT(out_file_id, var_bin_id, 'units', 1, '-')
+  ierr = NF_DEF_VAR(out_file_id, 'nbin', NF_FLOAT, nvdims, vdims, var_bin_id)
+  ierrs = ierrs + ierr
+  ierr = NF_PUT_ATT_TEXT(out_file_id, var_bin_id, 'units', 5, 'hPa/j')
   ierrs = ierrs + ierr
   ierr = NF_PUT_ATT_TEXT(out_file_id, var_bin_id, 'long_name', 15, &
@@ -326 +329 @@
     stop
   endif
+!
+! y bidon pour grads
+  ierr = NF_DEF_DIM(out_file_id, 'lat', 1, latid)
+  nvdims = 1
+  vdims(1) = latid
+  ierr = NF_DEF_VAR(out_file_id, 'lat', NF_FLOAT, nvdims, vdims,latid)
+  ierr = NF_PUT_ATT_TEXT(out_file_id, latid, 'units', 13, 'degrees_north')
+
 !
 ! presnivs
@@ -569 +580 @@
       endif ! lat
     enddo ! j
-    write(lunout,*)'nbre de points ou w500 < ',min_bin,' = ', nbr_inf
-    write(lunout,*)'nbre de points ou w500 > ',max_bin,' = ', nbr_sup 
+    if (nbr_inf /=0) &
+    & write(lunout,*)'nbre de points ou w500 < ',min_bin,' = ', nbr_inf
+    if (nbr_sup /=0) &
+    & write(lunout,*)'nbre de points ou w500 > ',max_bin,' = ', nbr_sup 
   enddo ! m
 
@@ -620 +633 @@
   deallocate(var_dim)
 !
+! ecriture y bidon
+!
+  ierr = NF_PUT_VAR1_REAL(out_file_id, latid, 1, 0.)
+  
+!
 ! ecriture niveaux verticaux
 !