new SST upwelling animations

IDL_PROGRAMS/soe_animations.pro

@@ -74,21 +74,21 @@
   FOR E=0, N_ELEMENTS(EVENTS)-1 DO BEGIN
     EVENT = EVENTS[E]
     EDIR = ANIMATION_DIR + EVENT + SL & PDIR = EDIR + 'PNGS' + SL & DIR_TEST, PDIR
-    SST_SCALE = [0,30]
-    SST_ASCALE = [-5,5]
-    CHLOR_SCALE = [0.1,10]
-    CHL_ASCALE = [0.2,6]
-    PPD_SCALE = [0.0,2.2]
-    PPD_ASCALE = [0.3,3]
+    SST_SCALE = [0,30] & SST_TICKS = [0,5,10,15,20,25,30] & SST_CONTOUR = []
+    SST_ASCALE = [-5,5] & SST_ATICKS = [-5,-2.5,0,2.5,5] 
+    CHL_SCALE = [0.1,10] & CHL_TICKS = [0.1,0.3,1,3,10]   & CHL_CONTOUR = []
+    CHL_ASCALE = [0.2,6] & CHL_ATICKS = [0.2,0.6,1,3,6]
+    PPD_SCALE = [0.1,3] & PPD_TICKS = []                  & PPD_CONTOUR = []
+    PPD_ASCALE = [0.2,6] & PPD_ATICKS = [0.2,0.6,1,3,6]
     CASE EVENT OF
       '2024_UPWELLING': BEGIN
         DTR = GET_DATERANGE([20240601,20240830])
         PRODS = LIST('SST','CHLOR_A',['SST','CHLOR_A']);,'SST-ANOMALY',['SST','SST-ANOMALY'])
         MPS = ['MAB'];,'NES','NJ_COAST']
-        SST_SCALE = [10,30]
+        SST_SCALE = [12,28] & SST_TICKS = [12,16,20,24,28] & SST_CONTOUR = 15
         SST_ASCALE = [-5,5]
         CHL_SCALE = [0.1,30]
-
+        BUFFER = 1
         PERIOD = 'D'
       END  
     ENDCASE
@@ -99,13 +99,13 @@
       IF N_ELEMENTS(PRODS[R]) GT 2 THEN MESSAGE, 'ERROR: Currently this program is only set up to work with 2 products'
 
       CASE VALIDS('PRODS',APROD) OF
-        'SST':     BEGIN & SYTITLE=UNITS('TEMP')    & AYTITLE=UNITS('TEMP',/NO_UNIT)+' Anomaly ' + UNITS('TEMP',/NO_NAME) & SRNG=SST_SCALE & ARNG=SST_ASCALE & PSTATS='MEAN' & ALG=0 & AMID=0 & ASTATS='AMEAN' & AYTICKS=[] & END
-        'CHLOR_A': BEGIN & SYTITLE=UNITS('CHLOR_A') & AYTITLE=UNITS('CHLOR_A',/NO_UNIT) + ' Ratio Anomaly'                & SRNG=CHL_SCALE & ARNG=CHL_ASCALE & PSTATS='MED'  & ALG=1 & AMID=1 & ASTATS='AMEAN' & AYTICKS=[] & END
-        'PPD':     BEGIN & SYTITLE=UNITS('PPD')     & AYTITLE=UNITS('PPD',/NO_UNIT) + ' Ratio Anomaly'                    & SRNG=PPD_SCALE & ARNG=PPD_ASCALE & PSTATS='MED'  & ALG=1 & AMID=1 & ASTATS='AMEAN' & AYTICKS=[] & END
+        'SST':     BEGIN & SYTITLE=UNITS('TEMP')    & AYTITLE=UNITS('TEMP',/NO_UNIT)+' Anomaly ' + UNITS('TEMP',/NO_NAME) & SRNG=SST_SCALE & ARNG=SST_ASCALE & STICKS=SST_TICKS & ATICKS=SST_ATICKS & SCONTOUR=SST_CONTOUR & PSTATS='MEAN' & ALG=0 & AMID=0 & ASTATS='AMEAN' & AYTICKS=[] & END
+        'CHLOR_A': BEGIN & SYTITLE=UNITS('CHLOR_A') & AYTITLE=UNITS('CHLOR_A',/NO_UNIT) + ' Ratio Anomaly'                & SRNG=CHL_SCALE & ARNG=CHL_ASCALE & STICKS=CHL_TICKS & ATICKS=CHL_ATICKS & SCONTOUR=CHL_CONTOUR & PSTATS='MED'  & ALG=1 & AMID=1 & ASTATS='AMEAN' & AYTICKS=[] & END
+        'PPD':     BEGIN & SYTITLE=UNITS('PPD')     & AYTITLE=UNITS('PPD',/NO_UNIT) + ' Ratio Anomaly'                    & SRNG=PPD_SCALE & ARNG=PPD_ASCALE & STICKS=PPD_TICKS & ATICKS=PPD_ATICKS & SCONTOUR=PPD_CONTOUR & PSTATS='MED'  & ALG=1 & AMID=1 & ASTATS='AMEAN' & AYTICKS=[] & END
       ENDCASE
       IF ATYPE EQ [] THEN APROD_SCALE = APROD + '_' + STRJOIN(NUM2STR(SRNG),'_') ELSE APROD_SCALE = APROD + '-' + STRJOIN(NUM2STR(ARNG),'_')
 
-      PROJECT_MAKE_ANIMATION, VERSTR, PRODS=APROD,PROD_SCALE=APROD_SCALE,CB_TITLE=ATITLE,MAPP=MPS, DATERANGE=DTR, /ADD_COLORBAR, /ADD_BATHY, /ADD_LONLAT,/ADD_BORDER,/ADD_DATEBAR
+      PROJECT_MAKE_ANIMATION, VERSTR, PRODS=APROD,PROD_SCALE=APROD_SCALE,CB_TITLE=ATITLE,CB_TICKVALUES=STICKS,C_LEVELS=SCONTOUR,C_COLOR=1,MAPP=MPS,DATERANGE=DTR,/ADD_COLORBAR,/ADD_BATHY,/ADD_LONLAT,/ADD_BORDER,/ADD_DATEBAR, BUFFER=BUFFER
 
 
       ASTR = VERSTR.PROD_INFO.(WHERE(TAG_NAMES(VERSTR.PROD_INFO) EQ APROD,/NULL))

IDL_PROGRAMS/soe_main.pro

@@ -145,7 +145,10 @@
       IF KEYWORD_SET(SST_PNGS)           THEN SOE_SST,VER, BUFFER=BUFFER
     ENDELSE
 
+    SST_DATASET = VERSTR.PROD_INFO.SST.DATASET
+    VERSTR.PROD_INFO.SST.DATASET = 'MUR'
     SOE_ANIMATIONS, VERSTR, EVENTS=['2024_UPWELLING'], BUFFER=1
+    VERSTR.PROD_INFO.SST.DATASET = SST_DATASET
  STOP   
     SOE_TIMESERIES_ANIMATION, VERSTR, PRODS=['CHLOR_A','PPD','SST'], BUFFER=1