Changeset 48 for trunk/forcagequimarche.pro

Timestamp:

03/16/14 20:38:39 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules

File:

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

trunk/forcagequimarche.pro

@@ -2 +2 @@
 ;
 ; @file_comments
-; calcul d'une matrice 3D de conductivité sigma=f(T,S) en fonction de T et S
+; calcul d'une matrice 3D de conductivité sigma=f(T,S) en fonction de T et S
 ; de ORCA, sur la grille T
 ;
@@ -67 +67 @@
 ;
 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
+    ;
+    ;----
+    ; 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
+    ;
+    ; 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 = 'eee : invalid orcares parameter = ' + orcares
-                 ras = report(msg)
-                 msg = 'eee : orcares must be ORCA2 or ORCA025++'
-                 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 iyear
-; ++ check ian
-
-;++@init2
-; init grid
-initocemeshmask, orcares, DRAKKAR_EXP = drakkar_exp
-
-@common
-;
-; 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'
-key_portrait = 0
-; 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'
-jpt = 73 ;time counter des fichiers ci-dessus
-     END
-  ENDCASE
-  file_Sed= geomag_id_env+'cond_sed_'+orcares+'.nc'
-  file_Br= geomag_id_env+'Br_'+orcares + '.nc'
-;
-; title
-     t_exp= e_exp
-     t_bt  = 'bar_transp'
-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')
-e3u3d=read_ncdf('e3u_ps',0,/timestep,iodir=ioORLN2,/nostruct,/tout,filename='meshmask_bab.nc')
-SIGMAsed=read_ncdf('cond_sed',0,/timestep,/nostruct,/tout,filename=file_Sed,/cont_nofill)
-;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)
-        for k=0, jpk-1 do begin              &$
-          for j=0,jpj-1 do begin              &$
-            for i=0,jpi-1 do begin             &$
-              e3t3d(i,j,k) = e3t(k)    &$
-            endfor                     &$
-          endfor                      &$
-        endfor
-;
-;vud = make_array(jpi,jpj,jpt)
-;vvd = make_array(jpi, jpj, jpt)
-divBustar = make_array(jpi, jpj, jpt)
-diver3=replicate(0.,182,149,1,jpt)
-sigma3=replicate(0.,182,149,1,jpt)
-;
-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)  &$
+  ;
+  ; 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'
+  key_portrait = 0
+  ; 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'
+  jpt = 73 ;time counter des fichiers ci-dessus
+       END
+    ENDCASE
+    file_Sed= geomag_id_env+'cond_sed_'+orcares+'.nc'
+    file_Br= geomag_id_env+'Br_'+orcares + '.nc'
+  ;
+  ; title
+       t_exp= e_exp
+       t_bt  = 'bar_transp'
+  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')
+  e3u3d=read_ncdf('e3u_ps',0,/timestep,iodir=ioORLN2,/nostruct,/tout,filename='meshmask_bab.nc')
+  SIGMAsed=read_ncdf('cond_sed',0,/timestep,/nostruct,/tout,filename=file_Sed,/cont_nofill)
+  ;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)
+          for k=0, jpk-1 do begin              &$
+            for j=0,jpj-1 do begin              &$
+              for i=0,jpi-1 do begin             &$
+                e3t3d(i,j,k) = e3t(k)    &$
+              endfor                     &$
+            endfor                      &$
+          endfor
+  ;
+  ;vud = make_array(jpi,jpj,jpt)
+  ;vvd = make_array(jpi, jpj, jpt)
+  divBustar = make_array(jpi, jpj, jpt)
+  diver3=replicate(0.,182,149,1,jpt)
+  sigma3=replicate(0.,182,149,1,jpt)
+  ;
+  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)  &$
+  ;stop
+  ; lecture salinite & temperature
+    temp= read_ncdf('votemper',jt,jt,/timestep,iodir=ioDATA,/nostruct,/tout,filename=file_T)
+  ;stop
+    salin=read_ncdf('vosaline',jt,jt,/timestep,iodir=ioDATA,/nostruct,/tout,filename=file_T)
+  ;stop
+    conduct=0.02047780622061 + 0.00273147624197*temp + 0.00035133182334*temp*temp + 0.09139808809909*salin + 0.00241425798890*salin*temp -0.00023998958774*salin*salin
+    mask_t=where(conduct GT 1.e+19)
+    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
+  
+  ; Transport horizontal
+    T_u=total( vu*e3u3d*umask(),3)
+    T_v=total( vv*e3v3d*vmask(),3)
+  
+  ;
+  ; Divergence du champ
+    Diver=divfred(Bu_star,Bv_star)
+  
+  Diver=Diver*1e-6
+  ;stop
+  lecontinent=where(Diver GT 1.E08)
+  Diver(lecontinent)=0.
+
 ;stop
-  vv=read_ncdf('vomecrty',jt,jt, /timestep,iodir=ioDATA,/nostruct,/TOUT,filename=file_V)  &$
-;stop
-; lecture salinite & temperature
-  temp= read_ncdf('votemper',jt,jt,/timestep,iodir=ioDATA,/nostruct,/tout,filename=file_T)
-;stop
-  salin=read_ncdf('vosaline',jt,jt,/timestep,iodir=ioDATA,/nostruct,/tout,filename=file_T)
-;stop
-  conduct=0.02047780622061 + 0.00273147624197*temp + 0.00035133182334*temp*temp + 0.09139808809909*salin + 0.00241425798890*salin*temp -0.00023998958774*salin*salin
-  mask_t=where(conduct GT 1.e+19)
-  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
-
-; Transport horizontal
-  T_u=total( vu*e3u3d*umask(),3)
-  T_v=total( vv*e3v3d*vmask(),3)
-
-;
-; Divergence du champ
-  Diver=divfred(Bu_star,Bv_star)
-
-Diver=Diver*1e-6
-;stop
-lecontinent=where(Diver GT 1.E08)
-Diver(lecontinent)=0.
-
-;stop
-;bande de recouvrement::
+;bande de recouvrement:
 
 diver3(1:180,0:147,*,jt)=Diver(*,*)
@@ -342 +342 @@
 
    NCDF_CLOSE, idout
-; ++ l'équivalent avec forcage_output
+; ++ l'équivalent avec forcage_output
 forcage_output, orcares, variable, 'Divergence', long_name, jpiglo, jpjglo, gphit, glamt,diver3