[RegCNET] Time step problem or something else ?

Dilip Ganguly dganguly at prl.res.in
Thu Jan 4 11:08:13 CET 2007


Many thanks Jagadish for your willingness to help.

As you are telling that changing domain boundary, convective scheme or 
increasing resolution may help. But my problem is I have already completed 
a long simulation and all this problem is occurring when there is just a 
small time remaing compared to the total time of simulation.

Can anyone suggest a solution to my problem ?

Thanks

Dilip.

On Thu, 4 Jan 2007, jagadish karmacharya wrote:

> Dear Dilip,
> 
> You have tried all the way down to dt=1 sec but as you have seen decreasing time step doesn't always solve the problem. I generally decrease  dt upto 1.5*resolution in km (instead of dt=3*res) to check if if overcome the problem.
> When the problem persists with decreased timestep as well it is difficult to figure out a  solution.  
> 
> Changing domain boundary, convective scheme, increasing resolution may help but if one has already completed long simulation then it become a problem.  One possibility is that your ICBC or driving data file for that month may be corrupt  but its less likely.
> 
> Can anyone give some suggestion?
> 
> Regards
> Jagadish
> 
> 
> Dilip Ganguly <dganguly at prl.res.in> wrote: Dear all,
> 
> I am facing problem in running RegCM3. The program hangs and does not seem 
> to work after a particular time and I am getting messages on the screen 
> like the following.
> 
>  BATS variables written at   2004120121   180.000000000000
>      at day =  396.8762, ktau =     210200 :  1st, 2nd time deriv of ps =  
> 0.10915E-02 0.49528E-05,  no. of points w/convection =   81
>      at day =  396.8849, ktau =     210250 :  1st, 2nd time deriv of ps =  
> 0.97549E-03 0.41641E-05,  no. of points w/convection =   71
>      at day =  396.8936, ktau =     210300 :  1st, 2nd time deriv of ps =  
> NaN        NaN        ,  no. of points w/convection =    0
> 
> -------------------------------------------------------------------------
> 
> It do not know if it is due to improper choice of time steps. I have 
> restarted the program several times with lower and lower values of time 
> step dt but every time the program hangs around a certain period of time 
> and not able to cross this period. I have tried with numerous combinations 
> of gradually decreasing time steps. Some are given below.
> 
> radfrq= 30 30 30 30 15 15 15
> abemh = 18 18 18 18 9 9 9 
> abatm = 300 120 90 45 30 9 3
> dt    = 100 60 30 15 10 3 1
> 
> But none of them worked. Everytime it shows the same problem possibly due 
> to convective instability as told by Jagadish. I am using a 90 km 
> horizontal resolution in my simulation.
> 
> Can anyone suggest a way to solve this problem ?
> 
> Thanking you in advance.
> 
> Dilip. 
>   
>      ---------------------------------------------------------------------
>         DILIP GANGULY                               
>         Postdoctoral Fellow                     
>         Space And Atmospheric Sciences Division      
>         Physical Research Laboratory                
>         Ahmedabad 380009
>  Gujarat, India
>         Phone: +91 79 26314558 (Laboratory)    
>                +91 79 26300505 (Hostel)
>  Fax:   +91 79 26314659
>  Alternate Email: ganguly.dilip at gmail.com          
>             Visit me at   http://www.prl.res.in/~dganguly
>      ---------------------------------------------------------------------
> 
> 
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
> 
> 
>  __________________________________________________
> Do You Yahoo!?
> Tired of spam?  Yahoo! Mail has the best spam protection around 
> http://mail.yahoo.com 

Dilip. 
  
     ---------------------------------------------------------------------
        DILIP GANGULY                               
        Postdoctoral Fellow                     
        Space And Atmospheric Sciences Division      
        Physical Research Laboratory                
        Ahmedabad 380009
	Gujarat, India
        Phone: +91 79 26314558 (Laboratory)    
               +91 79 26300505 (Hostel)
	Fax:   +91 79 26314659
	Alternate Email: ganguly.dilip at gmail.com          
            Visit me at   http://www.prl.res.in/~dganguly
     ---------------------------------------------------------------------





More information about the RegCNET mailing list