Changeset 4514 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
 Timestamp:
 20140225T12:41:51+01:00 (7 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
r4320 r4514 203 203 ENDELSE 204 204 ; 205 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 206 caldat, time, mm 207 march = where(mm EQ 3, cnt) 208 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 209 sept = where(mm EQ 9, cnt) 210 Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 211 undefine, Iage_1 212 ; 213 IF std_file2_I NE std_file1_I THEN BEGIN 214 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 215 caldat, time, mm 216 march = where(mm EQ 3, cnt) 217 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 218 sept = where(mm EQ 9, cnt) 219 Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 220 undefine, Iage_2 221 ENDIF ELSE BEGIN 222 Iage_march_2 = {arr:1} 223 Iage_sept_2 = {arr:1} 224 ENDELSE 225 ; 205 226 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 206 227 ;; help, Ifra_1.arr, jpi,jpj … … 225 246 ENDELSE 226 247 ; 227 ;; SF commentato perché non funziona snow 228 ;; SF Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 229 ;; SF help, Isnow_1.arr, jpi,jpj 230 ;; SF caldat, time, mm 231 ;; SF april = where(mm EQ 4, cnt) 232 ;; SF Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 233 ;; SF sept = where(mm EQ 9, cnt) 234 ;; SF Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 235 ;; SF ;undefine, Isnow_1 236 ;; SF ; 237 ;; SF IF std_file2_I NE std_file1_I THEN BEGIN 238 ;; SF Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 239 ;; SF caldat, time, mm 240 ;; SF april = where(mm EQ 4, cnt) 241 ;; SF Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 242 ;; SF sept = where(mm EQ 9, cnt) 243 ;; SF Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 244 ;; SF ; undefine, Isnow_2 245 ;; SF ENDIF ELSE BEGIN 246 ;; SF Isnow_april_2 = {arr:1} 247 ;; SF Isnow_sept_2 = {arr:1} 248 ;; SF ENDELSE 248 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 249 caldat, time, mm 250 march = where(mm EQ 3, cnt) 251 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 252 sept = where(mm EQ 9, cnt) 253 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 254 ;undefine, Isnow_1 255 ; 256 IF std_file2_I NE std_file1_I THEN BEGIN 257 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 258 caldat, time, mm 259 march = where(mm EQ 3, cnt) 260 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 261 sept = where(mm EQ 9, cnt) 262 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 263 ; undefine, Isnow_2 264 ENDIF ELSE BEGIN 265 Isnow_march_2 = {arr:1} 266 Isnow_sept_2 = {arr:1} 267 ENDELSE 268 ; 269 Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 270 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 271 msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 272 help, Ifra_1 273 help, mask 274 caldat, time, mm 275 march = where(mm EQ 3, cnt) 276 febr = where(mm EQ 2, cnt) 277 ;SF Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 278 Isal_1.arr = Isal_1.arr * msk 279 help, Isal_1 280 Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 281 help, Isal_march_1 282 Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 283 sept = where(mm EQ 9, cnt) 284 Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 285 ;undefine, Isal_1 286 ; 287 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 288 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 289 msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 290 IF std_file2_I NE std_file1_I THEN BEGIN 291 Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 292 Isal_2.arr = Isal_2.arr * msk 293 caldat, time, mm 294 march = where(mm EQ 3, cnt) 295 febr = where(mm EQ 2, cnt) 296 Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 297 Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 298 sept = where(mm EQ 9, cnt) 299 Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 300 ; undefine, Isal_2 301 ENDIF ELSE BEGIN 302 Isal_febr_2 = {arr:1} 303 Isal_march_2 = {arr:1} 304 Isal_sept_2 = {arr:1} 305 ENDELSE 249 306 ; 250 307 jpt = 1 … … 282 339 cnt = cnt+1 & blabla = 'Erp salinity damping term' 283 340 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 284 341 ; 285 342 cnt = cnt+1 & blabla = 'Evaporation  Precipitation  Runoff term' 286 343 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 287 344 ; 288 345 cnt = cnt+1 & blabla = 'Net heat flux' 289 346 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 290 347 ; 291 348 cnt = cnt+1 & blabla = 'Meridionnal Heat Transport' 292 349 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 293 350 ; 294 351 cnt = cnt+1 & blabla = 'Global Barotropic stream Function' 295 352 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 296 353 ; 297 354 cnt = cnt+1 & blabla = 'mean Temperature diff with New Reynolds' 298 355 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 299 356 ; 300 357 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus' 301 358 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 302 359 ; 303 360 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus' 304 361 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex 305 362 ; 306 363 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2' 307 364 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 308 365 ; 309 366 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 310 367 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex 311 368 ; 312 369 cnt = cnt+1 & blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 313 370 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 314 371 ; 315 372 cnt = cnt+1 & blabla = 'mean Temperature diff with Levitus at z=100 meters' 316 373 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 317 374 ; 318 375 cnt = cnt+1 & blabla = 'mean Salinity diff with Levitus at z=100 meters' 319 376 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 320 321 ;;SF need monthly climatology to plot difference of MAX of MLD 322 cnt = cnt+1 & blabla = 'Mixed layer depth & difference with DeBoyer' 323 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex 324 377 ; 378 cnt = cnt+1 & blabla = 'Mixed layer depth' 379 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 380 ; 325 381 cnt = cnt+1 & blabla = 'Mixed layer depth differences' 326 382 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex … … 339 395 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Global' 340 396 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 341 397 ; 342 398 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 343 399 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 344 400 ; 345 401 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Indian' 346 402 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 347 403 ; 348 404 cnt = cnt+1 & blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 349 405 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 350 406 ; 351 407 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Global' 352 408 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 353 409 ; 354 410 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 355 411 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 356 412 ; 357 413 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Indian' 358 414 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex 359 415 ; 360 416 cnt = cnt+1 & blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 361 417 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 362 418 ; 363 419 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 364 420 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 365 421 ; 366 422 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 367 423 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 368 424 ; 369 425 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 370 426 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 371 427 ; 372 428 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 373 429 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 374 430 ; 375 431 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 376 432 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 377 433 ; 378 434 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' 379 435 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 380 436 ; 437 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 438 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 439 ; 440 cnt = cnt+1 & blabla = 'Arctic Ice Age: SEPT' 441 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 442 ; 443 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 444 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 445 ; 446 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' 447 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 448 ; 381 449 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: MARCH' 382 450 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 383 451 ; 384 452 cnt = cnt+1 & blabla = 'Arctic Ice Fraction: SEPT' 385 453 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 386 454 ; 387 455 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 388 456 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 389 457 ; 390 458 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 391 459 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 392 393 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 394 ;;SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 395 ;;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 396 ;; SF 397 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 398 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 399 ;; SF 400 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 401 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 402 ;; SF 403 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 404 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 405 ; 406 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 407 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 408 409 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 410 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 411 460 ; 461 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: MARCH' 462 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 463 ; 464 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 465 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 466 ; 467 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 468 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 469 ; 470 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 471 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 472 ; 473 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: MARCH' 474 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 475 ; 476 cnt = cnt+1 & blabla = 'Arctic Ice Salinity: SEPT' 477 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 478 ; 479 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: FEBRUARY' 480 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 481 ; 482 cnt = cnt+1 & blabla = 'Antarctic Ice Salinity: SEPT' 483 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 484 ; 412 485 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)' 413 486 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 414 487 ; 415 488 cnt = cnt+1 & blabla = 'Meridional stream Function: Atlantic' 416 489 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 417 490 ; 418 491 cnt = cnt+1 & blabla = 'Meridional stream Function: Indian' 419 492 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 420 493 ; 421 494 cnt = cnt+1 & blabla = 'Meridional stream Function: IndoPacific' 422 495 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 423 496 ; 424 497 cnt = cnt+1 & blabla = 'Equatorial Temperature' 425 498 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 426 499 ; 427 500 cnt = cnt+1 & blabla = 'Equatorial Salinity' 428 501 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 429 502 ; 430 503 cnt = cnt+1 & blabla = 'Equatorial zonal velocity' 431 504 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 432 505 ; 433 506 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=700' 434 507 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 435 508 ; 436 509 cnt = cnt+1 & blabla = 'Mediterranean salt tongue at depth=1000' 437 510 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 438 511 ; 439 512 cnt = cnt+1 & blabla = 'Mediterranean water at lat=40°N' 440 513 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 441 514 ; 442 515 cnt = cnt+1 & blabla = 'Mediterranean water at lat=38°N' 443 516 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 444 517 ; 445 518 cnt = cnt+1 & blabla = 'Vertical Global mean T & S' 446 519 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 447 520 ; 448 521 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 449 522
Note: See TracChangeset
for help on using the changeset viewer.