Changeset 28 for trunk

Timestamp:

07/27/07 15:20:37 (17 years ago)

Author:

pinsard

Message:

forcagequimarche now run on my account in aedon

File:

• trunk/forcagequimarche.pro (modified) (9 diffs)

trunk/forcagequimarche.pro

@@ -1 @@
-PRO forcagequimarche,iyear,ian 
-@init2
-@initorca2_bab
-
-
-
+; +
+;
+; @file_comments
+; calcul d'une matrice 3D de conductivité sigma=f(T,S) en fonction de T et S
+; de ORCA, sur la grille T
+;
+; @categories
+; orca grid
+;
+; @param orcares {in}{required}{type=string}
+; code of ORCA grid to use for output results
+; must be 'ORCA2'
+;
+; @param iyear {in}{required}{type=integer}
+; real year
+;
+; @param ian {in}{required}{type=integer}
+; simulation year
+;
+; @keyword DRAKKAR_EXP {type=string}
+; code for Drakkar experiment
+; only used when orcares = ORCA025
+; must be G42 ++ G70
+;
+; @restrictions
+;  - Requires SAXO tools
+;
+; @todo
+; validation of ORCA2 processing (now some nan in output files !)
+; add ORCA025
+; provide tools to plot output files
+; produce a NetCDF GDT or CF compliant
+;
+; @pre
+; see geomag_profile.sh
+; be sure to have ++ in the directory defined in ${GEOMAG_ID}/
+;
+; @post
+; see geomag_profile.sh
+; ++ in the directory defined in ${GEOMAG_OD}/
+;
+; @examples
+; for ORCA2
+; IDL> .run forcagequimarche
+; IDL> forcagequimarche, 'ORCA2', '1993', '01'
+;
+; @history
+; fplod 2007-07-27T10:29:18Z cerbere.locean-ipsl.upmc.fr (Linux)
+; add orcares parameter
+; add optional keyword drakkar_exp
+; add header for idldoc
+; add check of io directories
+; replace io directories user dependant to $GEOMAG_ID and $GEOMAG_OD
+; replace call of divfred by call of div
+;
+; @version
+; $Id$
+;
+;-
+;
+PRO forcagequimarche,orcares,iyear,ian,DRAKKAR_EXP = drakkar_exp
+;
+;----
+; check input parameters
+;----
+;
+; check orcares definition
+;
+  CASE orcares OF
+     'ORCA2': BEGIN
+                 msg = 'iii : valid orcares parameter = ' + orcares
+                 ras = report(msg)
+                 filename_oce = 'meshmask_bab.nc'
+                 IF keyword_set(DRAKKAR_EXP) THEN BEGIN
+                    msg = 'www : unused DRAKKAR_EXP keyword = ' + drakkar_exp
+                    ras = report(msg)
+                 END
+              END
+      ELSE  : BEGIN
+                 msg = 'eee : invalid orcares parameter = ' + orcares
+                 ras = report(msg)
+                 msg = 'eee : orcares must be ORCA2 or ORCA025++'
+                 ras = report(msg)
+                 RETURN
+              END
+  ENDCASE
+
+; ++ check iyear
+; ++ check ian
+
+;++@init2
+; init grid
+initocemeshmask, orcares, DRAKKAR_EXP = drakkar_exp
 
 @common
-
-
-
-;ian='01'
-;iyear='1993'
-
+;
+; check for input files
+;
+; test if ${GEOMAG_ID} defined
+  geomag_id_env=GETENV('GEOMAG_ID')
+  CASE geomag_id_env OF
+     ''  :  BEGIN
+              msg = 'eee : ${GEOMAG_ID} is not defined'
+              ras = report(msg)
+              RETURN
+            END
+     ELSE: BEGIN
+             msg = 'iii : ${GEOMAG_ID} is ' + geomag_id_env
+             ras = report(msg)
+           END
+  ENDCASE
+;
+  iodirin = isadirectory(geomag_id_env)
+;
+; existence and protection of ${GEOMAG_ID}
+  IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN
+     msg = 'eee : the directory' + iodirin  + ' is not accessible.'
+     ras = report(msg)
+     RETURN
+  ENDIF
+;
+; test if ${GEOMAG_OD} defined
+  geomag_od_env=GETENV('GEOMAG_OD')
+  CASE geomag_od_env OF
+     '' : BEGIN
+             msg = 'eee : ${GEOMAG_OD} is not defined'
+             ras = report(msg)
+             RETURN
+          END
+     ELSE: BEGIN
+             msg = 'iii : ${GEOMAG_OD} is ' + geomag_od_env
+             ras = report(msg)
+           END
+  ENDCASE
+;
+; check if output data will be possible
+  iodirout = isadirectory(geomag_od_env)
+;
+; existence and protection
+  IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN
+     msg = 'eee : the directory' + iodirout  + ' was not found.'
+     ras = report(msg)
+     RETURN
+  ENDIF
+;
 e_exp='ESS'
-rep_fred='/usr/work/sur/fvi/OPA/geomag/'
 key_portrait = 0
-; stockage des fichiers brut 
-  ioDATA='/usr/work/sur/fvi/OPA/ORCA2/'
-  file_U=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_U.nc' 
-  print, file_U
-  file_V=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_V.nc' 
-  print, file_V
-  file_T=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_T.nc' 
-  print, file_T
-  file_Sed= rep_fred+'cond_sed_ORCA2.nc'
-  file_Br= rep_fred+'Br_ORCA2.nc'
-
-
+; stockage des fichiers brut
+  ioDATA=geomag_id_env
+
+  CASE orcares OF
+     'ORCA2': BEGIN
+  file_U=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_U.nc'
+  file_V=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_V.nc'
+  file_T=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_T.nc'
+  file_Sed= geomag_id_env+'cond_sed_ORCA2.nc'
+  file_Br= geomag_id_env+'Br_ORCA2.nc'
+     END
+  ENDCASE
+;
 ; title
      t_exp= e_exp
-;     t_exp0= e_exp0
      t_bt  = 'bar_transp'
-ioORLN2 = '/usr/work/sur/fvi/OPA/ORCA2'
+ioORLN2 = geomag_id_env
+;
 ;facteur d'echelle vertical  for partial steps
 e3v3d=read_ncdf('e3v_ps',0,/timestep,iodir=ioORLN2,/nostruct,/tout,filename='meshmask_bab.nc')
@@ -39 +170 @@
 ;BR=read_ncdf('Br',0,/timestep,/nostruct,/tout,filename=file_Br)
 BR=read_ncdf('Br',0,/timestep,/nostruct,/tout,filename=file_Br,/cont_nofill)
-
-
-
-
+;
 ; vertical integration:
-e3t3d=make_array(jpi,jpj,jpk)   
+e3t3d=make_array(jpi,jpj,jpk)
         for k=0, jpk-1 do begin              &$
-          for j=0,jpj-1 do begin              &$      
+          for j=0,jpj-1 do begin              &$
             for i=0,jpi-1 do begin             &$
               e3t3d(i,j,k) = e3t(k)    &$
@@ -53 +181 @@
         endfor
 jpt = 73
-
-;vud = make_array(jpi,jpj,jpt)   
+;
+;vud = make_array(jpi,jpj,jpt)
 ;vvd = make_array(jpi, jpj, jpt)
 divBustar = make_array(jpi, jpj, jpt)
 diver3=replicate(0.,182,149,1,73)
 sigma3=replicate(0.,182,149,1,73)
-
-
-FOR jt = 0, jpt-1 DO BEGIN &$   
-
+;
+FOR jt = 0, jpt-1 DO BEGIN &$
+;
 ; ouverture des fichiers et stockage en memoire partial steps
-  vu=read_ncdf('vozocrtx',jt,jt, /timestep, iodir=ioDATA,/nostruct,/TOUT,filename=file_U)  &$   
-;stop
-  vv=read_ncdf('vomecrty',jt,jt, /timestep,iodir=ioDATA,/nostruct,/TOUT,filename=file_V)  &$   
+  vu=read_ncdf('vozocrtx',jt,jt, /timestep, iodir=ioDATA,/nostruct,/TOUT,filename=file_U)  &$
+;stop
+  vv=read_ncdf('vomecrty',jt,jt, /timestep,iodir=ioDATA,/nostruct,/TOUT,filename=file_V)  &$
 ;stop
 ; lecture salinite & temperature
@@ -77 +204 @@
   conduct(mask_t)=0.
 ; Somme conduct au point T
-
-
+;
 ; Calcul SIGMA
-
+;
   SIGMAoc=total(conduct*e3t3d*tmask,3)
   SIGMA=SIGMAsed+SIGMAoc
-
-
-
+;
   SIGMA_u=(SIGMA+shift(SIGMA,-1,0))/2.
-
   SIGMA_v=(SIGMA+shift(SIGMA,0,1))/2.
-
+;
 ; Calcul B en points u et v
-
+;
   BR_u=(BR+shift(BR,-1,0))/2.
-
   BR_v=(BR+shift(BR,0,1))/2.
-
-
+;
 ; Calcul integrale conduct
-
+;
   conduct_u=(conduct+shift(conduct,-1,0,0))/2.
-
   conduct_v=(conduct+shift(conduct,0,1,0))/2.
-  
+;
   u_cond_u=total( vu*conduct_u*e3u3d*umask(),3)
-
   v_cond_v=total( vv*conduct_v*e3v3d*vmask(),3)
-
-
+;
   Bu_star= BR_u*u_cond_u/SIGMA_u
-
   Bv_star= BR_v*v_cond_v/SIGMA_v
-
+;
 ; Divergence du champ
-
-  Diver=divfred(Bu_star,Bv_star)
+  Diver=div(Bu_star,Bv_star)
 
 Diver=Diver*1e-6
@@ -120 +236 @@
 lecontinent=where(Diver GT 1.E08)
 Diver(lecontinent)=0.
-
-
 
 ;stop
@@ -134 +248 @@
 sigma3(181,*,*,*)=sigma3(1,*,*,*)
 
-
 print,  jt
 
-
-ENDFOR  
+ENDFOR
 ; on ferme le NetCDF
 ;NCDF_CLOSE,id3
 ;NCDF_CLOSE,id4
-
-
-
 
 temps=fltarr(73)
@@ -154 +263 @@
 
    vargrid = 'T'
-   iodir = '/usr/work/sur/fvi/OPA/ORCA2/'
+   iodir = geomag_od_env
 ; Nom
    idout = NCDF_CREATE(iodir+'DivBustar_5d_'+iyear+'_grid_T.nc',/clobber)
@@ -203 +312 @@
 ; Creation de la latitude
    NCDF_VARPUT, idout, id1, gphit
-; Creation de la profondeur 
-   NCDF_VARPUT, idout, id2, gdept 
+; Creation de la profondeur
+   NCDF_VARPUT, idout, id2, gdept
 ; Creation du calendrier
 
    NCDF_VARPUT, idout, id3, temps
 
-   
+
 ; Ecriture des donnees
 
-; ecriture des glam_8 
+; ecriture des glam_8
    NCDF_VARPUT, idout, id4 , diver3
 
@@ -267 +376 @@
 ; Creation de la latitude
    NCDF_VARPUT, idout, id1, gphit
-; Creation de la profondeur 
-   NCDF_VARPUT, idout, id2, gdept 
+; Creation de la profondeur
+   NCDF_VARPUT, idout, id2, gdept
 ; Creation du calendrier
 
    NCDF_VARPUT, idout, id3, temps
 
-   
+
 ; Ecriture des donnees
 
-; ecriture des glam_8 
+; ecriture des glam_8
    NCDF_VARPUT, idout, id4 , sigma3