[MITgcm-devel] change to ctrl_set_pack_xy for default weighttype
gael forget
gforget at mit.edu
Thu Oct 23 09:27:54 EDT 2014
> To be honest I never understood what "wuinit" is. Does that set the weight array uniformly to a predefined value?
Like I said, in my opinion, it is a feature that is only
useful for early development and testing, which you
don’t really want to use in the long run.
> &CTRL_NML_GENARR
> genarr2dPrecond = 20.,
> gentim2dPrecond = .1,.1,.1,.1,
> xx_gentim2d_startdate1 = 0,
> xx_gentim2d_startdate2 = 0,
> xx_gentim2d_period = 63072000, 63072000, 63072000, 63072000,
> /
>
> and I did NOT create any input files corresponding to weights. To get the same functionality currently I must add the lines
>
> xx_gentim2d_weight = 'wunit', 'wunit', 'wunit','wunit'
> xx_genarr2d_weight = 'wunit'
>
> to data.ctrl, and include a single binary file named
>
> wunit.bin
Looks about right at first glance.
I would use a different name though, to be safer, avoid confusion
(pkg/ctrl does generate a wunit.data output…) and indicate what
is in the file. Maybe wt_one.data in analogy with the names
I use in the ocean state estimate, which provides an example @
http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/gael/verification/global_oce_llc90/input_itXX/data.ctrl?revision=1.7&view=markup
Also it would be better to add indices in lines such
as xx_gentim2d_weight (1:4)= 'wunit', 'wunit', 'wunit',’wunit'
> (i think) with the appropriate dimensions -- in this case Nx by Ny, as these are 2D controls -- and set uniformly to 1.0, in the run directory. (Note i have ALLOW_NONDIMENSIONAL_CONTROL_IO #define'd in ECCO_CPPOPTIONS.h)
It would be better to put ALLOW_NONDIMENSIONAL_CONTROL_IO
and control related compile options in CTRL_OPTIONS.h
One more thing I should mention : not right now, but sooner rather
than later, we want to remove custom versions of ctrl_map_ini_genarr.F
that are in verification experiments. This should boil down to adding a
few lines in pkg/ctrl_map_ini_genarr.F
But let’s wait until after you update your weights/pre-conditioner run
time settings, before we start looking into those two aspects.
Please give us a heads up when you get there.
Cheers,
Gael
> Thanks
> Dan
>
>
> Cheers,
> Gael
>
>
>> Thanks
>> Dan
>>
>>
>>
>> --
>>
>> Daniel Goldberg, PhD
>> Lecturer in GlaciologySchool of Geosciences, University of Edinburgh
>> Geography Building, Drummond Street, Edinburgh EH8 9XP
>>
>>
>> em: Dan.Goldberg at ed.ac.uk
>> web: http://ocean.mit.edu/~dgoldberg_______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
>
>
>
> --
>
> Daniel Goldberg, PhD
> Lecturer in Glaciology
> School of Geosciences, University of Edinburgh
> Geography Building, Drummond Street, Edinburgh EH8 9XP
>
>
> em: Dan.Goldberg at ed.ac.uk
> web: http://ocean.mit.edu/~dgoldberg
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20141023/5c820c3f/attachment.htm>
More information about the MITgcm-devel
mailing list