[MITgcm-support] Error while creating bathymery : ini_masks_etc
Camille Mazoyer
mazoyer at univ-tln.fr
Mon Jan 30 17:18:22 EST 2017
Hi David,
Thanks for your reply.
My overlap are OLx = OLy = 3.
In fact, I wanted to have a netcdf grid file with mask values to prepare
my OBC and my init state.
So usually I first run a simu without obc and init file to obtain a
netcdf grid file with all informations (mask hfacc, ....), I then
prepare my OBC + initial file.
Maybe It's not the right way to do it?
I understand that i=-2 and j=-2 can exist, because of OLx and OLy. But I
think I have a problem with my depth or something else there because I
have a floating invalid problem there.
If you have an idea of something to watch, or to verify?
Thanks a lot,
Camille
Le 30/01/2017 17:42, David Ferreira a écrit :
> Hi Camille,
> You probably don;t need to worry about i=-2, j=-2 if your overlaps are Olx=Oly=2 (in SIZE.h).
> That said, I'm not sure to follow what your are trying to do. You already have a binary bathymetry file and you want the MITgcm to output a netcdf version of the same file?
>
> cheers,
> david
>
> ________________________________________
> From: Camille Mazoyer [mazoyer at univ-tln.fr]
> Sent: Saturday, January 28, 2017 5:58 PM
> To: mitgcm-support at mitgcm.org
> Subject: [MITgcm-support] Error while creating bathymery : ini_masks_etc
>
> Dear all,
>
> I am creating a new bathymetry for a Bay (Toulon Bay, France). My grid
> is a 200x100x130 points grid.
> I want to generate a netcdf bathymetry file for step one: for that, I
> have no OBC, no initial state files.
> I have a binary bathymetry file and binary levels file.
> I use debugging options.
>
> I have a floating invalid problem that I don't understand.
> It's in ini_masks_etc.f, so I decided to modify the code with some
> "write toto" lines:
>
> 65 C o Non-dimensional distance between grid bound. and domain
> lower_R bound.
> 66 hFacCtmp = (rF(k)-R_low(i,j,bi,bj))*recip_drF(k)
> 67 ! cmzy
> 68 write(*,*) 'i,j,k',i,j,k
> 69 write(*,*) 'rF(k)',rF(k)
> 70 write(*,*) 'hFacCtmpR_low(i,j,bi,bj)',R_low(i,j,bi,bj)
> 71 write(*,*) 'recip_drF(k)',recip_drF(k)
> 72 write(*,*) 'hFacCtmp',hFacCtmp
> 73 C o Select between, closed, open or partial (0,1,0-1)
> 74 hFacCtmp=min( max( hFacCtmp, 0. _d 0) , 1. _d 0)
>
> And here is the output:
>
> ** Warning ** CAL_READPARMS: ignores "data.cal" file since
> ** Warning ** CAL_READPARMS: useCAL= F (set from "data.pkg")
> i,j,k -2 -2 1
> rF(k) 0.000000000000000E+000
> hFacCtmpR_low(i,j,bi,bj) NaN
> recip_drF(k) 1.00000000000000
> hFacCtmp NaN
> forrtl: error (65): floating invalid
> Image PC Routine Line Source
> mitgcmuv 0000000001730315 Unknown Unknown Unknown
> mitgcmuv 000000000172DF37 Unknown Unknown Unknown
> mitgcmuv 00000000016E4382 Unknown Unknown Unknown
> mitgcmuv 00000000016E41D6 Unknown Unknown Unknown
> mitgcmuv 000000000167CE25 Unknown Unknown Unknown
> mitgcmuv 00000000016805B7 Unknown Unknown Unknown
> libpthread.so.0 0000003AC5C0F710 Unknown Unknown Unknown
> mitgcmuv 0000000000410D92 ini_masks_etc_ 3151
> ini_masks_etc.f
> mitgcmuv 00000000014EDC49 initialise_fixed_ 1857
> initialise_fixed.f
> mitgcmuv 00000000015F89EA the_model_main_ 2382
> the_model_main.f
> mitgcmuv 000000000113CEAF MAIN__ 3919 main.f
> mitgcmuv 000000000040F41E Unknown Unknown Unknown
> libc.so.6 0000003AC541ED5D Unknown Unknown Unknown
> mitgcmuv 000000000040F329 Unknown Unknown Unknown
>
> Do you have any idea to help me? any test?
> I don't really understand what's going on....
> I must admit that I don't know why there is a point i=-2, j=-2 , and why
> the rfk = 0 there.... I thought my Levels file was ok.
> If somebody has already encounter this error?
>
> Thanks for any advice,
> Camille
>
> --
> ------------------------------------------
> Camille Mazoyer
> Phd Student
> Mediterranean Institute of Oceanography (MIO)
> Institut de Mathématiques de Toulon (IMATH)
> Université de TOULON
> Bat X - CS 60584
> 83041 TOULON cedex 9
> France
> tel: +33.4.94.14.24.50
> Email: mazoyer at univ-tln.fr
> http://mio.pytheas.univ-amu.fr/
> http://imath.fr/
>
>
> _______________________________________________
> 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
--
------------------------------------------
Camille Mazoyer
Phd Student
Mediterranean Institute of Oceanography (MIO)
Institut de Mathématiques de Toulon (IMATH)
Université de TOULON
Bat X - CS 60584
83041 TOULON cedex 9
France
tel: +33.4.94.14.24.50
Email: mazoyer at univ-tln.fr
http://mio.pytheas.univ-amu.fr/
http://imath.fr/
More information about the MITgcm-support
mailing list