[MITgcm-support] Script for computing CFL conditions

Habib Toye Mahamadou Kele habib.toyemahamadoukele at kaust.edu.sa
Thu Dec 3 08:37:14 EST 2015


Thanks Ryan for your reply,

Is 'advcfl_uvel_max' the same as the parameter S_a in equation (3.9) at
http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html ?



I would like a script which will compute conditions like the ones given in:

http://mitgcm.org/public/r2_manual/latest/online_documents/node86.html#SECTION003211600000000000000
(equations (2.217) and (2.218))


http://mitgcm.org/public/r2_manual/latest/online_documents/node131.html#SECTION004122100000000000000
(from (3.43) to (3.48))

http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html
(from (3.24) to (3.29))

http://mitgcm.org/public/r2_manual/latest/online_documents/node120.html
(from (3.6) to (3.9))


The result would look like:

(for the values at
http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html )


NUMERICAL STABILITY PARAMETERS COMPUTATION FROM STDOUT.0000:

A_h = 400 ms^(-1)
delta_t = 1200 s
Delta_x = ...


Condition Name  |       Condition     |    User          |  Is user
condition fulfilled?  |

numerical stability for inertial oscillations         |    S_l < 0.3    (S_l
= 4* A_h * delta_t / Delta_x^2)                                  |   0.012 <
0.3   |                     Yes                 |

stability parameter to the horizontal Laplacian friction         |  S_i <
0.5     (S_i = f^2 * delta_t^2)  |   0.0144 < 0.5       |
 Yes                  |

...




Regards,
Habib

On Tue, Dec 1, 2015 at 5:18 PM, Ryan Abernathey <ryan.abernathey at gmail.com>
wrote:

> Hi Habib,
>
> MITgcm automatically computes the CFL for you as part of the monitor
> statistics.
>
> $ grep advcfl STDOUT.0000
> (PID.TID 0000.0001) %MON advcfl_uvel_max              =
> 1.1796518160724E-01
> (PID.TID 0000.0001) %MON advcfl_vvel_max              =
> 7.8377125552173E-02
> (PID.TID 0000.0001) %MON advcfl_wvel_max              =
> 3.8263958212481E+00
> (PID.TID 0000.0001) %MON advcfl_W_hf_max              =
> 1.8899738930714E+00
> ...
>
> Cheers,
> Ryan
>
> On Tue, Dec 1, 2015 at 8:47 AM, Habib Toye Mahamadou Kele <
> habib.toyemahamadoukele at kaust.edu.sa> wrote:
>
>> Greetings,
>>
>> Is there any script or routine for computing a mitgcmuv run CFL
>> conditions for example from the "output.txt" file obtained by "./mitgcmuv >
>> output.txt" ?
>> If not, could you please provide us with one which output would look like
>> the following table?
>>
>> Condition Name  |       Condition                               |    User
>>          |  Is user condition fulfilled?  |
>>
>> condition_1         |    C1 <= 1.5                                  |
>> 5.4 <= 1.5  |                     No                   |
>>
>> condition_2         |  var_2_1*var_2_2 / var_2_3 < 1  |   0.7 < 1
>>  |                    Yes                  |
>>
>> ...
>>
>>
>> Thanks and Regards,
>> Habib TOYE
>>
>> ------------------------------
>> This message and its contents, including attachments are intended solely
>> for the original recipient. If you are not the intended recipient or have
>> received this message in error, please notify me immediately and delete
>> this message from your computer system. Any unauthorized use or
>> distribution is prohibited. Please consider the environment before printing
>> this email.
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>

-- 

------------------------------
This message and its contents, including attachments are intended solely 
for the original recipient. If you are not the intended recipient or have 
received this message in error, please notify me immediately and delete 
this message from your computer system. Any unauthorized use or 
distribution is prohibited. Please consider the environment before printing 
this email.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20151203/0fde1ef2/attachment.htm>


More information about the MITgcm-support mailing list