Hi Gustavo, Please the see attached plots from four runs I did for Katrina (start date: 2005082700) with four different LC/Ring initialization files Rich provided. I labeled the runs as G1, G2, G3 and G4 in the plots. The plots are intensity and track plots, SST and surface currents, temperature and currents the beginning of the forecast for G1, G2, G3 and G4. Let me know if I you have any questions about these runs or if you want me to do more experiments. Cheers, Biju Richard M. Yablonsky wrote: > Hello Gustavo (and Biju and Isaac), > > I am not sure if Biju has had a chance yet to run the coupled model based > on the LC/Ring initialization files I provided. I have again cc'ed Biju > and Isaac so they can follow up with you as appropriate. I will be out of > town and/or committed to other activities until June 1st, but based on my > previous email below, I believe Biju now has all of the necessary > information to provide you with the forecasts when he gets a chance to do > so. > > Take care, > Rich > > On Wed, April 22, 2009 7:50 am, Gustavo.Goni@noaa.gov wrote: > >> Hello Richard, >> >> I wonder if you had any new development on this. >> >> Thank you very much, >> >> Gustavo >> >> ----- Original Message ----- >> From: "Richard M. Yablonsky" >> Date: Saturday, March 28, 2009 10:42 am >> Subject: Re: Hello, question (fwd) >> >> >>> Biju, >>> >>> As I believe we discussed last week, Gustavo Goni has asked our >>> group to >>> help him with a project to test whether having multiple altimeters >>> makes a >>> difference in our coupled GFDL model intensity forecasts for storms >>> thattraverse the LC and/or WCR/CCRs (see detailed correspondence >>> below). In >>> summary, he has provided me with OHC maps of the GoM for Katrina >>> with 1, >>> 2, 3, and 4 altimeters. Based on these maps, I have created >>> LC/ring text >>> files for use in our ocean initalization. These text files are in >>> separate directories that I tarred together in the following >>> tarfile on >>> cyclone: >>> >>> /hurricane5/richard/oceanspin/mcsd/lcringprogram_fortpc/lcringprogram2_fortpc/lcring4goni.tar >>> >>> When you get a chance (I know you are busy with other high-priority >>> tasks), please use the files to initialize four separate coupled model >>> Katrina runs. Track, of course, will be important for intensity, so >>> please choose an initialization date/time for which we have a >>> reasonablygood track forecast in the GoM up to landfall. >>> >>> I have cc'ed Isaac and Gustavo on this email so they are kept in >>> the loop. >>> >>> Thanks, >>> Rich >>> >>> ****************************************** >>> * Richard M. Yablonsky * >>> * Ph.D. Student - Physical Oceanography * >>> * University of Rhode Island * >>> * http://www.gso.uri.edu/~richard * >>> * URI Work #: (401) 874-6589 * >>> * Email: ryablonsky@gso.uri.edu * >>> ****************************************** >>> >>> On Wed, March 25, 2009 7:10 pm, Gustavo.Goni wrote: >>> >>>> Hello Richard, >>>> >>>> I placed the new figures with a 0.5x0.5 degree grid at >>>> the same ftp location. >>>> >>>> The data used correspond to Aug 15-25, 2005, for sea >>>> height, and Aug 20 for SST. >>>> >>>> I hope this helps. >>>> >>>> Best, >>>> >>>> Gustavo >>>> >>>> >>>> >>>> On Tue, 24 Mar 2009, Richard M. Yablonsky wrote: >>>> >>>> >>>>> Date: Tue, 24 Mar 2009 21:22:34 -0400 (EDT) >>>>> From: Richard M. Yablonsky >>>>> To: Gustavo.Goni >>>>> Cc: Isaac Ginis >>>>> Subject: Re: Hello, question (fwd) >>>>> >>>>> Gustavo, >>>>> >>>>> These versions are much better. If you happen to be able to >>>>> >>> easily make >>> >>>>> versions with a grid of 1/2 degree, that would be even better. >>>>> >>> Also,>> please let me know which date these fields are valid. >>> >>>>> Thanks, >>>>> Rich >>>>> >>>>> On Tue, March 24, 2009 7:28 pm, Gustavo.Goni wrote: >>>>> >>>>>> Hello Richard, >>>>>> >>>>>> I placed the new versions in the same ftp site. The grid >>>>>> is now of 1degree and the altimeter groundtracks are gone. >>>>>> >>>>>> The dif files are more or less to see that although some >>>>>> fields look alike, they are still different from the one >>>>>> that used only one altimeter. >>>>>> >>>>>> Thank you very much, >>>>>> >>>>>> Gustavo >>>>>> >>>>>> >>>>>> On Tue, 24 Mar 2009, Richard M. Yablonsky wrote: >>>>>> >>>>>> >>>>>>> Date: Tue, 24 Mar 2009 18:40:00 -0400 (EDT) >>>>>>> From: Richard M. Yablonsky >>>>>>> To: Gustavo.Goni >>>>>>> Cc: Isaac Ginis >>>>>>> Subject: Re: Hello, question (fwd) >>>>>>> >>>>>>> Gustavo, >>>>>>> >>>>>>> The TCHP fields themselves work fine, but a 5-degree increment >>>>>>> >>> between>>>> lat/lon grid lines makes it very difficult to pin down >>> the lat/lon >>> >>>>>>> positions of the features. If you could add more lat/lon grid >>>>>>> >>> lines>>>> so >>> >>>>>>> that the increment is perhaps every 1/2 degree (or at the very >>>>>>> >>> least>>>> every >>> >>>>>>> 1 degree), that would help greatly. Also, the satellite >>>>>>> >>> tracks are >>> >>>>>>> unnecessary, and removing them would help to unclutter the plots. >>>>>>> >>>>>>> Thanks, >>>>>>> Rich >>>>>>> >>>>>>> On Tue, March 24, 2009 5:48 pm, Gustavo.Goni wrote: >>>>>>> >>>>>>>> Hello Richard, >>>>>>>> >>>>>>>> I wonder if you could take a look at the fields that I placed >>>>>>>> at: ftp://ftp.aoml.noaa.gov/pub/phod/goni/RICHARD/ >>>>>>>> >>>>>>>> The *dif* files are the differences between the fields, with >>>>>>>> respect to the one created using only data from one >>>>>>>> altimeter (tchp1.ps). >>>>>>>> >>>>>>>> Please let me know if this is what you more or less need. >>>>>>>> >>>>>>>> Thank you very much for your time and help. >>>>>>>> >>>>>>>> Gustavo >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> On Wed, 18 Mar 2009, Richard M. Yablonsky wrote: >>>>>>>> >>>>>>>> >>>>>>>>> Date: Wed, 18 Mar 2009 10:35:30 -0400 (EDT) >>>>>>>>> From: Richard M. Yablonsky >>>>>>>>> To: Gustavo.Goni >>>>>>>>> Cc: iginis@gso.uri.edu >>>>>>>>> Subject: Re: Hello, question (fwd) >>>>>>>>> >>>>>>>>> Gustavo, >>>>>>>>> >>>>>>>>> No problem... I am happy to help. Projection is not >>>>>>>>> >>> important, as >>> >>>>>>>>> long >>>>>>>>> as >>>>>>>>> the correct TCHP values are located at the correct lat/lon >>>>>>>>> positions. >>>>>>>>> Also, the more lat/lon grid lines you can draw on the map, the >>>>>>>>> better >>>>>>>>> (as >>>>>>>>> long as the TCHP color and/or contours are still visible). >>>>>>>>> >>>>>>>>> Basically, our methodology requires us to look at the map for >>>>>>>>> regions >>>>>>>>> of >>>>>>>>> strong gradient (TCHP gradient, in this case), thereby >>>>>>>>> >>> signifying>>>>>> the >>> >>>>>>>>> perimeter of the features (LC, WCRs, CCRs), and then we >>>>>>>>> >>> enter the >>> >>>>>>>>> lat/lon >>>>>>>>> position of points along the perimeter (to the nearest 0.1 >>>>>>>>> >>> degree)>>>>>> into >>> >>>>>>>>> our initialization procedure. So, for example, if the map >>>>>>>>> >>> shows a >>> >>>>>>>>> strong >>>>>>>>> TCHP gradient at the northern end of the LC at (-86.4N, >>>>>>>>> >>> 25.5N), I >>> >>>>>>>>> might >>>>>>>>> use this lat/lon pair to represent the northern end of the >>>>>>>>> >>> LC in our >>> >>>>>>>>> initialization procedure. >>>>>>>>> >>>>>>>>> Hope this helps! >>>>>>>>> >>>>>>>>> Rich >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> On Wed, March 18, 2009 10:01 am, Gustavo.Goni wrote: >>>>>>>>> >>>>>>>>>> Hello Richard, >>>>>>>>>> >>>>>>>>>> Thank you very much for your email and prompt reply. >>>>>>>>>> >>>>>>>>>> Sure, I can create these maps. I wonder if you need any >>>>>>>>>> >>> specific>>>>>>> type of projection (Mercator ?) >>> >>>>>>>>>> Good luck with the preparation of your dissertation and >>>>>>>>>> >>> defense ! >>> >>>>>>>>>> Gustavo >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> On Wed, 18 Mar 2009, Richard M. Yablonsky wrote: >>>>>>>>>> >>>>>>>>>> >>>>>>>>>>> Date: Wed, 18 Mar 2009 09:37:38 -0400 (EDT) >>>>>>>>>>> From: Richard M. Yablonsky >>>>>>>>>>> To: Gustavo.Goni >>>>>>>>>>> Cc: iginis@gso.uri.edu >>>>>>>>>>> Subject: Re: Hello, question (fwd) >>>>>>>>>>> >>>>>>>>>>> Gustavo, >>>>>>>>>>> >>>>>>>>>>> I would be happy to help out with this mini-project. >>>>>>>>>>> >>> However, it >>> >>>>>>>>>>> would >>>>>>>>>>> be >>>>>>>>>>> a bit time-consuming for me to develop code to process and >>>>>>>>>>> >>> plot>>>>>>>> the >>> >>>>>>>>>>> tchp >>>>>>>>>>> data files you provided on the ftp site. If you already >>>>>>>>>>> >>> have a >>> >>>>>>>>>>> program >>>>>>>>>>> that can do it easily, could you please create tchp maps >>>>>>>>>>> >>> based on >>> >>>>>>>>>>> the >>>>>>>>>>> data >>>>>>>>>>> you provided? Ideally, it would be nice if the lat/lon >>>>>>>>>>> >>> range of >>> >>>>>>>>>>> the >>>>>>>>>>> maps >>>>>>>>>>> covered only the Gulf of Mexico... perhaps something like >>>>>>>>>>> >>> 98W to >>> >>>>>>>>>>> 81W >>>>>>>>>>> and >>>>>>>>>>> 18N to 31N. >>>>>>>>>>> >>>>>>>>>>> Thanks, >>>>>>>>>>> Rich >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> On Tue, March 17, 2009 12:50 pm, Isaac Ginis wrote: >>>>>>>>>>> >>>>>>>>>>>> Gustavo, >>>>>>>>>>>> >>>>>>>>>>>> Thanks. I will ask Rich Yablonsky to initialize the ocean >>>>>>>>>>>> >>> model>>>>>>>>> based >>> >>>>>>>>>>>> on >>>>>>>>>>>> your fields and rerun Katrina forecasts. He is finishing >>>>>>>>>>>> >>> up his >>> >>>>>>>>>>>> dissertation this month, so this would probably need to wait >>>>>>>>>>>> until >>>>>>>>>>>> he >>>>>>>>>>>> is >>>>>>>>>>>> done. >>>>>>>>>>>> >>>>>>>>>>>> Isaac >>>>>>>>>>>> >>>>>>>>>>>> Gustavo.Goni wrote: >>>>>>>>>>>> >>>>>>>>>>>>> Hello Isaac, >>>>>>>>>>>>> >>>>>>>>>>>>> I just replaced the fields that I placed last night int >>>>>>>>>>>>> the ftp because I had forgotten to multiply their values by >>>>>>>>>>>>> rho*cp. >>>>>>>>>>>>> >>>>>>>>>>>>> Gustavo >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> ---------- Forwarded message ---------- >>>>>>>>>>>>> Date: Mon, 16 Mar 2009 19:49:14 -0400 (EDT) >>>>>>>>>>>>> From: Gustavo.Goni >>>>>>>>>>>>> To: Isaac Ginis >>>>>>>>>>>>> Subject: Re: Hello, question >>>>>>>>>>>>> >>>>>>>>>>>>> Hello Isaac, >>>>>>>>>>>>> >>>>>>>>>>>>> Thank you very much for your prompt response. I really >>>>>>>>>>>>> >>> like>>>>>>>>>> the manuscript, as I was not very familiar on what >>> the process >>> >>>>>>>>>>>>> was to incorporate the heat content maps. >>>>>>>>>>>>> >>>>>>>>>>>>> I placed in my anonymous ftp: >>>>>>>>>>>>> ftp://ftp.aoml.noaa.gov/pub/phod/goni/ISAAC/ >>>>>>>>>>>>> four files of upper ocean heat content. I wonder if you >>>>>>>>>>>>> could run your analysis in the Gulf for Katrina for each >>>>>>>>>>>>> of them ? They were computed using data from 1,2,3, >>>>>>>>>>>>> >>> and 4 >>> >>>>>>>>>>>>> altimeters. >>>>>>>>>>>>> >>>>>>>>>>>>> I am very curious to see if intensity (I believe you >>>>>>>>>>>>> >>> compute>>>>>>>>>> minimum pressure) gets closer to obs with more >>> than one >>> >>>>>>>>>>>>> altimeter. >>>>>>>>>>>>> The fields correspond to August 15-25, 2005. >>>>>>>>>>>>> >>>>>>>>>>>>> Thank you very much, >>>>>>>>>>>>> >>>>>>>>>>>>> Gustavo >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> On Fri, 13 Mar 2009, Isaac Ginis wrote: >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>>> Date: Fri, 13 Mar 2009 13:00:10 -0400 >>>>>>>>>>>>>> From: Isaac Ginis >>>>>>>>>>>>>> To: Gustavo.Goni >>>>>>>>>>>>>> Subject: Re: Hello, question >>>>>>>>>>>>>> >>>>>>>>>>>>>> Gustavo, >>>>>>>>>>>>>> >>>>>>>>>>>>>> I have attached a paper that describes the methodology >>>>>>>>>>>>>> >>> we use >>> >>>>>>>>>>>>>> to >>>>>>>>>>>>>> initialize >>>>>>>>>>>>>> the LC and eddies in the GoM. If you can provide OHC >>>>>>>>>>>>>> >>> maps in >>> >>>>>>>>>>>>>> the >>>>>>>>>>>>>> region, >>>>>>>>>>>>>> shown in Fig. 2, that would be sufficient. The grid >>>>>>>>>>>>>> >>> spacing in >>> >>>>>>>>>>>>>> our >>>>>>>>>>>>>> model is >>>>>>>>>>>>>> 9 km, but we can deal with any grid resolution you can >>>>>>>>>>>>>> provide. >>>>>>>>>>>>>> We >>>>>>>>>>>>>> are >>>>>>>>>>>>>> going to start the forecasts on Aug. 26 00Z, 2005, >>>>>>>>>>>>>> >>> shortly>>>>>>>>>>> after >>> >>>>>>>>>>>>>> Katrina >>>>>>>>>>>>>> left Florida. So we will need the OHC maps at or >>>>>>>>>>>>>> >>> before that >>> >>>>>>>>>>>>>> time. >>>>>>>>>>>>>> >>>>>>>>>>>>>> Regards, >>>>>>>>>>>>>> >>>>>>>>>>>>>> Isaac >>>>>>>>>>>>>> >>>>>>>>>>>>>> Gustavo.Goni wrote: >>>>>>>>>>>>>> >>>>>>>>>>>>>>> Hello Isaac, >>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Thank you very much. We can produce the heat content >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>> fields using 1-2-3 altimeters. When you increase the >>>>>>>>>>>>>>> number of altimeters the mesoscale features show >>>>>>>>>>>>>>> >>> better.>>>>>>>>>>>> Some mesoscale features do not even show up >>> when you >>> >>>>>>>>>>>>>>> have one altimeter only. >>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Could you please let me know the region, the >>>>>>>>>>>>>>>> >>> horizontal>>>>>>>>>>>> grid scale and the temporal resolution >>> (daily) from what >>> >>>>>>>>>>>>>>> day to what day, you need to make these runs ? >>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Thank you very much again, >>>>>>>>>>>>>>>> Gustavo >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> On Wed, 11 Mar 2009, Isaac Ginis wrote: >>>>>>>>>>>>>>>>> Date: Wed, 11 Mar 2009 10:27:59 -0400 >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> From: Isaac Ginis >>>>>>>>>>>>>>>> To: Gustavo.Goni@noaa.gov >>>>>>>>>>>>>>>> Subject: Re: Hello, question >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> Gustavo, >>>>>>>>>>>>>>>>>> It should not be hard for us to rerun the Katrina >>>>>>>>>>>>>>>>>> >>> case.>>>>>>>>>>>>>>> In >>> >>>>>>>>>>>>>> the GFDL > > coupled system, we can initialize the ocean >>>>>>>>>>>>>> model >>>>>>>>>>>>>> based on different > > ocean analysis. What we will >>>>>>>>>>>>>> >>> need from >>> >>>>>>>>>>>>>> you >>>>>>>>>>>>>> is the OHC maps that are > > based on one and two >>>>>>>>>>>>>> >>> altimeter.>>>>>>>>>>>>>>> Of course, one relevant question would >>> be: are there any >>> >>>>>>>>>>>>>> material > > differences in the position and >>>>>>>>>>>>>> >>> structure of the >>> >>>>>>>>>>>>>> GS >>>>>>>>>>>>>> and the warm-core > > eddy between the two maps? >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> Regards, >>>>>>>>>>>>>>>>>> Isaac >>>>>>>>>>>>>>>>>> Gustavo.Goni@noaa.gov wrote: >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> Hello Isaac, >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>> I wonder if I could ask you a question, maybe >>>>>>>>>>>>>>>>>>>> >>> a favor. >>> >>>>>>>>>>>>>>>>> You constructed the figure that has the GFDL model >>>>>>>>>>>>>>>>> results >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>> with the >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> Minimum sea level pressure in function of time for >>>>>>>>>>>>>>>>> Katrina, >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>> with > > > the plot >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> containing the actual observations and the output >>>>>>>>>>>>>>>>> >>> from the >>> >>>>>>>>>>>>>> GFDL model > > > with >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> and without heat content. >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>> I was wondering if you could make this run >>>>>>>>>>>>>>>>>>>> >>> again if I >>> >>>>>>>>>>>>>> provide to you > > > the >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> ocean heat content computed with >>>>>>>>>>>>>>>>> one altimetry and with two altimeters. The goal >>>>>>>>>>>>>>>>> >>> is to >>> >>>>>>>>>>>>>>>>> see >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>> if the > > > pressure >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> obtained with >>>>>>>>>>>>>>>>> the ocean heat content input using two altimeters is >>>>>>>>>>>>>>>>> better >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>> than > > > the one >>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> obtained using >>>>>>>>>>>>>>>>> the OHC with only one altimeter. >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>> Thank you very much, >>>>>>>>>>>>>>>>>>>> Gustavo >>>>>>>>>>>>>>>>>>>>