[MITgcm-support] EXF and Calendar Package
Martin Losch
Martin.Losch at awi.de
Tue Aug 5 02:27:09 EDT 2014
Hi Andres,
dz(k+1)/dz(k) < sqrt(2) is a rule of thumb (that you should follow). The finite volume discretization requires some averaging to grid cell interfaces and currently (and I don’t believe that this will change) this averaging requires that the grid spacing varies smoothly to avoid large discretrization errors. The same is true for the horizontal grid spacing, but in the horizontal very irregular grid appear unnatural, right?
Can you gzip/tar your code and input directories and put them somewhere to download and look at (please don’t send large files to this list)?
For problems like yours, it makes sense to start from a working example and change it file by file to the actual configuration that you want to run, so that you can identify the problems individually.
Martin
On Aug 5, 2014, at 12:46 AM, Andres Fernandez <andres.fernandez at nyu.edu> wrote:
> Data file, STDOUT.0000 and the .e file that is generated are attached...
> the STDERR.0000 is blank, so i did not attach it.
>
> @Nicolas, I was not aware of this limitation of irregular spacing. Any comments from other users?
>
>
> Thank you all!
>
>
> Andres Fernandez
> Undergraduate Chemical Engineering
> New York University Abu Dhabi - 2015
> Mobile (USA): +1 646 717 2231
>
> *Think before printing. The environment is in our hands
>
>
>
> On Mon, Aug 4, 2014 at 5:23 PM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
> Hi Andres,
>
> It would be helpful for you to share your data file…plus the last few lines from STDOUT.0000 and the output from STDERR.0000.
>
> Dustin
>
> On Aug 4, 2014, at 7:15 PM, Andres Fernandez <andres.fernandez at nyu.edu> wrote:
>
>> Dear Dustin,
>> I realized I had a comma missing in the list when copying you the delR declaration. I have an irregular spacing declared as follows: So, thank you.
>>
>> delR= 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 5.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 10.0,
>> 50.0,
>> 50.0,
>> 50.0,
>> 50.0,
>> 50.0,
>> &
>>
>>
>>
>>
>> ** I do still get an error though!
>>
>> forrtl: severe (24): end-of-file during read, unit 11, file .../fortzlTn72 line 8, position 0
>> "" ... / fort5xkez2
>> "" .../ fort2YWOn2fort5xkez2
>> and a bunch more files..
>>
>> Unsure if you could assist with this error?
>>
>> @Jean- Michel, I have ran a couple models before and the struggle I have right now is with external forcing, I want to add precipitation and evaporation and I am not succeeding. I have not found any comprehensive information on the manual.
>>
>>
>>
>> Andres Fernandez
>> Undergraduate Chemical Engineering
>> New York University Abu Dhabi - 2015
>> Mobile (USA): +1 646 717 2231
>>
>> *Think before printing. The environment is in our hands
>>
>>
>>
>> On Mon, Aug 4, 2014 at 3:58 PM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
>> Hi Andres,
>>
>> The vertical grid spacing (delZ, delP, or delR) is set in your data file. You can specific with a *.bin file, i.e delRfile=‘dz. bin’
>>
>> or without a bin file...delR=100*10
>>
>> Cheers,
>> Dustin
>>
>>
>>
>>
>> On Aug 4, 2014, at 5:43 PM, Andres Fernandez <andres.fernandez at nyu.edu> wrote:
>>
>>> Dear support,
>>>
>>> I get the (PID.TID 0001.0001) *** ERROR *** S/R INI_PARMS: No value for delZ/delP/delR at k = 16 ....25
>>> I am unsure where to set these initial conditions... or is this when I build it?
>>>
>>> Did I build it incorrectly?
>>>
>>> Thank you
>>>
>>> Andres Fernandez
>>> Undergraduate Chemical Engineering
>>> New York University Abu Dhabi - 2015
>>> Mobile (USA): +1 646 717 2231
>>>
>>> *Think before printing. The environment is in our hands
>>>
>>>
>>>
>>> On Thu, Jul 31, 2014 at 8:34 PM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
>>> Hi Andres,
>>>
>>> You can debug this further by looking at your STDERR.0000 file for more specific error information... also take a look at STDOUT.0000 to see where your model crashed. You probably have a syntax error in one of your data files.
>>>
>>> Cheers,
>>> Dustin
>>>
>>>
>>> On Jul 31, 2014, at 9:21 PM, Andres Fernandez <andres.fernandez at nyu.edu> wrote:
>>>
>>>> So I recompiled the executable; and placed it in my input folder;
>>>>
>>>> I have the following folders:
>>>>
>>>> data
>>>> data.cal
>>>> data.exf
>>>> data.obcs
>>>> data.pkg
>>>> data.rbcs
>>>> eedata
>>>> eedata.mth
>>>> executable
>>>> (some binary files)
>>>>
>>>> And I get this error when I run:
>>>> ABNORMAL END: S/R INI_PARMS
>>>>
>>>> Thank you for your assistance.
>>>>
>>>>
>>>> Best,
>>>> Andres
>>>>
>>>> Andres Fernandez
>>>> Mobile (USA): +1 646 717 2231
>>>>
>>>> *Think before printing. The environment is in our hands
>>>>
>>>>
>>>>
>>>> On Thu, Jul 31, 2014 at 3:46 PM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
>>>> Hi Andres,
>>>>
>>>> also, under code/packages.conf add the line:
>>>>
>>>> cal
>>>>
>>>> and recompile. Then you should be fine with useEXF=.TRUE., in your data.pkg file.
>>>>
>>>> Cheers,
>>>>
>>>>
>>>> Dustin
>>>>
>>>> On Jul 31, 2014, at 3:53 PM, Andres Fernandez <andres.fernandez at nyu.edu> wrote:
>>>>
>>>>> Thank you!
>>>>> So building a code with the following packages should allow me to use external forcing correctly?
>>>>>
>>>>> # Packages
>>>>> &PACKAGES
>>>>> useOBCS=.TRUE.,
>>>>> usePtracers=.FALSE.,
>>>>> useRBCS=.FALSE.,
>>>>> useEXF=.TRUE.,
>>>>> #useLayers=.TRUE.,
>>>>> &
>>>>>
>>>>> Best,
>>>>> Andres
>>>>>
>>>>> Andres Fernandez
>>>>> Undergraduate Chemical Engineering
>>>>> New York University Abu Dhabi - 2015
>>>>> Mobile (USA): +1 646 717 2231
>>>>>
>>>>> *Think before printing. The environment is in our hands
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Jul 31, 2014 at 12:08 PM, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:
>>>>> Hi Andres,
>>>>>
>>>>> Since pkg/exf requires the use of calendar package (pkg/cal),
>>>>> this later (pkg/cal) is always compiled (MITgcm/pkg/pkg_depend)
>>>>> and used (in: model/src/packages_boot.F: IF (useEXF) useCAL = .TRUE.)
>>>>> anytime pkg/exf is compiled and used.
>>>>>
>>>>> Cheers,
>>>>> Jean-Michel
>>>>>
>>>>> On Thu, Jul 31, 2014 at 11:30:08AM -0400, Andres Fernandez wrote:
>>>>> > Dear MITgcm support,
>>>>> >
>>>>> > I have been trying to run the external forcing package; and I have yet to
>>>>> > succeed, but I was wondering if a reason why it is not working is because I
>>>>> > do not have the calendar package?
>>>>> >
>>>>> > Thank you,
>>>>> >
>>>>> > *Andres Fernandez*
>>>>> > Mobile (USA): +1 646 717 2231
>>>>> >
>>>>> > *Think before printing. The environment is in our hands
>>>>>
>>>>> > _______________________________________________
>>>>> > 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
>>>>>
>>>>> _______________________________________________
>>>>> 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
>>>>
>>>>
>>>> _______________________________________________
>>>> 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
>>>
>>>
>>> _______________________________________________
>>> 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
>>
>>
>> _______________________________________________
>> 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
>
>
> <STDOUT.0000.odt><DATA.odt><dotEfile.odt>_______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list