Changeset 4320 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM
- Timestamp:
- 2013-11-28T12:41:27+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 23 added
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_main.sh
r2751 r4320 173 173 pdflist=${pdflist}' '${PDF_DIR}/${file}pdf 174 174 done 175 #pdfjam needed to use pdfjoin 176 pdfjoin ${pdflist} --outfile $PDF_DIR/${filepdf} 175 177 texexec --pdfarrange --result=$PDF_DIR/$filepdf $pdflist 176 178 rm -f $PDF_DIR/$( basename $filepdf .pdf ).aux $PDF_DIR/$( basename $filepdf .pdf ).log 177 echo ${PDF_DIR}/$filepdf done 179 if [ ! -e ${PDF_DIR}/$filepdf ] 180 then 181 echo "\n ATTENTION !!!!! \n" 182 echo "${PDF_DIR}/$filepdf do not exist \n" 183 else 184 echo "${PDF_DIR}/$filepdf done" 185 fi 178 186 #commented because convert gives images of bad quality 179 187 # convert ${PS_DIR}/*.ps /tmp/all_${PLOTTYPE}.pdf -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.pro
r2751 r4320 26 26 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 27 27 ; 28 domdef, 20, 380, 60, 9028 domdef, 20, 380, 50, 90 29 29 ; 30 30 varunit = S1.unit -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_GlobMeanTS.pro
r2751 r4320 25 25 title = 'Salinity Global mean!C'+std_file1_T+' - Levitus (Black)' 26 26 if std_file1_T NE std_file2_T THEN title = title+'!C'+std_file2_T+' - Levitus (Red)' 27 plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = . 2, boxzoom = [4., 5300.], /KEEPBOTTOM $27 plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = .4, boxzoom = [4., 5300.], /KEEPBOTTOM $ 28 28 , small = [1, 2, 2], XGRIDSTYLE = 2, YGRIDSTYLE = 2, TITLE = title, /NOERASE, _extra = ex 29 29 if std_file1_T NE std_file2_T then begin 30 plt1d, S2.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -2., MAX = 2., boxzoom = [4., 5300.], /KEEPBOTTOM $30 plt1d, S2.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -2., MAX = 4., boxzoom = [4., 5300.], /KEEPBOTTOM $ 31 31 , /ov1d, COLOR = 250, TITLE = title, _extra = ex 32 32 endif -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r2751 r4320 40 40 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 41 41 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 42 endif 43 42 ;;SF finire, aggiungere lettura obs 43 ;;SF endif else begin 44 ;;SF title = title+ ' - Observations' 45 ;;SF plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $ 46 ;;SF , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $ 47 ;;SF , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 48 ;;SF endelse 49 endif 50 44 51 domdef 45 52 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro
r2751 r4320 1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 12 12 IF keyword_set(arc) THEN var = var+'_Arc_' 13 13 IF keyword_set(ant) THEN var = var+'_Ant_' 14 IF keyword_set( march) THEN var = var+'March'14 IF keyword_set(april) THEN var = var+'April' 15 15 IF keyword_set(sept) THEN var = var+'Sept' 16 16 17 17 filename = cdti3 + '_'+var+'_'+std_file1_I 18 ifstd_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I19 ifKEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 118 IF std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I 19 IF KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 20 20 21 21 IF keyword_set(arc) THEN BEGIN 22 22 domdef, 20, 380, 50, 90 23 23 vmin = 0. 24 vmax = 10.24 vmax = 6. 25 25 vint = 0.5 26 fmt = '(i2)'26 fmtd = '(f4.1)' 27 27 vmind = -3. 28 28 vmaxd = 3. 29 29 vintd = 0.5 30 fmtd = '(i2)' 30 fmtd = '(f4.1)' 31 div = 6 31 32 map = [90, 0, 0] 32 33 ENDIF … … 41 42 vintd = 0.2 42 43 fmtd = '(f4.1)' 44 div = 5 43 45 map = [-90, 0, 0] 44 46 ENDIF … … 48 50 title = var+'!C'+std_file1_I 49 51 plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 50 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $52 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div $ 51 53 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 52 54 ; ; … … 54 56 title = title + std_file2_I 55 57 plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 56 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $58 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 57 59 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 58 60 endif -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r2755 r4320 15 15 PRINT, '' 16 16 ; 17 std_iodir_data = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format')18 std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data')19 std_iodir_mask = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)')17 std_iodir_data = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format') 18 std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data') 19 std_iodir_mask = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)') 20 20 ; meshmask 21 std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask)21 std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask) 22 22 std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask) 23 23 24 ; Levitus 98 25 std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 26 std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 27 std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 28 std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 29 std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 24 ; climatologies 25 std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 26 std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 27 std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 28 std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 29 std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 30 std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato) 31 std_file_snow_arc = isafile(getenv('FILE_SNOW_ARC'), title = 'SNOW_ARC', iodir = std_iodir_climato) 32 std_file_snow_ant = isafile(getenv('FILE_SNOW_ANT'), title = 'SNOW_ANT', iodir = std_iodir_climato) 30 33 31 34 IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $ … … 155 158 Ithi_1 = read_ncdf(getenv('VAR1_Ithick'), allrecords = allrec, filename = std_file1_I ) 156 159 caldat, time, mm 157 march = where(mm EQ 3, cnt)158 Ithi_ march_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(march)], 3), unit:Ithi_1.unit}160 april = where(mm EQ 4, cnt) 161 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 159 162 sept = where(mm EQ 9, cnt) 160 Ithi_sept_1 = {arr:1./float(cnt) * total( Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit}163 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 161 164 undefine, Ithi_1 162 165 ; 163 166 IF std_file2_I NE std_file1_I THEN BEGIN 164 Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )165 caldat, time, mm166 march = where(mm EQ 3, cnt)167 Ithi_march_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(march)], 3), unit:Ithi_2.unit}168 sept = where(mm EQ 9, cnt)169 Ithi_sept_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(sept)], 3), unit:Ithi_2.unit}170 undefine, Ithi_2167 Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I ) 168 caldat, time, mm 169 april = where(mm EQ 4, cnt) 170 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 171 sept = where(mm EQ 9, cnt) 172 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 173 undefine, Ithi_2 171 174 ENDIF ELSE BEGIN 172 Ithi_march_2 = {arr:-1}173 Ithi_sept_2 = {arr:-1}175 Ithi_april_2 = {arr:-1} 176 Ithi_sept_2 = {arr:-1} 174 177 ENDELSE 175 178 ; 176 179 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 180 help, Ifra_1.arr, jpi,jpj 177 181 caldat, time, mm 178 182 march = where(mm EQ 3, cnt) 179 Ifra_march_1 = {arr:1./float(cnt) * total( Ifra_1.arr[*, *, temporary(march)], 3), unit:Ifra_1.unit}183 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 180 184 sept = where(mm EQ 9, cnt) 181 Ifra_sept_1 = {arr:1./float(cnt) * total( Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit}185 Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 182 186 undefine, Ifra_1 183 187 ; … … 186 190 caldat, time, mm 187 191 march = where(mm EQ 3, cnt) 188 Ifra_march_2 = {arr:1./float(cnt) * total( Ifra_2.arr[*, *, temporary(march)], 3), unit:Ifra_2.unit}192 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 189 193 sept = where(mm EQ 9, cnt) 190 Ifra_sept_2 = {arr:1./float(cnt) * total( Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit}194 Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 191 195 undefine, Ifra_2 192 196 ENDIF ELSE BEGIN … … 194 198 Ifra_sept_2 = {arr:-1} 195 199 ENDELSE 200 ; 201 ;; Sf commentato perché non funziona snow 202 ;; Sf Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 203 ;; Sf help, Isnow_1.arr, jpi,jpj 204 ;; Sf caldat, time, mm 205 ;; Sf april = where(mm EQ 4, cnt) 206 ;; Sf Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 207 ;; Sf sept = where(mm EQ 9, cnt) 208 ;; Sf Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 209 ;; Sf ;undefine, Isnow_1 210 ;; Sf ; 211 ;; Sf IF std_file2_I NE std_file1_I THEN BEGIN 212 ;; Sf Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 213 ;; Sf caldat, time, mm 214 ;; Sf april = where(mm EQ 4, cnt) 215 ;; Sf Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 216 ;; Sf sept = where(mm EQ 9, cnt) 217 ;; Sf Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 218 ;; Sf ; undefine, Isnow_2 219 ;; Sf ENDIF ELSE BEGIN 220 ;; Sf Isnow_april_2 = {arr:-1} 221 ;; Sf Isnow_sept_2 = {arr:-1} 222 ;; Sf ENDELSE 196 223 ; 197 224 jpt = 1 … … 299 326 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 300 327 301 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 302 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 303 328 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 329 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 330 331 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 332 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 333 304 334 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 305 335 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 306 336 307 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 308 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 337 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 338 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 339 340 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 341 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 309 342 310 343 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' … … 322 355 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 323 356 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 357 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 358 ;;SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 359 ;;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 360 ;; SF 361 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 362 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 363 ;; SF 364 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 365 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 366 ;; SF 367 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 368 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 369 ; 370 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 371 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 372 373 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 374 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 324 375 325 376 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)' -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_bsf.pro
r2751 r4320 14 14 domdef, 0, 6000 15 15 16 bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0] )16 bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0], /transport ) 17 17 title = 'Barotropic Stream Function!C'+std_file1_U 18 18 plt, bb1, min = -200., max = 200., int = 10., /portrait, FORMAT = '(I4)', STYLE = 'so0so' $ … … 20 20 21 21 if std_file1_U NE std_file2_U then BEGIN 22 bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0] )22 bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0], /transport ) 23 23 title = 'Barotropic Stream Function!C'+std_file1_U+' - '+std_file2_U 24 24 plt, bb1.arr - bb2.arr, min = -20., max = 20., int = 2., FORMAT = '(I2)', STYLE = 'so0so' $ -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mld.pro
r2751 r4320 32 32 ENDIF ELSE BEGIN 33 33 title = titleorg+std_file1_T 34 plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 35 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 34 plt, MLD1,label=4, cb_label=[0, 15, 25, 50, 75, 100, 125, 150, 250, 400, 600, 850],/NOCONTOUR,FORMAT = '(I3)' $ 35 , small = [1, 2, 1],COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 36 ;ORI SF plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 37 ;ORI SF , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 36 38 ENDELSE 37 39 ; -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro
r2760 r4320 33 33 domdef, 0, 6000 34 34 35 mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma )35 mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma, /transport ) 36 36 title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V 37 37 pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $ … … 39 39 40 40 if std_file1_V NE std_file2_V then begin 41 mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma )41 mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma, /transport) 42 42 title = title+' - '+std_file2_V 43 43 pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)' $ -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1
r2753 r4320 16 16 #===================== User PATHS ===================== 17 17 # 18 idl_command=/Applications/itt/idl64/bin/idl 18 #idl_command=/Applications/itt/idl64/bin/idl 19 idl_command=/Applications/itt/idl71/bin/idl 20 #idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl 21 #idl_command=idl71 19 22 # 20 PS_DIR=$( pwd )/ps 21 PDF_DIR=$( pwd )/pdf 22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/s masson/SAXO_DIR23 PS_DIR=$( pwd )/ps_plot 24 PDF_DIR=$( pwd )/pdf_plot 25 HTML_DIR=$( pwd )/html_plot 26 SAXO_DIR=/Users/sflod/SAXO_DIR 24 27 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 31 # 29 32 #===================== Model GRID ===================== 30 33 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask34 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 35 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 36 # 34 37 #===================== DATA ===================== 35 38 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 39 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 40 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 41 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 42 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 43 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 44 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 45 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 46 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 47 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 48 VAR_SNOW_NH=SnowDepth ; FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc # Snow arctic 49 VAR_SNOW_SH=sno_dept_mm ; FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc # Snow antartic 41 50 # 42 51 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 60 69 VAR1_Ithick=sit 61 70 VAR1_Ifrac=sic 71 VAR1_SNOW=isnowthi 62 72 # 63 73 #===================== EXP2 ===================== … … 77 87 VAR2_Ithick=sit 78 88 VAR2_Ifrac=sic 89 VAR2_SNOW=isnowthi 79 90 # 80 91 ######################### Export Variables ############################### … … 86 97 export FILE_MESH_MASK FILE_MASK_SUBDOMAIN 87 98 #===================== DATA ===================== 88 export FILE_TEMP_3D VAR_TEMP_3D 89 export FILE_SAL_3D VAR_SAL_3D 90 export FILE_SST VAR_SST 91 export FILE_GEOHEAT VAR_GEOHEAT 92 export FILE_FLUX VAR_FLUX 93 export FILE_MLD VAR_MLD 99 export FILE_TEMP_3D VAR_TEMP_3D 100 export FILE_SAL_3D VAR_SAL_3D 101 export FILE_SST VAR_SST 102 export FILE_GEOHEAT VAR_GEOHEAT 103 export FILE_FLUX VAR_FLUX 104 export FILE_MLD VAR_MLD 105 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 106 export FILE_SNOW_ARC VAR_SNOW_NH 107 export FILE_SNOW_ANT VAR_SNOW_SH 94 108 #===================== EXP1 ===================== 95 109 export READ_ONLY_FIRST_RECORD 96 export FILE1_T FILE1_U FILE1_VFILE1_I97 export VAR1_T VAR1_SVAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD98 export VAR1_U VAR1_V99 export VAR1_Ithick VAR1_Ifrac 110 export FILE1_T FILE1_U FILE1_V FILE1_I 111 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 112 export VAR1_U VAR1_V 113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 100 114 #===================== EXP2 ===================== 101 export FILE2_T FILE2_U FILE2_VFILE2_I102 export VAR2_T VAR2_SVAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD103 export VAR2_U VAR2_V104 export VAR2_Ithick VAR2_Ifrac 115 export FILE2_T FILE2_U FILE2_V FILE2_I 116 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 117 export VAR2_U VAR2_V 118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 105 119 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2
r2753 r4320 16 16 #===================== User PATHS ===================== 17 17 # 18 idl_command=/Applications/itt/idl64/bin/idl 18 #idl_command=/Applications/itt/idl64/bin/idl 19 idl_command=/Applications/itt/idl71/bin/idl 20 #idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl 21 #idl_command=idl71 19 22 # 20 PS_DIR=$( pwd )/ ps21 PDF_DIR=$( pwd )/ pdf22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/s masson/SAXO_DIR23 PS_DIR=$( pwd )/500yfwb0_ps_plot_500y 24 PDF_DIR=$( pwd )/500yfwb0_pdf_plot_500y 25 HTML_DIR=$( pwd )/html_plot 26 SAXO_DIR=/Users/sflod/SAXO_DIR 24 27 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 31 # 29 32 #===================== Model GRID ===================== 30 33 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask34 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 35 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 36 # 34 37 #===================== DATA ===================== 35 38 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 39 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 40 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 41 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 42 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 43 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 44 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 45 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 46 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 47 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 48 VAR_SNOW_NH=SnowDepth ; FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc # Snow arctic 49 VAR_SNOW_SH=sno_dept_mm ; FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc # Snow antartic 41 50 # 42 51 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 60 69 VAR1_Ithick=sit 61 70 VAR1_Ifrac=sic 71 VAR1_SNOW=isnowthi 62 72 # 63 73 #===================== EXP2 ===================== … … 77 87 VAR2_Ithick=sit 78 88 VAR2_Ifrac=sic 89 VAR2_SNOW=isnowthi 79 90 # 80 91 ######################### Export Variables ############################### … … 86 97 export FILE_MESH_MASK FILE_MASK_SUBDOMAIN 87 98 #===================== DATA ===================== 88 export FILE_TEMP_3D VAR_TEMP_3D 89 export FILE_SAL_3D VAR_SAL_3D 90 export FILE_SST VAR_SST 91 export FILE_GEOHEAT VAR_GEOHEAT 92 export FILE_FLUX VAR_FLUX 93 export FILE_MLD VAR_MLD 99 export FILE_TEMP_3D VAR_TEMP_3D 100 export FILE_SAL_3D VAR_SAL_3D 101 export FILE_SST VAR_SST 102 export FILE_GEOHEAT VAR_GEOHEAT 103 export FILE_FLUX VAR_FLUX 104 export FILE_MLD VAR_MLD 105 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH, VAR_ICE_area_NH, VAR_ICE_area_SH 106 export FILE_SNOW_ARC VAR_SNOW_NH 107 export FILE_SNOW_ANT VAR_SNOW_SH 94 108 #===================== EXP1 ===================== 95 109 export READ_ONLY_FIRST_RECORD 96 export FILE1_T FILE1_U FILE1_VFILE1_I97 export VAR1_T VAR1_SVAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD98 export VAR1_U VAR1_V99 export VAR1_Ithick VAR1_Ifrac 110 export FILE1_T FILE1_U FILE1_V FILE1_I 111 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 112 export VAR1_U VAR1_V 113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 100 114 #===================== EXP2 ===================== 101 export FILE2_T FILE2_U FILE2_VFILE2_I102 export VAR2_T VAR2_SVAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD103 export VAR2_U VAR2_V104 export VAR2_Ithick VAR2_Ifrac 115 export FILE2_T FILE2_U FILE2_V FILE2_I 116 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 117 export VAR2_U VAR2_V 118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 105 119 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_zonal_S.pro
r2751 r4320 33 33 ; 34 34 title = 'Salinity, '+subname+'!C'+std_file1_T 35 pltz, S1, MININ = 3 2., MAXIN = 37., INTER = .2, typein = 'yz', FORMAT = '(f4.1)' $35 pltz, S1, MININ = 31., MAXIN = 36.8, INTER = .2, typein = 'yz', FORMAT = '(f4.1)' $ 36 36 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT, _extra = ex 37 37 ; -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro
r2751 r4320 14 14 ts_Time = 0. 15 15 16 trans = strlowcase(var_name) EQ 'vocetr_eff' 16 ;SF trans = strlowcase(var_name) EQ 'vocetr_eff' 17 trans = strlowcase(var_name) EQ 'vomeeftr' 17 18 19 ; ADD definition of mask_atl_nomed to compute AMOC (with no Mediterranean Sea) 20 msk_atl_nomed = read_ncdf( 'atlmsk_nomed', filename = std_file_msksub, /nostruct, _extra = ex ) 18 21 FOR i = 0, nfiles-1 DO BEGIN 19 22 var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct) … … 21 24 ; 22 25 FOR t = 0, jpt-1 DO BEGIN 26 msfatl = msf(var[*, *, *, t], msk_atl_nomed, TRANSPORT = trans, /nostruct, indexboxzoom = ind) 23 27 msftot = msf(var[*, *, *, t], TRANSPORT = trans, /nostruct, indexboxzoom = ind) 24 28 yaxis = gphit[ind[0], ind[2]:ind[3]] … … 27 31 indy = where(yaxis gt 40 and yaxis le 50) 28 32 domdef, 0, 3500 29 t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 33 ;SF commented because we've to compute msftot starting from msf atlantic; i.e. without mediterranean sea. 34 ;SF t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 35 t45 = [t45, max(msfatl[indy, firstzw:lastzw], /NaN)] 30 36 31 37 ; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S … … 78 84 79 85 ; compute the MSF 80 86 ; 81 87 std_ts_read_AMOC, vV1, date1, date2, prefix, suffix, t45, t70, t50 82 88 if prefix NE prefix2 then begin … … 114 120 if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 115 121 116 title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2122 title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 117 123 pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $ 118 124 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, _extra = ex -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Drake.pro
r2751 r4320 32 32 u1 = rseries_ncdf(vU1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) 33 33 drk1 = fltarr(jpt) 34 trans = strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' 34 ;SF trans = strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' 35 trans = strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' 35 36 FOR t = 0, jpt-1 DO BEGIN 36 37 tmp = bsf(u1[*, *, *, t], TRANSPORT = trans, /nostruct) … … 48 49 u2 = rseries_ncdf(vU2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, /nostruct) 49 50 drk2 = fltarr(jpt) 50 trans = strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' 51 ;SF trans = strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' 52 trans = strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' 51 53 FOR t = 0, jpt-1 DO BEGIN 52 54 tmp = bsf(u2[*, *, *, t], TRANSPORT = trans, /nostruct) -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro
r2751 r4320 27 27 ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce ) & ts_EMP.unit = 'Sv' 28 28 title = prefix+' '+d1_d2+'!C'+blabla 29 pltt, ts_EMP, 't', - .001, .001, date1, date2, /REMPLI $29 pltt, ts_EMP, 't', -2., 2., date1, date2, /REMPLI $ 30 30 , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 31 31 … … 39 39 40 40 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', - .001, .001, date1, date2, /REMPLI $41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', -2., 2., date1, date2, /REMPLI $ 42 42 , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 43 43 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro
r2751 r4320 1 pro std_ts_ICE, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 1 function read_arr2d, filename, varname, t1, t2 2 ;; function that read input file and return 2d array with monthly timecounter 3 nyear = (t2-t1+1)/12 4 arr2d = ncdf_lec(filename, VAR=varname) 5 arr2d = arr2d[t1:t2] 6 arr2d = reform(arr2d,12,nyear) ; put in 2D array 7 arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 8 9 return, arr2d 10 end 11 12 ;; here start procedure that use function read_arr2d 13 pro std_ts_ICE, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 2 14 3 15 compile_opt idl2, strictarrsubs … … 10 22 ; get exp2 info 11 23 vICE2 = getenv('VAR2_ICE') & prefix2 = getenv('V2ICE_PREF') & suffix2 = getenv('V2ICE_SUFF') 24 ; get ice climatology info 25 std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE Extent Climatology', iodir = std_iodir_climato) 12 26 ; 27 time_ice = ncdf_lec( std_file_ice, VAR='time' ) 28 time_ice = (time_ice - FLOOR(time_ice) ) * 12 29 time_ice = (round(time_ice) + 11) mod 12; round to nearest integer 30 t1 = where(time_ice eq 0) 31 t1 = t1[0] ; jannuary 32 t2 = where(time_ice eq 11, count) 33 t2 = t2[count-1] ; last day of december 34 35 vICE_ext_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_NH'), t1, t2 ) 36 vICE_ext_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_SH'), t1, t2 ) 37 ; 38 vICE_area_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_NH'), t1, t2 ) 39 vICE_area_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_SH'), t1, t2 ) 40 ; 13 41 cdti3 = string(cnt, format = '(i3.3)') 14 42 print, cdti3 + ') ' + blabla 15 filename = cdti3 + '_ts_ AMOC_'+prefix43 filename = cdti3 + '_ts_ICE_'+prefix 16 44 if prefix NE prefix2 then filename = filename + '_'+prefix2 17 if KEYWORD_SET(postscript) then openps, filename+' _1.ps', portrait = 145 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 18 46 ; 19 47 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 20 48 ; 21 49 iodir = std_iodir_data 50 ; ICE Area in NORTH Hemisphere 51 domdef, 0, jpi-1, 30, 90, /xindex 52 ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 53 ICE_N = grossemoyenne(ICE_N.arr, 'xy', /integration, mask2d = masknp) 22 54 23 domdef, 0, jpi-1, 30, 90, /xindex 24 ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 25 ICE_N.arr = ICE_N.arr * 1.e-12 & ICE_N.unit = '10^12 m^2' 26 domdef, 0, jpi-1, -90, -30, /xindex 27 ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 28 ICE_S.arr = ICE_S.arr * 1.e-12 & ICE_S.unit = '10^12 m^2' 55 if jpt mod 12 ne 0 then stop 56 nyr=jpt/12. 57 ICE_N = reform(ice_n, 12,nyr) 58 ICE_n = total(ice_n,2)/nyr 59 ICE_N = {arr:ICE_N * 1.e-12, unit : '10^12 m^2'} 29 60 30 title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 31 pltt, ICE_N, 't', 0., 15., date1, date2, /REMPLI, /PORTRAIT $ 61 ; ICE EXTENT (Area minus 15%) in NORTH Hemisphere 62 ICE_N_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 63 msk = ICE_N_15.arr gt 0.15 ; remove 0.15% for observations 64 ICE_N_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 65 if jpt mod 12 ne 0 then stop 66 nyr=jpt/12. 67 ICE_N_15 = reform(ice_n_15, 12,nyr) 68 ICE_n_15 = total(ice_n_15,2)/nyr 69 ICE_N_15 = {arr:ICE_N_15 * 1.e-12, unit : '10^12 m^2'} 70 ; 71 ;ICE Area in SOUTH Hemisphere 72 domdef, 0, jpi-1, -90, -30, /xindex 73 ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 74 ICE_S = grossemoyenne(ICE_S.arr, 'xy', /integration, mask2d = masknp) 75 if jpt mod 12 ne 0 then stop 76 nyr=jpt/12. 77 ICE_S = reform(ice_S, 12,nyr) 78 ICE_S = total(ice_S,2)/nyr 79 ICE_S = {arr:ICE_S * 1.e-12, unit : '10^12 m^2'} 80 ; ICE EXTENT (Area minus 15%) in SOUTH Hemisphere 81 ICE_S_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 82 msk = ICE_S_15.arr gt 0.15 ; remove 0.15% for observations 83 ICE_S_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 84 if jpt mod 12 ne 0 then stop 85 nyr=jpt/12. 86 ICE_S_15 = reform(ice_S_15, 12,nyr) 87 ICE_S_15 = total(ice_S_15,2)/nyr 88 ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 89 ; 90 title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 91 jpt=12 92 time=julday(1,15,1900)+30*lindgen(12) 93 pltt, ICE_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 4., MAX = 16. $ 32 94 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 33 34 title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 35 pltt, ICE_S, 't', 0., 15., date1, date2, /REMPLI, /NOERASE $ 95 pltt, ICE_N_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; dashed lines is LINESTYLE=2 $ 96 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 97 pltt, vICE_area_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 98 , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 99 pltt, vICE_ext_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 100 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 101 ; 102 title ='Southern Hemisphere !C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 103 pltt, ICE_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 20. $ 36 104 , small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 37 38 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png /> ' ] 105 pltt, ICE_S_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ 106 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 107 pltt, vICE_area_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 108 , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 109 pltt, vICE_ext_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 110 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 111 ; 112 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 39 113 if KEYWORD_SET(postscript) then closeps 40 114 … … 42 116 43 117 d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 44 45 118 tsave = time 46 119 domdef, 0, jpi-1, 30, 90, /xindex 47 ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 48 ICE_N2.arr = ICE_N2.arr * 1.e-12 & ICE_N2.unit = '10^12 m^2' 49 domdef, 0, jpi-1, -90, -30, /xindex 50 ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 51 ICE_S2.arr = ICE_S2.arr * 1.e-12 & ICE_S2.unit = '10^12 m^2' 52 time = tsave & IF n_elements(time) NE jpt THEN stop 120 ;ICE Extent in NORTH Hemisphere 121 ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 122 ICE_N2 = grossemoyenne(ICE_N2.arr, 'xy', /integration, mask2d = masknp) 123 if jpt mod 12 ne 0 then stop 124 nyr=jpt/12. 125 ICE_N2 = reform(ICE_N2, 12,nyr) 126 ICE_N2 = total(ICE_N2,2)/nyr 127 ICE_N2 = {arr:ICE_N2 * 1.e-12, unit : '10^12 m^2'} 128 ;ICE Extent minus 15% in NORTH Hemisphere 129 ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 130 msk = ICE_N2_15.arr gt 0.15 ; remove 0.15% for observations 131 ICE_N2_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 132 if jpt mod 12 ne 0 then stop 133 nyr=jpt/12. 134 ICE_N2_15 = reform(ICE_N2_15, 12,nyr) 135 ICE_N2_15 = total(ICE_N2_15,2)/nyr 136 ICE_N2_15 = {arr:ICE_N2_15 * 1.e-12, unit : '10^12 m^2'} 137 ;ICE Extent in SOUTH Hemisphere 138 domdef, 0, jpi-1, -90, -30, /xindex 139 ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 140 ICE_S2 = grossemoyenne(ICE_S2.arr, 'xy', /integration, mask2d = masknp) 141 if jpt mod 12 ne 0 then stop 142 nyr=jpt/12. 143 ICE_S2 = reform(ICE_S2, 12,nyr) 144 ICE_S2 = total(ICE_S2,2)/nyr 145 ICE_S2 = {arr:ICE_S2 * 1.e-12, unit : '10^12 m^2'} 146 ;ICE Extent minus 15% in SOUTH Hemisphere 147 ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 148 msk = ICE_S2_15.arr gt 0.15 ; remove 0.15% for observations 149 ICE_S2_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 150 if jpt mod 12 ne 0 then stop 151 nyr=jpt/12. 152 ICE_S2_15 = reform(ICE_S2_15, 12,nyr) 153 ICE_S2_15 = total(ICE_S2_15,2)/nyr 154 ICE_S2_15 = {arr:ICE_S2_15 * 1.e-12, unit : '10^12 m^2'} 155 ; 156 ; time = tsave & IF n_elements(time) NE jpt THEN stop 53 157 54 158 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 55 159 56 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 57 pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 58 , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 59 60 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 61 pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 62 , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 160 161 title = 'Northern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+'!C'+' Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 162 jpt=12 163 time=julday(1,15,1900)+30*lindgen(12) 164 pltt, ICE_N.arr - ICE_N2.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 165 , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 166 pltt, ICE_N_15.arr - ICE_N2_15.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 167 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 168 ;SF pltt, vICE_area_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 169 ;SF , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 170 ;SF pltt, vICE_ext_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 171 ;SF , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 172 ; 173 ; title ='Southern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+' - '+'!C'+'Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 174 title ='Southern Hemisphere'+'!C' 175 pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE $ 176 , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 177 pltt, ICE_S_15.arr - ICE_S2_15.arr , 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 178 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 179 ;SF pltt, vICE_area_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 180 ;SF , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 181 ;SF pltt, vICE_ext_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 182 ;SF , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 183 ; 63 184 64 185 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png /> ' ] -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro
r2751 r4320 22 22 ; 23 23 ts_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 24 ts_SSH.arr = ts_SSH.arr * 100 & ts_SSH.unit = 'cm' 24 ;SF ts_SSH.arr = ts_SSH.arr * 100 & ts_SSH.unit = 'cm' 25 ts_SSH.arr = ts_SSH.arr & ts_SSH.unit = 'm' 25 26 title = prefix+' '+d1_d2+'!C'+blabla 26 pltt, ts_SSH, 't', -.1, .1, date1, date2, /REMPLI $ 27 , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 27 pltt, ts_SSH, 't', -12., 12., date1, date2, /REMPLI $ 28 ;SF , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 29 , small = [1, 2, 1], YTITLE = 'm', TITLE = title, /PORTRAIT, _extra = ex 28 30 29 31 IF prefix NE prefix2 THEN BEGIN … … 32 34 tsave = time 33 35 ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 34 ts_SSH2.arr = ts_SSH2.arr * 100 & ts_SSH2.unit = 'cm' 36 ;SF ts_SSH2.arr = ts_SSH2.arr * 100 & ts_SSH2.unit = 'cm' 37 ts_SSH2.arr = ts_SSH2.arr & ts_SSH2.unit = 'm' 35 38 time = tsave & IF n_elements(time) NE jpt THEN stop 36 39 37 40 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 38 pltt, ts_SSH.arr - ts_SSH2.arr, 't', -.1, .1, date1, date2, /REMPLI $ 39 , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 41 pltt, ts_SSH.arr - ts_SSH2.arr, 't', -10., 10., date1, date2, /REMPLI $ 42 ;SF , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 43 , color = 250, small = [1, 2, 2], YTITLE = 'm', TITLE = title, /noerase, _extra = ex 40 44 41 45 ENDIF -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro
r3610 r4320 47 47 48 48 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 49 pltt, ts_z, 'zt', -2., 2., inter = . 2, date1, date2, /REMPLI, style = 'so0so' $49 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $ 50 50 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 51 51 52 52 title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 53 pltt, ts_z, 'zt', -2., 2., inter = . 2, date1, date2, /REMPLI, style = 'so0so' $53 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $ 54 54 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 55 55 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r2751 r4320 24 24 ; reading variables 25 25 masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill) 26 ; climatologies 27 std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 28 std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 29 std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 30 std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 31 std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 32 std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato) 33 std_file_snow_arc = isafile(getenv('FILE_SNOW_ARC'), title = 'SNOW_ARC', iodir = std_iodir_climato) 34 std_file_snow_ant = isafile(getenv('FILE_SNOW_ANT'), title = 'SNOW_ANT', iodir = std_iodir_climato) 35 26 36 27 37 date1 = long(getenv('DATE1')) & date2 = long(getenv('DATE2')) … … 64 74 cnt = cnt+1 & blabla = 'Sea-Ice cover' 65 75 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 76 ; 77 cnt = cnt+1 & blabla = 'Sea-Ice Volume' 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 66 79 80 ;SF cnt = cnt+1 & blabla = 'SNOW Volume' 81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 82 ; 67 83 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] 68 84 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example1
r2753 r4320 18 18 idl_command=/Applications/itt/idl64/bin/idl 19 19 # 20 PS_DIR=$( pwd )/ps_ core_tmx21 PDF_DIR=$( pwd )/pdf_ core_tmx22 HTML_DIR=$( pwd )/html _core_tmx23 SAXO_DIR=/Users/s masson/SAXO_DIR20 PS_DIR=$( pwd )/ps_tag33_tag331 21 PDF_DIR=$( pwd )/pdf_tag33_tag331 22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/sflod/SAXO_DIR 24 24 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 28 # 29 29 #===================== Model GRID ===================== 30 30 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask31 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 32 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 33 # 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 41 45 # 42 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 48 52 READ_ONLY_FIRST_RECORD=0 # if 0 then read all records in files else read only the first reacord in each file 49 53 # 50 DATE1= 10101 ; DATE2=2000123154 DATE1=20010101 ; DATE2=21001231 51 55 # 52 VAR1_T=thetao ; V1T_PREF=tmx_2000 ; V1T_SUFF=_1Y_grid_T.nc 53 VAR1_S=so ; V1S_PREF=tmx_2000 ; V1S_SUFF=_1Y_grid_T.nc 54 VAR1_SSH=zos ; V1SSH_PREF=tmx_2000 ; V1SSH_SUFF=_1Y_grid_T.nc 55 VAR1_Q=qns+qsr ; V1Q_PREF=tmx_2000 ; V1Q_SUFF=_1Y_grid_T.nc 56 VAR1_EMP=wfo ; V1EMP_PREF=tmx_2000 ; V1EMP_SUFF=_1Y_grid_T.nc 57 VAR1_U=uoce_eff ; V1U_PREF=tmx_2000 ; V1U_SUFF=_1Y_grid_U.nc 58 VAR1_V=voce_eff ; V1V_PREF=tmx_2000 ; V1V_SUFF=_1Y_grid_V.nc 59 VAR1_ICE=sic ; V1ICE_PREF=tmx_2000 ; V1ICE_SUFF=_1Y_icemod.nc 56 VAR1_T=thetao ; V1T_PREF=tmx_2000 ; V1T_SUFF=_1Y_grid_T.nc 57 VAR1_S=so ; V1S_PREF=tmx_2000 ; V1S_SUFF=_1Y_grid_T.nc 58 VAR1_SSH=zos ; V1SSH_PREF=tmx_2000 ; V1SSH_SUFF=_1Y_grid_T.nc 59 VAR1_Q=qns+qsr ; V1Q_PREF=tmx_2000 ; V1Q_SUFF=_1Y_grid_T.nc 60 VAR1_EMP=wfo ; V1EMP_PREF=tmx_2000 ; V1EMP_SUFF=_1Y_grid_T.nc 61 VAR1_U=uoce_eff ; V1U_PREF=tmx_2000 ; V1U_SUFF=_1Y_grid_U.nc 62 VAR1_V=voce_eff ; V1V_PREF=tmx_2000 ; V1V_SUFF=_1Y_grid_V.nc 63 VAR1_ICE=sic ; V1ICE_PREF=tmx_2000 ; V1ICE_SUFF=_1M_icemod.nc 64 VAR1_Ithick=sit ; V1It_PREF=tmx_2000 ; V1It_SUFF=_1M_icemod.nc 65 VAR1_SNOW=isnowthi ; V1SNOW_PREF=tmx_2000 ; V1SNOW_SUFF=_1M_icemod.nc 60 66 # 61 67 #===================== EXP2 ===================== … … 63 69 DATE1_2=10101 ; DATE2_2=20001231 64 70 # 65 VAR2_T=thetao ; V2T_PREF=core2000 ; V2T_SUFF=_1Y_grid_T.nc 66 VAR2_S=so ; V2S_PREF=core2000 ; V2S_SUFF=_1Y_grid_T.nc 67 VAR2_SSH=zos ; V2SSH_PREF=core2000 ; V2SSH_SUFF=_1Y_grid_T.nc 68 VAR2_Q=qns+qsr ; V2Q_PREF=core2000 ; V2Q_SUFF=_1Y_grid_T.nc 69 VAR2_EMP=wfo ; V2EMP_PREF=core2000 ; V2EMP_SUFF=_1Y_grid_T.nc 70 VAR2_U=uoce_eff ; V2U_PREF=core2000 ; V2U_SUFF=_1Y_grid_U.nc 71 VAR2_V=voce_eff ; V2V_PREF=core2000 ; V2V_SUFF=_1Y_grid_V.nc 72 VAR2_ICE=sic ; V2ICE_PREF=core2000 ; V2ICE_SUFF=_1Y_icemod.nc 71 VAR2_T=thetao ; V2T_PREF=core2000 ; V2T_SUFF=_1Y_grid_T.nc 72 VAR2_S=so ; V2S_PREF=core2000 ; V2S_SUFF=_1Y_grid_T.nc 73 VAR2_SSH=zos ; V2SSH_PREF=core2000 ; V2SSH_SUFF=_1Y_grid_T.nc 74 VAR2_Q=qns+qsr ; V2Q_PREF=core2000 ; V2Q_SUFF=_1Y_grid_T.nc 75 VAR2_EMP=wfo ; V2EMP_PREF=core2000 ; V2EMP_SUFF=_1Y_grid_T.nc 76 VAR2_U=uoce_eff ; V2U_PREF=core2000 ; V2U_SUFF=_1Y_grid_U.nc 77 VAR2_V=voce_eff ; V2V_PREF=core2000 ; V2V_SUFF=_1Y_grid_V.nc 78 VAR2_ICE=sic ; V2ICE_PREF=core2000 ; V2ICE_SUFF=_1M_icemod.nc 79 VAR2_Ithick=sit ; V2It_PREF=core2000 ; V2It_SUFF=_1M_icemod.nc 80 VAR2_SNOW=isnowthi ; V2ICE_SNOW=core2000 ; V2SNOW_SUFF=_1M_icemod.nc 73 81 # 74 82 ######################### Export Variables ############################### … … 86 94 export FILE_FLUX VAR_FLUX 87 95 export FILE_MLD VAR_MLD 96 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 97 export FILE_SNOW VAR_SNOW_NH VAR_SNOW_SH VAR_SNOW_area_NH VAR_SNOW_area_SH 88 98 #===================== EXP1 ===================== 89 99 export READ_ONLY_FIRST_RECORD 90 100 # 91 export DATE1 DATE2 92 export VAR1_T V1T_PREF V1T_SUFF 93 export VAR1_S V1S_PREF V1S_SUFF 94 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 95 export VAR1_Q V1Q_PREF V1Q_SUFF 96 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 97 export VAR1_U V1U_PREF V1U_SUFF 98 export VAR1_V V1V_PREF V1V_SUFF 99 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 101 export DATE1 DATE2 102 export VAR1_T V1T_PREF V1T_SUFF 103 export VAR1_S V1S_PREF V1S_SUFF 104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 105 export VAR1_Q V1Q_PREF V1Q_SUFF 106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 107 export VAR1_U V1U_PREF V1U_SUFF 108 export VAR1_V V1V_PREF V1V_SUFF 109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 110 export VAR1_Ithick V1It_PREF V1It_SUFF 111 export VAR1_SNOW V1SNOW_PREF V1SNOW_SUFF 112 100 113 #===================== EXP2 ===================== 101 export DATE1_2 DATE2_2 102 export VAR2_T V2T_PREF V2T_SUFF 103 export VAR2_S V2S_PREF V2S_SUFF 104 export VAR2_SSH V2SSH_PREF V2SSH_SUFF 105 export VAR2_Q V2Q_PREF V2Q_SUFF 106 export VAR2_EMP V2EMP_PREF V2EMP_SUFF 107 export VAR2_U V2U_PREF V2U_SUFF 108 export VAR2_V V2V_PREF V2V_SUFF 109 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 114 export DATE1_2 DATE2_2 115 export VAR2_T V2T_PREF V2T_SUFF 116 export VAR2_S V2S_PREF V2S_SUFF 117 export VAR2_SSH V2SSH_PREF V2SSH_SUFF 118 export VAR2_Q V2Q_PREF V2Q_SUFF 119 export VAR2_EMP V2EMP_PREF V2EMP_SUFF 120 export VAR2_U V2U_PREF V2U_SUFF 121 export VAR2_V V2V_PREF V2V_SUFF 122 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 123 export VAR2_Ithick V2It_PREF V2It_SUFF 124 export VAR2_SNOW V2SNOW_PREF V2SNOW_SUFF 110 125 # -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example2
r2753 r4320 20 20 PS_DIR=$( pwd )/ps_tag33_tag331 21 21 PDF_DIR=$( pwd )/pdf_tag33_tag331 22 HTML_DIR=$( pwd )/html _tag33_tag33123 SAXO_DIR=/Users/s masson/SAXO_DIR22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/sflod/SAXO_DIR 24 24 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 28 # 29 29 #===================== Model GRID ===================== 30 30 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask31 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 32 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 33 # 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 41 45 # 42 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 50 54 DATE1=20010101 ; DATE2=21001231 51 55 # 52 VAR1_T=thetao ; V1T_PREF=tag33 ; V1T_SUFF=_1Y_grid_T.nc 53 VAR1_S=so ; V1S_PREF=tag33 ; V1S_SUFF=_1Y_grid_T.nc 54 VAR1_SSH=zos ; V1SSH_PREF=tag33 ; V1SSH_SUFF=_1Y_grid_T.nc 55 VAR1_Q=qt ; V1Q_PREF=tag33 ; V1Q_SUFF=_1Y_grid_T.nc 56 VAR1_EMP=wfo ; V1EMP_PREF=tag33 ; V1EMP_SUFF=_1Y_grid_T.nc 57 VAR1_U=uocetr_eff ; V1U_PREF=tag33 ; V1U_SUFF=_1Y_grid_U.nc 58 VAR1_V=vocetr_eff ; V1V_PREF=tag33 ; V1V_SUFF=_1Y_grid_V.nc 59 VAR1_ICE=sic ; V1ICE_PREF=tag33 ; V1ICE_SUFF=_1Y_icemod.nc 56 VAR1_T=thetao ; V1T_PREF=tag33 ; V1T_SUFF=_1Y_grid_T.nc 57 VAR1_S=so ; V1S_PREF=tag33 ; V1S_SUFF=_1Y_grid_T.nc 58 VAR1_SSH=zos ; V1SSH_PREF=tag33 ; V1SSH_SUFF=_1Y_grid_T.nc 59 VAR1_Q=qt ; V1Q_PREF=tag33 ; V1Q_SUFF=_1Y_grid_T.nc 60 VAR1_EMP=wfo ; V1EMP_PREF=tag33 ; V1EMP_SUFF=_1Y_grid_T.nc 61 VAR1_U=uocetr_eff ; V1U_PREF=tag33 ; V1U_SUFF=_1Y_grid_U.nc 62 VAR1_V=vocetr_eff ; V1V_PREF=tag33 ; V1V_SUFF=_1Y_grid_V.nc 63 VAR1_ICE=sic ; V1ICE_PREF=tag33 ; V1ICE_SUFF=_1M_icemod.nc 64 VAR1_Ithick=sit ; V1It_PREF=tag33 ; V1It_SUFF=_1M_icemod.nc 65 VAR1_SNOW=isnowthi ; V1SNOW_PREF=tag33 ; V1SNOW_SUFF=_1M_icemod.nc 60 66 # 61 67 #===================== EXP2 ===================== … … 63 69 DATE1_2=20010101 ; DATE2_2=21001231 64 70 # 65 VAR2_T=thetao ; V2T_PREF=trunk331 ; V2T_SUFF=_1Y_grid_T.nc 66 VAR2_S=so ; V2S_PREF=trunk331 ; V2S_SUFF=_1Y_grid_T.nc 67 VAR2_SSH=zos ; V2SSH_PREF=trunk331 ; V2SSH_SUFF=_1Y_grid_T.nc 68 VAR2_Q=qt ; V2Q_PREF=trunk331 ; V2Q_SUFF=_1Y_grid_T.nc 69 VAR2_EMP=wfo ; V2EMP_PREF=trunk331 ; V2EMP_SUFF=_1Y_grid_T.nc 70 VAR2_U=uocetr_eff ; V2U_PREF=trunk331 ; V2U_SUFF=_1Y_grid_U.nc 71 VAR2_V=vocetr_eff ; V2V_PREF=trunk331 ; V2V_SUFF=_1Y_grid_V.nc 72 VAR2_ICE=sic ; V2ICE_PREF=trunk331 ; V2ICE_SUFF=_1Y_icemod.nc 71 VAR2_T=thetao ; V2T_PREF=trunk331 ; V2T_SUFF=_1Y_grid_T.nc 72 VAR2_S=so ; V2S_PREF=trunk331 ; V2S_SUFF=_1Y_grid_T.nc 73 VAR2_SSH=zos ; V2SSH_PREF=trunk331 ; V2SSH_SUFF=_1Y_grid_T.nc 74 VAR2_Q=qt ; V2Q_PREF=trunk331 ; V2Q_SUFF=_1Y_grid_T.nc 75 VAR2_EMP=wfo ; V2EMP_PREF=trunk331 ; V2EMP_SUFF=_1Y_grid_T.nc 76 VAR2_U=uocetr_eff ; V2U_PREF=trunk331 ; V2U_SUFF=_1Y_grid_U.nc 77 VAR2_V=vocetr_eff ; V2V_PREF=trunk331 ; V2V_SUFF=_1Y_grid_V.nc 78 VAR2_ICE=sic ; V2ICE_PREF=trunk331 ; V2ICE_SUFF=_1M_icemod.nc 79 VAR2_Ithick=sit ; V2It_PREF=trunk331 ; V2It_SUFF=_1M_icemod.nc 80 VAR2_SNOW=isnowthi ; V2ICE_SNOW=trunk331 ; V2SNOW_SUFF=_1M_icemod.nc 73 81 # 74 82 ######################### Export Variables ############################### … … 86 94 export FILE_FLUX VAR_FLUX 87 95 export FILE_MLD VAR_MLD 96 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 97 export FILE_SNOW VAR_SNOW_NH VAR_SNOW_SH VAR_SNOW_area_NH VAR_SNOW_area_SH 88 98 #===================== EXP1 ===================== 89 99 export READ_ONLY_FIRST_RECORD 90 100 # 91 export DATE1 DATE2 92 export VAR1_T V1T_PREF V1T_SUFF 93 export VAR1_S V1S_PREF V1S_SUFF 94 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 95 export VAR1_Q V1Q_PREF V1Q_SUFF 96 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 97 export VAR1_U V1U_PREF V1U_SUFF 98 export VAR1_V V1V_PREF V1V_SUFF 99 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 101 export DATE1 DATE2 102 export VAR1_T V1T_PREF V1T_SUFF 103 export VAR1_S V1S_PREF V1S_SUFF 104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 105 export VAR1_Q V1Q_PREF V1Q_SUFF 106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 107 export VAR1_U V1U_PREF V1U_SUFF 108 export VAR1_V V1V_PREF V1V_SUFF 109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 110 export VAR1_Ithick V1It_PREF V1It_SUFF 111 export VAR1_SNOW V1SNOW_PREF V1SNOW_SUFF 112 100 113 #===================== EXP2 ===================== 101 export DATE1_2 DATE2_2 102 export VAR2_T V2T_PREF V2T_SUFF 103 export VAR2_S V2S_PREF V2S_SUFF 104 export VAR2_SSH V2SSH_PREF V2SSH_SUFF 105 export VAR2_Q V2Q_PREF V2Q_SUFF 106 export VAR2_EMP V2EMP_PREF V2EMP_SUFF 107 export VAR2_U V2U_PREF V2U_SUFF 108 export VAR2_V V2V_PREF V2V_SUFF 109 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 114 export DATE1_2 DATE2_2 115 export VAR2_T V2T_PREF V2T_SUFF 116 export VAR2_S V2S_PREF V2S_SUFF 117 export VAR2_SSH V2SSH_PREF V2SSH_SUFF 118 export VAR2_Q V2Q_PREF V2Q_SUFF 119 export VAR2_EMP V2EMP_PREF V2EMP_SUFF 120 export VAR2_U V2U_PREF V2U_SUFF 121 export VAR2_V V2V_PREF V2V_SUFF 122 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 123 export VAR2_Ithick V2It_PREF V2It_SUFF 124 export VAR2_SNOW V2SNOW_PREF V2SNOW_SUFF 110 125 #
Note: See TracChangeset
for help on using the changeset viewer.