Changeset 48 for trunk/FIG_bar_transp_write.pro
- Timestamp:
- 03/16/14 20:38:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FIG_bar_transp_write.pro
r2 r48 1 PRO FIG_bar_transp_write,iyear,ian 1 PRO FIG_bar_transp_write,iyear,ian 2 2 @init2 3 3 @initorca2_bab … … 15 15 ;e_year0='10' 16 16 key_portrait = 0 17 ; stockage des fichiers brut 17 ; stockage des fichiers brut 18 18 ioDATA='/usr/work/sur/fvi/OPA/ORCA2/'+e_exp 19 file_U=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_U.nc' 19 file_U=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_U.nc' 20 20 print, file_U 21 file_V=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_V.nc' 21 file_V=e_exp+'_5d_'+ian+'0101_'+ian+'1231_grid_V.nc' 22 22 print, file_V 23 23 … … 27 27 t_bt = 'bar_transp' 28 28 ioORLN2 = '/usr/work/sur/fvi/OPA/ORCA2' 29 ;facteur d'echelle vertical 29 ;facteur d'echelle vertical for partial steps 30 30 e3v3d=read_ncdf('e3v_ps',0,/timestep,iodir=ioORLN2,/nostruct,/tout,filename='meshmask_bab.nc') 31 31 e3u3d=read_ncdf('e3u_ps',0,/timestep,iodir=ioORLN2,/nostruct,/tout,filename='meshmask_bab.nc') … … 34 34 35 35 ; vertical integration: 36 e3t3d=make_array(jpi,jpj,jpk) 36 e3t3d=make_array(jpi,jpj,jpk) 37 37 for k=0, jpk-1 do begin &$ 38 for j=0,jpj-1 do begin &$ 38 for j=0,jpj-1 do begin &$ 39 39 for i=0,jpi-1 do begin &$ 40 40 e3t3d(i,j,k) = e3t(k) &$ … … 44 44 jpt = 73 45 45 46 ;vud = make_array(jpi,jpj,jpt) 46 ;vud = make_array(jpi,jpj,jpt) 47 47 ;vvd = make_array(jpi, jpj, jpt) 48 48 49 ; ouverture des fichiers dans lesquels on va écrire49 ; ouverture des fichiers dans lesquels on va écrire 50 50 id3=NCDF_OPEN('/usr/work/sur/fvi/geomag/U_5d_'+iyear+'_grid_T.nc',/write) 51 51 id4=NCDF_OPEN('/usr/work/sur/fvi/geomag/V_5d_'+iyear+'_grid_T.nc',/write) 52 52 53 FOR jt = 0, jpt-1 DO BEGIN &$ 53 FOR jt = 0, jpt-1 DO BEGIN &$ 54 54 55 55 ; ouverture des fichiers et stockage en memoire partial steps 56 vu=read_ncdf('vozocrtx',jt,jt, /timestep, iodir=ioDATA,/nostruct,/TOUT,filename=file_U) &$ 57 vv=read_ncdf('vomecrty',jt,jt, /timestep,iodir=ioDATA,/nostruct,/TOUT,filename=file_V) &$ 56 vu=read_ncdf('vozocrtx',jt,jt, /timestep, iodir=ioDATA,/nostruct,/TOUT,filename=file_U) &$ 57 vv=read_ncdf('vomecrty',jt,jt, /timestep,iodir=ioDATA,/nostruct,/TOUT,filename=file_V) &$ 58 58 59 59 60 60 61 61 ; Somme sur la verticale partial steps 62 vum=total( vu*e3u3d*umask(),3 ) &$ 63 vvm=total( vv*e3v3d*vmask(),3 ) &$ 62 vum=total( vu*e3u3d*umask(),3 ) &$ 63 vvm=total( vv*e3v3d*vmask(),3 ) &$ 64 64 65 65 ; Shit sur la grille T partial steps 66 vut= (vum+shift(vum,1,0) )*0.5 &$ 66 vut= (vum+shift(vum,1,0) )*0.5 &$ 67 67 vvt= (vvm+shift(vvm,0,1) )*0.5 &$ 68 ; Bande de recouvrement 68 ; Bande de recouvrement 69 69 vut(0, *) = vut(jpi-2, *) 70 70 vvt(*, 0) = 0. 71 ; stockage dans le fichier de sortie 71 ; stockage dans le fichier de sortie 72 72 NCDF_VARPUT, id3,'sossheig',vut, offset = [0, 0, jt] 73 73 NCDF_VARPUT, id4,'sossheig',vvt, offset = [0, 0, jt] … … 75 75 print, jt 76 76 77 ENDFOR 77 ENDFOR 78 78 ; on ferme le NetCDF 79 79 NCDF_CLOSE,id3 … … 107 107 ; tidout = NCDF_DIMDEF(idout, 'time_counter', /unlimited) 108 108 109 109 110 110 ; Attributs globaux 111 111 ; id0 = NCDF_VARDEF(idout, 'nav_lon' , [xidout, yidout ], /FLOAT) … … 141 141 ; NCDF_ATTPUT, idout, id5, 'title', 'vomecrty' 142 142 ; NCDF_ATTPUT, idout, id5, 'long_name', 'meridian current summmed on vertical' 143 143 144 144 ; NCDF_CONTROL, idout, /ENDEF 145 145 … … 148 148 ; Creation de la latitude 149 149 ; NCDF_VARPUT, idout, id1, gphit 150 ; Creation de la profondeur 150 ; Creation de la profondeur 151 151 ; prof=1. 152 ; NCDF_VARPUT, idout, id2, prof 152 ; NCDF_VARPUT, idout, id2, prof 153 153 ; Creation du calendrier 154 154 ; temps = findgen(jpt)*5.*86400. +julday(1, 1, 1) … … 158 158 ; NCDF_VARPUT, idout, id3, temps 159 159 160 160 161 161 ; Ecriture des donnees 162 162
Note: See TracChangeset
for help on using the changeset viewer.