Changeset 48 for trunk/forcage.pro

Timestamp:

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

Author:

pinsard

Message:

fix thanks to coding rules

File:

• trunk/forcage.pro (modified) (11 diffs)

trunk/forcage.pro

@@ -1 @@
-PRO forcagequimarche,iyear,ian 
+PRO forcagequimarche,iyear,ian
 @init2
 @initorca2_bab
 
-
-
-
 @common
-
-
 
 ;ian='01'
@@ -16 +11 @@
 rep_fred='/usr/work/sur/fvi/OPA/geomag/'
 key_portrait = 0
-; stockage des fichiers brut 
+; stockage des fichiers brut
   ioDATA='/usr/work/sur/fvi/OPA/ORCA2/'
-  file_U=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_U.nc' 
+  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' 
+  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' 
+  file_T=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_T.nc'
   print, file_T
   file_Sed= rep_fred+'cond_sed_ORCA2.nc'
@@ -33 +28 @@
      t_bt  = 'bar_transp'
 ioORLN2 = '/usr/work/sur/fvi/OPA/ORCA2'
-;facteur d'echelle vertical  for partial steps
+;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')
@@ -44 +39 @@
 
 ; 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)    &$
@@ -54 +49 @@
 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)
 diver2=fltarr(jpi,jpj,1,73)
 
-; ouverture des fichiers dans lesquels on va écrire 
+; ouverture des fichiers dans lesquels on va écrire
 ;id3=NCDF_OPEN('/usr/work/sur/fvi/OPA/geomag/U_5d_'+iyear+'_grid_T.nc',/write)
 ;id4=NCDF_OPEN('/usr/work/sur/fvi/OPA/geomag/V_5d_'+iyear+'_grid_T.nc',/write)
@@ -66 +61 @@
 ;id4=NCDF_OPEN('/usr/work/sur/fvi/OPA/ORCA2/DivBustar_5d_'+iyear+'_grid_T.nc',/write)
 
-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
@@ -107 +102 @@
 
   conduct_v=(conduct+shift(conduct,0,1,0))/2.
-  
+
   u_cond_u=total( vu*conduct_u*e3u3d*umask(),3)
 
@@ -126 +121 @@
 ;for jj=0,147 do begin &$
 ;  for ji=0,179 do begin &$
-;   if( Diver(ji,jj) GT 1e10 ) then begin &$ 
-;       Diver(ji,jj) = 0. &$
+;   if( Diver(ji,jj) GT 1e10 ) then begin &$
+;      Diver(ji,jj) = 0. &$
 ;   endif &$
 ;  endfor &$
-;endfor 
+;endfor
 
 ;stop
@@ -137 +132 @@
 ;stop
 ; Somme sur la verticale partial steps
-;                vum=total( vu*e3u3d*umask(),3 )  &$   
-;                vvm=total( vv*e3v3d*vmask(),3 )  &$   
+;                vum=total( vu*e3u3d*umask(),3 )  &$
+;                vvm=total( vv*e3v3d*vmask(),3 )  &$
 
 ; Shift sur la grille T partial steps
-;             vut= (vum+shift(vum,1,0) )*0.5   &$   
+;             vut= (vum+shift(vum,1,0) )*0.5   &$
 ;             vvt= (vvm+shift(vvm,0,1) )*0.5   &$
-; Bande de recouvrement 
+; Bande de recouvrement
 ;             vut(0, *) = vut(jpi-2, *)
 ;             vvt(*, 0) = 0.
-; stockage dans le fichier de sortie 
+; stockage dans le fichier de sortie
 ;NCDF_VARPUT, id3,'sossheig',vut, offset = [0, 0, jt]
 ;NCDF_VARPUT, id4,'sossheig',Diver, offset = [0, 0, jt]
@@ -153 +148 @@
 
 
-ENDFOR  
+ENDFOR
 ; on ferme le NetCDF
 ;NCDF_CLOSE,id3
@@ -218 +213 @@
 ; 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 , diver2