Changeset 4559
- Timestamp:
- 2014-03-21T10:22:46+01:00 (10 years ago)
- Location:
- trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 17 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro
r4519 r4559 1 pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceAge, Iage1, Iage2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 50 50 ; 51 51 title = var+'!C'+std_file1_I 52 53 Iage1.arr = Iage1.arr * ( Ifra1.arr gt 0.15 ) 54 52 55 plt, Iage1.arr - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 53 56 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div $ … … 56 59 if std_file1_I NE std_file2_I then begin 57 60 title = title + std_file2_I 61 62 Iage2.arr = Iage2.arr * ( Ifra2.arr gt 0.15 ) 63 58 64 plt, Iage1.arr - Iage2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 59 65 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r4518 r4559 32 32 ; 33 33 title = var+'!C'+std_file1_I 34 plt, (Ifra1.arr < 1. ) - 1.E-04, MIN = 0., MAX = 1., INTER = .1, /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 34 Ifra1.arr = Ifra1.arr * ( Ifra1.arr gt 0.15 ) 35 plt, Ifra1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = .1, /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 35 36 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 36 37 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex … … 38 39 if std_file1_I NE std_file2_I then begin 39 40 title = title + std_file2_I 41 Ifra2.arr = Ifra2.arr * ( Ifra2.arr gt 0.15 ) 40 42 plt, Ifra1.arr - Ifra2.arr, MIN = -1, MAX = 1, INTER = .1, STYLE = 'so0so', format = '(f4.1)' $ 41 43 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro
r4519 r4559 1 pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceSal, Isal1, Isal2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 32 32 ; 33 33 title = var+'!C'+std_file1_I 34 plt, Isal1.arr - 1.E-04, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 34 35 Isal1.arr = Isal1.arr * ( Ifra1.arr gt 0.15 ) 36 plt, Isal1.arr - 1.E-04, MIN = 0., MAX = 10., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 35 37 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 36 38 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex … … 38 40 if std_file1_I NE std_file2_I then begin 39 41 title = title + std_file2_I 40 plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 42 43 Isal2.arr = Isal2.arr * ( Ifra2.arr gt 0.15 ) 44 45 plt, Isal1.arr - Isal2.arr, MIN = -2., MAX = 2., INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 41 46 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 42 47 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro
r4518 r4559 1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_IceThick, Ithi1, Ithi2, Ifra1, Ifra2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 50 50 ; 51 51 title = var+'!C'+std_file1_I 52 53 Ithi1.arr = Ithi1.arr * ( Ifra1.arr gt 0.15 ) 54 52 55 plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 53 56 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div $ … … 56 59 if std_file1_I NE std_file2_I then begin 57 60 title = title + std_file2_I 61 62 Ithi2.arr = Ithi2.arr * ( Ifra2.arr gt 0.15 ) 63 58 64 plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 59 65 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
r4518 r4559 1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex1 pro std_plot_SnowThick, Isnow1, Isnow2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 2 2 3 3 compile_opt idl2, strictarrsubs … … 50 50 ; 51 51 title = var+'!C'+std_file1_I 52 53 Isnow1.arr = Isnow1.arr * ( Ifra1.arr gt 0.15 ) 54 52 55 plt, (Isnow1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 53 56 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ … … 56 59 if std_file1_I NE std_file2_I then begin 57 60 title = title + std_file2_I 61 62 Isnow2.arr = Isnow2.arr * ( Ifra2.arr gt 0.15 ) 63 58 64 plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 59 65 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4518 r4559 36 36 ; 37 37 allrec = 1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD'))) 38 38 39 ; Output run experience1 39 40 std_file1_T = isafile(getenv('FILE1_T'), title = 'exp1 grid T input file', iodir = std_iodir_data) … … 192 193 ENDELSE 193 194 ; 194 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )195 caldat, time, mm196 febr = where(mm EQ 2, cnt)197 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit}198 march = where(mm EQ 3, cnt)199 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit}200 sept = where(mm EQ 9, cnt)201 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}202 undefine, Iage_1203 ;204 IF std_file2_I NE std_file1_I THEN BEGIN205 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )206 caldat, time, mm207 febr = where(mm EQ 2, cnt)208 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}209 march = where(mm EQ 3, cnt)210 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}211 sept = where(mm EQ 9, cnt)212 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}213 undefine, Iage_2214 ENDIF ELSE BEGIN215 Iage_febr_2 = {arr:-1}216 Iage_march_2 = {arr:-1}217 Iage_sept_2 = {arr:-1}218 ENDELSE219 ;220 195 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 221 196 caldat, time, mm 197 jan = where(mm EQ 1, cnt) 198 Ifra_jan_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 222 199 febr = where(mm EQ 2, cnt) 223 200 Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 224 201 march = where(mm EQ 3, cnt) 225 202 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 203 april = where(mm EQ 4, cnt) 204 Ifra_april_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 226 205 sept = where(mm EQ 9, cnt) 227 206 Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} … … 231 210 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 232 211 caldat, time, mm 212 jan = where(mm EQ 1, cnt) 213 Ifra_jan_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 233 214 febr = where(mm EQ 2, cnt) 234 215 Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 235 216 march = where(mm EQ 3, cnt) 236 217 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 218 april = where(mm EQ 4, cnt) 219 Ifra_april_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 237 220 sept = where(mm EQ 9, cnt) 238 221 Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 239 222 undefine, Ifra_2 240 223 ENDIF ELSE BEGIN 224 Ifra_jan_2 = {arr:-1} 241 225 Ifra_febr_2 = {arr:-1} 242 226 Ifra_march_2 = {arr:-1} 227 Ifra_april_2 = {arr:-1} 243 228 Ifra_sept_2 = {arr:-1} 244 229 ENDELSE … … 271 256 ; 272 257 Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 273 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 274 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 275 help, Ifra_1 276 help, mask 258 ;SF ready for mask : remove 0.15% for observations 259 ;SF 260 ;SF Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 261 ;SF msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 277 262 caldat, time, mm 278 263 march = where(mm EQ 3, cnt) 279 264 febr = where(mm EQ 2, cnt) 280 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}281 Isal_1.arr = Isal_1.arr * msk282 help, Isal_1283 265 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 284 help, Isal_march_1 266 ;SF 267 ;SF Isal_1.arr = Isal_1.arr * msk 268 ;SF 269 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 285 270 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 286 271 sept = where(mm EQ 9, cnt) … … 289 274 ; 290 275 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 291 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )292 msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations293 276 IF std_file2_I NE std_file1_I THEN BEGIN 294 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )295 Isal_2.arr = Isal_2.arr * msk296 277 caldat, time, mm 297 278 march = where(mm EQ 3, cnt) … … 307 288 Isal_sept_2 = {arr:-1} 308 289 ENDELSE 290 ; 291 IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I ) 292 IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I ) 293 Ivelo_1 = read_ncdf(getenv('VAR1_Ivelo'), allrecords = allrec, filename = std_file1_I ) 294 caldat, time, mm 295 febr = where(mm EQ 2, cnt) 296 IvelU_febr_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 297 IvelV_febr_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 298 Ivelo_febr_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'} 299 march = where(mm EQ 3, cnt) 300 IvelU_march_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 301 IvelV_march_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 302 Ivelo_march_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'} 303 sept = where(mm EQ 9, cnt) 304 IvelU_sept_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 305 IvelV_sept_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 306 Ivelo_sept_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'} 307 ; 308 IF std_file2_I NE std_file1_I THEN BEGIN 309 IvelU_2 = read_ncdf(getenv('VAR2_IvelU'), allrecords = allrec, filename = std_file2_I ) 310 IvelV_2 = read_ncdf(getenv('VAR2_IvelV'), allrecords = allrec, filename = std_file2_I ) 311 Ivelo_2 = read_ncdf(getenv('VAR2_Ivelo'), allrecords = allrec, filename = std_file2_I ) 312 caldat, time, mm 313 febr = where(mm EQ 2, cnt) 314 IvelU_febr_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 315 IvelV_febr_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelV_2.unit} 316 Ivelo_febr_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 317 march = where(mm EQ 3, cnt) 318 IvelU_march_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 319 IvelV_march_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelV_2.unit} 320 Ivelo_march_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 321 sept = where(mm EQ 9, cnt) 322 IvelU_sept_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 323 IvelV_sept_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelV_2.unit} 324 Ivelo_sept_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 325 undefine, Ifra_2 326 ENDIF ELSE BEGIN 327 IvelU_febr_2 = {arr:-1} 328 IvelV_febr_2 = {arr:-1} 329 Ivelo_febr_2 = {arr:-1} 330 IvelU_march_2 = {arr:-1} 331 IvelV_march_2 = {arr:-1} 332 Ivelo_march_2 = {arr:-1} 333 IvelU_sept_2 = {arr:-1} 334 IvelV_sept_2 = {arr:-1} 335 Ivelo_sept_2 = {arr:-1} 336 ENDELSE 337 ; 338 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 339 caldat, time, mm 340 febr = where(mm EQ 2, cnt) 341 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 342 march = where(mm EQ 3, cnt) 343 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 344 sept = where(mm EQ 9, cnt) 345 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 346 undefine, Iage_1 347 ; 348 IF std_file2_I NE std_file1_I THEN BEGIN 349 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 350 caldat, time, mm 351 febr = where(mm EQ 2, cnt) 352 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 353 march = where(mm EQ 3, cnt) 354 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 355 sept = where(mm EQ 9, cnt) 356 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 357 undefine, Iage_2 358 ENDIF ELSE BEGIN 359 Iage_febr_2 = {arr:-1} 360 Iage_march_2 = {arr:-1} 361 Iage_sept_2 = {arr:-1} 362 ENDELSE 309 363 ; 310 364 jpt = 1 … … 442 496 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 443 497 ; 444 cnt = cnt+1 & blabla = 'Mediterranean water at lat=40 °N'498 cnt = cnt+1 & blabla = 'Mediterranean water at lat=40¡ N' 445 499 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 446 500 ; 447 cnt = cnt+1 & blabla = 'Mediterranean water at lat=38 °N'501 cnt = cnt+1 & blabla = 'Mediterranean water at lat=38¡ N' 448 502 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 449 503 ; 450 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 451 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 452 ; 453 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 454 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 455 ; 456 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 457 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 458 ; 459 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' 460 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 461 ; 462 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 463 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 464 ; 465 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 466 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 467 ; 468 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 469 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 470 ; 471 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 472 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 473 ; 504 ;; 505 ; all plot are done for ice fraction > 0.15% 506 ;; 474 507 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: MARCH' 475 508 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex … … 484 517 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 485 518 ; 519 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 520 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, Ifra_jan_1, Ifra_jan_2, /ARC, /JAN, _extra = ex 521 ; 522 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 523 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, Ifra_april_1, Ifra_april_2, /ARC, /APRIL, _extra = ex 524 ; 525 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 526 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, Ifra_april_1, Ifra_april_2, /ANT, /APRIL, _extra = ex 527 ; 528 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' 529 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 530 ; 486 531 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: MARCH' 487 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex532 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 488 533 ; 489 534 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 490 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex535 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 491 536 ; 492 537 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 493 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex538 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 494 539 ; 495 540 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 496 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex541 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 497 542 ; 498 543 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: MARCH' 499 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex544 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 500 545 ; 501 546 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: SEPT' 502 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex547 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 503 548 ; 504 549 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: FEBRUARY' 505 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex550 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 506 551 ; 507 552 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: SEPT' 508 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 509 ; 553 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 554 ; 555 cnt = cnt+1 & blabla = 'Arctic Ice Velocity: FEBRUARY' 556 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_febr_1, IvelU_febr_2, IvelV_febr_1, IvelV_febr_2, Ivelo_febr_1, Ivelo_febr_2, /ARC, /FEBR, _extra = ex 557 ; 558 cnt = cnt+1 & blabla = 'Arctic Ice Velocity: SEPT' 559 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ARC, /SEPT, _extra = ex 560 ; 561 cnt = cnt+1 & blabla = 'Antartic Ice Velocity: MARCH' 562 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_march_1, IvelU_march_2, IvelV_march_1, IvelV_march_2, Ivelo_march_1, Ivelo_march_2, /ANT, /MARCH, _extra = ex 563 ; 564 cnt = cnt+1 & blabla = 'Antartic Ice Velocity: SEPT' 565 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT, /SEPT, _extra = ex 566 ; 567 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 568 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 569 ; 570 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 571 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 572 ; 573 cnt = cnt+1 & blabla = 'Antarctic Ice Age: FEBR' 574 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 575 ; 576 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 577 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 578 ; 510 579 511 580 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro
r4518 r4559 273 273 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 274 274 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 275 help, Ifra_1276 help, mask277 275 caldat, time, mm 278 276 march = where(mm EQ 3, cnt) … … 280 278 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 281 279 Isal_1.arr = Isal_1.arr * msk 282 help, Isal_1283 280 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 284 help, Isal_march_1285 281 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 286 282 sept = where(mm EQ 9, cnt) -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
r4518 r4559 289 289 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 290 290 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 291 help, Ifra_1292 help, mask293 291 caldat, time, mm 294 292 march = where(mm EQ 3, cnt) … … 296 294 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 297 295 Isal_1.arr = Isal_1.arr * msk 298 help, Isal_1299 296 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 300 help, Isal_march_1301 297 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 302 298 sept = where(mm EQ 9, cnt) -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
r4516 r4559 72 72 VAR1_Isal=sisali 73 73 VAR1_Iage=siages 74 VAR1_IvelU=sivelu 75 VAR1_IvelV=sivelv 76 VAR1_Ivelo=sivelo 74 77 # 75 78 #===================== EXP2 ===================== … … 92 95 VAR2_Isal=sisali 93 96 VAR2_Iage=siages 97 VAR2_IvelU=sivelu 98 VAR2_IvelV=sivelv 99 VAR2_Ivelo=sivelo 94 100 # 95 101 ######################### Export Variables ############################### … … 115 121 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 116 122 export VAR1_U VAR1_V 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 118 124 #===================== EXP2 ===================== 119 125 export FILE2_T FILE2_U FILE2_V FILE2_I 120 126 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 121 127 export VAR2_U VAR2_V 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 123 129 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D
r4514 r4559 72 72 VAR1_Isal=iicesali 73 73 VAR1_Iage=iiceages 74 VAR1_IvelU=iicevelu 75 VAR1_IvelV=iicevelv 76 VAR1_Ivelo=iicevelo 74 77 # 75 78 #===================== EXP2 ===================== … … 92 95 VAR2_Isal=iicesali 93 96 VAR2_Iage=iiceages 97 VAR2_IvelU=iicevelu 98 VAR2_IvelV=iicevelv 99 VAR2_Ivelo=iicevelo 94 100 # 95 101 ######################### Export Variables ############################### … … 115 121 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 116 122 export VAR1_U VAR1_V 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 118 124 #===================== EXP2 ===================== 119 125 export FILE2_T FILE2_U FILE2_V FILE2_I 120 126 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 121 127 export VAR2_U VAR2_V 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 123 129 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names
r4514 r4559 72 72 VAR1_Isal=sisali 73 73 VAR1_Iage=siages 74 VAR1_IvelU=sivelu 75 VAR1_IvelV=sivelv 76 VAR1_Ivelo=sivelo 74 77 # 75 78 #===================== EXP2 ===================== … … 92 95 VAR2_Isal=sisali 93 96 VAR2_Iage=siages 97 VAR2_IvelU=sivelu 98 VAR2_IvelV=sivelv 99 VAR2_Ivelo=sivelo 94 100 # 95 101 ######################### Export Variables ############################### … … 115 121 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 116 122 export VAR1_U VAR1_V 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 118 124 #===================== EXP2 ===================== 119 125 export FILE2_T FILE2_U FILE2_V FILE2_I 120 126 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 121 127 export VAR2_U VAR2_V 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 123 129 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names
r4514 r4559 72 72 VAR1_Isal=iicesali 73 73 VAR1_Iage=iiceages 74 VAR1_IvelU=iicevelu 75 VAR1_IvelV=iicevelv 76 VAR1_Ivelo=iicevelo 74 77 # 75 78 #===================== EXP2 ===================== … … 94 97 VAR2_Isal=iicesali 95 98 VAR2_Iage=iiceages 99 VAR2_IvelU=iicevelu 100 VAR2_IvelV=iicevelv 101 VAR2_Ivelo=iicevelo 96 102 # 97 103 ######################### Export Variables ############################### … … 117 123 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 118 124 export VAR1_U VAR1_V 119 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 125 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 120 126 #===================== EXP2 ===================== 121 127 export FILE2_T FILE2_U FILE2_V FILE2_I 122 128 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 123 129 export VAR2_U VAR2_V 124 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 130 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 125 131 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1
r4514 r4559 72 72 VAR1_Isal=iicesali 73 73 VAR1_Iage=iiceages 74 VAR1_IvelU=iicevelu 75 VAR1_IvelV=iicevelv 76 VAR1_Ivelo=iicevelo 74 77 # 75 78 #===================== EXP2 ===================== … … 92 95 VAR2_Isal=iicesali 93 96 VAR2_Iage=iiceages 97 VAR2_IvelU=iicevelu 98 VAR2_IvelV=iicevelv 99 VAR2_Ivelo=iicevelo 94 100 # 95 101 ######################### Export Variables ############################### … … 115 121 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 116 122 export VAR1_U VAR1_V 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 118 124 #===================== EXP2 ===================== 119 125 export FILE2_T FILE2_U FILE2_V FILE2_I 120 126 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 121 127 export VAR2_U VAR2_V 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 123 129 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2
r4514 r4559 72 72 VAR1_Isal=iicesali 73 73 VAR1_Iage=iiceages 74 VAR1_IvelU=iicevelu 75 VAR1_IvelV=iicevelv 76 VAR1_Ivelo=iicevelo 74 77 # 75 78 #===================== EXP2 ===================== … … 92 95 VAR1_Isal=iicesali 93 96 VAR1_Iage=iiceages 97 VAR2_IvelU=iicevelu 98 VAR2_IvelV=iicevelv 99 VAR2_Ivelo=iicevelo 94 100 # 95 101 ######################### Export Variables ############################### … … 115 121 export VAR1_T VAR1_S VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 116 122 export VAR1_U VAR1_V 117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 118 124 #===================== EXP2 ===================== 119 125 export FILE2_T FILE2_U FILE2_V FILE2_I 120 126 export VAR2_T VAR2_S VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 121 127 export VAR2_U VAR2_V 122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 123 129 # -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vel.pro
r4514 r4559 17 17 @common 18 18 @std_common 19 19 ; get exp1 info 20 vICE1 = getenv('VAR1_ICE') & prefix = getenv('V1ICE_PREF') & suffix = getenv('V1ICE_SUFF') 21 ; get exp2 info 22 vICE2 = getenv('VAR2_ICE') & prefix2 = getenv('V2ICE_PREF') & suffix2 = getenv('V2ICE_SUFF') 20 23 ; get exp1 info 21 24 vICE_vel_1 = getenv('VAR1_Ivel') & prefix = getenv('V1Iv_PREF') & suffix = getenv('V1Iv_SUFF') … … 34 37 ; ICE Velocity in NORTH Hemisphere 35 38 domdef, 0, jpi-1, 30, 90, /xindex 36 Velo_N = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward 39 Velo_N = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) ;!! warning positive northward 40 ICE_N_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) 41 print, 'N15', max(ICE_N_15) 37 42 38 ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54 Array[180, 54, 120]39 ;;;;;; with the domain dimensions [jpi/nx, jpj/ny, jpk/nz, jpt] = [180/180, 149/50, 31/31, 12]40 43 41 ICE_velo_N = grossemoyenne(Velo_N, 'xy', /integration, mask2d = masknp) 42 help, ICE_velo_N 43 print, 'max ice velo', max(ICE_velo_N) 44 ICE_N_15[where(ICE_N_15 lt 0.15)] = 0. 45 46 ICE_velo_N = grossemoyenne( (Velo_N < 1.e10) * (ICE_N_15 < 1.e10), 'xy',/integration, mask2d = masknp) 47 ICE_N_15 = grossemoyenne(ICE_N_15, 'xy',/integration, mask2d = masknp) 48 49 ICE_velo_N = ICE_velo_N / ICE_N_15 44 50 45 51 if jpt mod 12 ne 0 then stop … … 47 53 ICE_velo_N = reform(ICE_velo_N, 12, nyr) 48 54 ICE_velo_N = total(ICE_velo_N,2)/nyr 49 ICE_velo_N = {arr:ICE_velo_N * 1.e-19 * 86400 * 365, unit : 'm/year'} 50 help, ICE_velo_N 51 print, 'max ice velo', max(ICE_velo_N.arr) 52 55 ICE_velo_N = {arr:ICE_velo_N , unit : 'm/s'} 56 57 53 58 ;ICE Velocity in SOUTH Hemisphere 54 59 domdef, 0, jpi-1, -90, -30, /xindex 55 60 Velo_S = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward 61 56 62 57 63 … … 59 65 jpt=12 60 66 time=julday(1,15,1900)+30*lindgen(12) 61 pltt, ICE_velo_N, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = 60000000. , XGRIDSTYLE = 1 $ 67 68 69 pltt, ICE_velo_N, 't', /REMPLI, /PORTRAIT, XGRIDSTYLE = 1 $ 62 70 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 63 71 ; 64 72 title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Ice Velocity (Black SOLID simulation)' 65 pltt, Velo_S, 't', /REMPLI, MIN = 0., MAX = 60000000.,/NOERASE, XGRIDSTYLE = 1 $73 pltt, Velo_S, 't', /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 66 74 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 67 75 ; … … 77 85 domdef, 0, jpi-1, 30, 90, /xindex 78 86 Velo_N_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct) 79 help, Velo_N_280 87 81 88 ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54 Array[180, 54, 120] … … 83 90 84 91 ICE_velo_N_2 = grossemoyenne(Velo_N_2, 'xy', /integration, mask2d = masknp) 85 help, ICE_velo_N_286 92 print, 'max ice velo', max(ICE_velo_N_2) 87 93 … … 91 97 ICE_velo_N_2 = total(ICE_velo_N_2,2)/nyr 92 98 ICE_velo_N_2 = {arr:ICE_velo_N_2 * 86400 * 365, unit : 'm/year'} 93 help, ICE_velo_N_294 99 print, 'max ice velo', max(ICE_velo_N_2.arr) 95 100 … … 98 103 domdef, 0, jpi-1, -90, -30, /xindex 99 104 Velo_S_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct) 100 ;SF help, Velo_S_2101 105 102 106 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 103 107 108 jpt=12 109 time=julday(1,15,1900)+30*lindgen(12) 110 104 111 title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) ' 105 pltt, ICE_velo_N , 't', MIN = 0., MAX = 60000000. , date1, date2,/REMPLI, XGRIDSTYLE = 1, window = 2 $112 pltt, ICE_velo_N , 't', /REMPLI, XGRIDSTYLE = 1, window = 2 $ 106 113 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 107 pltt, ICE_velo_N_2, 't', date1, date2,/REMPLI $114 pltt, ICE_velo_N_2, 't', /REMPLI $ 108 115 , /ov1d, color = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 109 116 110 117 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) ' 111 pltt, Velo_S , 't', MIN = -.04, MAX = .04, date1, date2,/REMPLI, XGRIDSTYLE = 1 $118 pltt, Velo_S , 't', /REMPLI, XGRIDSTYLE = 1 $ 112 119 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 113 pltt, Velo_S_2, 't', date1, date2,/REMPLI $120 pltt, Velo_S_2, 't', /REMPLI $ 114 121 , /ov1d, color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 115 122 -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r4514 r4559 78 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 79 79 ; 80 ;SFcnt = cnt+1 & blabla = 'ICE Velocity '81 ;SFIF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex80 cnt = cnt+1 & blabla = 'ICE Velocity ' 81 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 82 82 ; 83 83 cnt = cnt+1 & blabla = 'Sea-Ice Volume Export at Fram Strait' -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh
r4516 r4559 18 18 idl_command=/Applications/itt/idl64/bin/idl 19 19 # 20 PS_DIR=$( pwd )/ts_heatbug -sbcmodMV_ps21 PDF_DIR=$( pwd )/ts_heatbug -sbcmodMV_pdf20 PS_DIR=$( pwd )/ts_heatbug_ps 21 PDF_DIR=$( pwd )/ts_heatbug_pdf 22 22 HTML_DIR=$( pwd )/html 23 23 SAXO_DIR=/Users/sflod/SAXO_DIR … … 53 53 # 54 54 #DATE1=20010101 ; DATE2=25001231 55 DATE1=20 010101 ; DATE2=2050123155 DATE1=20310101 ; DATE2=20401231 56 56 #DATE1=00010101 ; DATE2=00101231 57 57 # … … 72 72 # 73 73 #DATE1_2=20010101 ; DATE2_2=20101231 74 DATE1_2=20 010101 ; DATE2_2=2050123174 DATE1_2=20310101 ; DATE2_2=20401231 75 75 #DATE1_2=20010101 ; DATE2_2=25001231 76 76 # 77 VAR2_T=thetao ; V2T_PREF= sbcmodMV; V2T_SUFF=_1Y_grid_T.nc78 VAR2_S=so ; V2S_PREF= sbcmodMV; V2S_SUFF=_1Y_grid_T.nc79 VAR2_SSH=zos ; V2SSH_PREF= sbcmodMV; V2SSH_SUFF=_1Y_grid_T.nc80 VAR2_Q=tohfls ; V2Q_PREF= sbcmodMV; V2Q_SUFF=_1Y_grid_T.nc81 VAR2_EMP=wfo ; V2EMP_PREF= sbcmodMV; V2EMP_SUFF=_1Y_grid_T.nc82 VAR2_U=uocetr_eff ; V2U_PREF= sbcmodMV; V2U_SUFF=_1Y_grid_U.nc83 VAR2_V=vocetr_eff ; V2V_PREF= sbcmodMV; V2V_SUFF=_1Y_grid_V.nc84 VAR2_ICE=siconc ; V2ICE_PREF= sbcmodMV; V2ICE_SUFF=_1M_icemod.nc85 VAR2_Ithick=sithic ; V2It_PREF= sbcmodMV; V2It_SUFF=_1M_icemod.nc86 VAR2_SNOW=sndept ; V2SNOW_PREF= sbcmodMV; V2SNOW_SUFF=_1M_icemod.nc87 VAR2_IvelV=sivelv ; V2IvV_PREF= sbcmodMV; V2IvV_SUFF=_1M_icemod.nc88 VAR2_Ivel=sivelo ; V2Iv_PREF= sbcmodMV; V2Iv_SUFF=_1M_icemod.nc77 VAR2_T=thetao ; V2T_PREF=heatbug ; V2T_SUFF=_1Y_grid_T.nc 78 VAR2_S=so ; V2S_PREF=heatbug ; V2S_SUFF=_1Y_grid_T.nc 79 VAR2_SSH=zos ; V2SSH_PREF=heatbug ; V2SSH_SUFF=_1Y_grid_T.nc 80 VAR2_Q=tohfls ; V2Q_PREF=heatbug ; V2Q_SUFF=_1Y_grid_T.nc 81 VAR2_EMP=wfo ; V2EMP_PREF=heatbug ; V2EMP_SUFF=_1Y_grid_T.nc 82 VAR2_U=uocetr_eff ; V2U_PREF=heatbug ; V2U_SUFF=_1Y_grid_U.nc 83 VAR2_V=vocetr_eff ; V2V_PREF=heatbug ; V2V_SUFF=_1Y_grid_V.nc 84 VAR2_ICE=siconc ; V2ICE_PREF=heatbug ; V2ICE_SUFF=_1M_icemod.nc 85 VAR2_Ithick=sithic ; V2It_PREF=heatbug ; V2It_SUFF=_1M_icemod.nc 86 VAR2_SNOW=sndept ; V2SNOW_PREF=heatbug ; V2SNOW_SUFF=_1M_icemod.nc 87 VAR2_IvelV=sivelv ; V2IvV_PREF=heatbug ; V2IvV_SUFF=_1M_icemod.nc 88 VAR2_Ivel=sivelo ; V2Iv_PREF=heatbug ; V2Iv_SUFF=_1M_icemod.nc 89 89 # 90 90 ######################### Export Variables ###############################
Note: See TracChangeset
for help on using the changeset viewer.