Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Tue Nov 21, 2017 12:50 am




Post new topic Reply to topic  [ 11 posts ] 

All times are UTC

Author Message
 Post subject: ecmwf2roms
PostPosted: Sun Dec 08, 2013 2:45 pm 
Offline

Joined: Sun Dec 30, 2012 2:58 pm
Posts: 77
Location: inio:Iranian National Institute for Oceanography
Dear all user
for using d_ecmwf2roms.m file all needed files exist in ECMWF_ERA Interim except shflux, swflux, sustr, svstr, and photosynthetically active radiation at surface. i have 2 question:

1- for computing swflux, evaporation is necessary while it have not defined in that m file? Why?

2- sustr and svstri are exactly just same Eastward turbulent surface stress and Eastward turbulent surface stress, respectivly whish are in ECMWF2ROMS?

Many thanks to receive guidance from someone who has experiment to use that routin.
Cheers
fereshte


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Mon Dec 09, 2013 2:16 pm 
Offline

Joined: Sun Dec 30, 2012 2:58 pm
Posts: 77
Location: inio:Iranian National Institute for Oceanography
sorry i can find them
but about evaporation i belive that it should be added to that m file,is not it?
thanks


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Mon Dec 09, 2013 9:45 pm 
Offline

Joined: Fri Nov 19, 2010 2:33 pm
Posts: 51
Location: University of Aegean
Dear fereshteh,

if you check the lines from 465 to 468 of d_ecmwf2roms.m you will see
that evaporation is used on case "swflux".

Have a pleasant evening
Giannis


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Tue Dec 10, 2013 6:21 am 
Offline

Joined: Sun Dec 30, 2012 2:58 pm
Posts: 77
Location: inio:Iranian National Institute for Oceanography
Dear ymamoutos
many thanks for your answer
i have got what is my problem :wink:
sincerely yours
fereshte


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Mon Sep 26, 2016 8:42 am 
Offline

Joined: Tue Nov 03, 2015 12:52 am
Posts: 17
Location: School of Space Science of USTC
fereshteh wrote:
Dear ymamoutos
many thanks for your answer
i have got what is my problem :wink:
sincerely yours
fereshte

hello , in the file 'd_ecmwf2roms.m',there have the lines as follow to change the 'previous accumulation'.
but I find that the variable 'field_previous' used without in initial value. How to correct it?
Code:
 step = rem(frc_time,0.5)*24;
        if step == 3
          fieldfinal = field;
        else
          fieldfinal = field - field_previous;  % At other steps subtract
        end                                     % the previous accumulation

        frc_time = frc_time - 1.5/24;           % Center forcing time on the
                                                % accumulation interval

        field_previous = field;                 % Save this accumulation
                                                % to on the next step


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Tue Sep 27, 2016 11:07 am 
Offline
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 421
Location: Rutgers University
I wrote the script assuming that 3-hour interval data were downloaded for processing. Hence the step==3 condition is true on the first pass through and field_previous gets initialized with the hour 3 data.

The comments at the start note that ...

Quote:
% The following parameters are used to extract ERA-Interim fields:
%
% Select time: 00:00:00 12:00:00
%
% Select step: 0 3 6 9 12
%

_________________
John Wilkin: IMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Thu Feb 23, 2017 5:51 am 
Offline

Joined: Thu Sep 29, 2016 4:03 pm
Posts: 2
Location: Nanjing University of Information Sci. &Technology
Code:
 if (F(n).scale < 0),
        step = rem(frc_time,0.5)*24;
        if step == 3
          fieldfinal = field;
        elseif step == 0
           
          fieldfinal = field - field_previous;  % At other steps subtract
        end                                     % the previous accumulation

        frc_time = frc_time - 1.5/24;           % Center forcing time on the
                                                % accumulation interval

        field_previous = field;                 % Save this accumulation
                                                % to on the next step
      end

i met this question. when i run this script ,it showed that
** Processing: sustr for 01-Jan-2015 **
Undefined function or variable 'field_previous'.

Error in d_ecmwf2roms (line 509)
fieldfinal = field - field_previous; % At other steps subtract

i want to know how to correct it .

thank you in advance.


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Thu Feb 23, 2017 7:36 am 
Offline
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 421
Location: Rutgers University
Did you download 3 hour interval data as discussed previously in this thread?

_________________
John Wilkin: IMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Fri Feb 24, 2017 8:36 am 
Offline

Joined: Mon Mar 03, 2014 1:57 pm
Posts: 10
Location: SKLEC, ECNU, China
Caffery wrote:
Code:
 if (F(n).scale < 0),
        step = rem(frc_time,0.5)*24;
        if step == 3
          fieldfinal = field;
        elseif step == 0
           
          fieldfinal = field - field_previous;  % At other steps subtract
        end                                     % the previous accumulation

        frc_time = frc_time - 1.5/24;           % Center forcing time on the
                                                % accumulation interval

        field_previous = field;                 % Save this accumulation
                                                % to on the next step
      end

i met this question. when i run this script ,it showed that
** Processing: sustr for 01-Jan-2015 **
Undefined function or variable 'field_previous'.

Error in d_ecmwf2roms (line 509)
fieldfinal = field - field_previous; % At other steps subtract

i want to know how to correct it .

thank you in advance.


I added another judgement to it, that is:
Code:
      if ( F(n).scale < 0 )
        step = rem(frc_time,0.5)*24;    % remainder after division
        if step == 3
          fieldfinal = field;
          frc_time = frc_time - 1.5/24; % Center forcing time on the
                                        % accumulation interval
        % we don't know what the 1st time record is (maybe 6, 9 ...)
        % then admit it as appropriate data despite some mistakes
        % introduced
        elseif MyRec(n) == 1           
            fieldfinal = field;         
            frc_time = frc_time - step/48;  % at center:   step/2 hours
                                            %              step/48 days
        else
            fieldfinal = field - field_previous;  % Other steps subtract
                                                  % the previous accumulation
            frc_time = mean([time( indx(MyRec(n) - 1)), frc_time]); % average
        end                                     

        field_previous = field;                 % Save this accumulation
                                                % to on the next step
      end


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Fri Feb 24, 2017 10:03 am 
Offline
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 421
Location: Rutgers University
That code modification is going to give you the wrong values.

For data that are integrated over time during the analysis interval, like precipitation, you must unravel the time accumulation. The only way to do this it to start the input data from the 3-hour step in ERA files. In this case, the loop never throws an error because the else case is never processed first.

_________________
John Wilkin: IMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu


Top
 Profile  
Reply with quote  
 Post subject: Re: ecmwf2roms
PostPosted: Fri Feb 24, 2017 3:55 pm 
Offline

Joined: Mon Mar 03, 2014 1:57 pm
Posts: 10
Location: SKLEC, ECNU, China
Thanks for your correction!

In my opinion, if the added judgement works, the 1st record is the only bad one and should be discarded, is that right?


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 11 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group