ad passive tracer

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

ad passive tracer

#1 Unread post by thaitlhp »

i'm trying ad passive tracer for upwelling test case step by step following:

viewtopic.php?f=14&t=2611

i added dye_01 in rst.nc and set ININAME= rst.nc then i set NPT = 1
but when i run uwelling.in i get this:
mpirun -np 3 ./romsM roms_upwelling.in
--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.7
Thursday - May 23, 2019 - 1:16:02 PM
--------------------------------------------------------------------------------
At line 2008 of file mod_ncparam.f90
Fortran runtime error: End of record

Error termination. Backtrace:
At line 2008 of file mod_ncparam.f90
Fortran runtime error: End of record

What is this error?
Error termination. Backtrace:
At line 2008 of file mod_ncparam.f90
Fortran runtime error: End of record

Error termination. Backtrace:
#0 0x7f54f7ed22da in ???
#1 0x7f54f7ed2ec5 in ???
#2 0x7f54f7ed368d in ???
#3 0x7f54f8049060 in ???
#4 0x7f54f80549e9 in ???
#5 0x7f54f804c2cc in ???
#6 0x7f54f804c46c in ???
#7 0x5643a8ed4552 in ???
#8 0x5643a8ce7c71 in ???
#9 0x5643a8c85a39 in ???
#10 0x5643a8c81fd4 in ???
#11 0x5643a8c818de in ???
#12 0x5643a8c8164e in ???
#13 0x7f54f7749b96 in ???
#14 0x5643a8c816a9 in ???
#15 0xffffffffffffffff in ???
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
#0 0x7f6c4b0a02da in ???
#1 0x7f6c4b0a0ec5 in ???
#2 0x7f6c4b0a168d in ???
#3 0x7f6c4b217060 in ???
#4 0x7f6c4b2229e9 in ???
#5 0x7f6c4b21a2cc in ???
#6 0x7f6c4b21a46c in ???
#7 0x5583bb66b552 in ???
#8 0x5583bb47ec71 in ???
#9 0x5583bb41ca39 in ???
#10 0x5583bb418fd4 in ???
#11 0x5583bb4188de in ???
#12 0x5583bb41864e in ???
#13 0x7f6c4a917b96 in ???
#14 0x5583bb4186a9 in ???
#15 0xffffffffffffffff in ???
#0 0x7ff523fea2da in ???
#1 0x7ff523feaec5 in ???
#2 0x7ff523feb68d in ???
#3 0x7ff524161060 in ???
#4 0x7ff52416c9e9 in ???
#5 0x7ff5241642cc in ???
#6 0x7ff52416446c in ???
#7 0x55dd42f1d552 in ???
#8 0x55dd42d30c71 in ???
#9 0x55dd42ccea39 in ???
#10 0x55dd42ccafd4 in ???
#11 0x55dd42cca8de in ???
#12 0x55dd42cca64e in ???
#13 0x7ff523861b96 in ???
#14 0x55dd42cca6a9 in ???
#15 0xffffffffffffffff in ???
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

Process name: [[64664,1],1]
Exit code: 2
what is
Attachments
upwelling.h
(1.73 KiB) Downloaded 202 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: ad passive tracer

#2 Unread post by kate »

You've got one clue here:

Code: Select all

At line 2008 of file mod_ncparam.f90
Can you tell us what is at line 2008 of your mod_ncparam.f90 file?

What is NRREC in your roms.in file? The restart file is not necessarily the same as a file needed for initialization, especially when you go changing settings.

thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

Re: ad passive tracer

#3 Unread post by thaitlhp »

kate wrote:You've got one clue here:

Code: Select all

At line 2008 of file mod_ncparam.f90
Can you tell us what is at line 2008 of your mod_ncparam.f90 file?

What is NRREC in your roms.in file? The restart file is not necessarily the same as a file needed for initialization, especially when you go changing settings.
thanks for your reply!
now, i don't use perfect_restart, i set NRREC == 0 run again and get:
At line 1860 of file mod_ncparam.f90
Fortran runtime error: End of record
my mod_ncparam.f90 here
Attachments
mod_ncparam.f90
(84.08 KiB) Downloaded 195 times

thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

Re: ad passive tracer

#4 Unread post by thaitlhp »

kate wrote:You've got one clue here:

Code: Select all

At line 2008 of file mod_ncparam.f90
Can you tell us what is at line 2008 of your mod_ncparam.f90 file?

What is NRREC in your roms.in file? The restart file is not necessarily the same as a file needed for initialization, especially when you go changing settings.


my input file
Attachments
roms_upwelling.in
(132.31 KiB) Downloaded 184 times

thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

Re: ad passive tracer

#5 Unread post by thaitlhp »

kate wrote:You've got one clue here:

Code: Select all

At line 2008 of file mod_ncparam.f90
Can you tell us what is at line 2008 of your mod_ncparam.f90 file?

What is NRREC in your roms.in file? The restart file is not necessarily the same as a file needed for initialization, especially when you go changing settings.
my header file
Looking forward to your help! :D :D :D :D
Attachments
upwelling.h
(1.09 KiB) Downloaded 195 times

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: ad passive tracer

#6 Unread post by wilkin »

Your roms_upwelling.in has 4 passive tracers, not one. You have set NPT=4.

Your roms_upwelling.in has also set "closed" as the options for open boundary conditions on passive tracers:

Code: Select all

   LBC(isTvar) ==   Per     Clo     Per     Clo \       ! temperature
                    Per     Clo     Per     Clo \
Clo     Clo     Clo     Clo \
Clo     Clo     Clo     Clo \
Clo     Clo     Clo     Clo \        ! salinity
		    Clo     Clo     Clo     Clo
when that should be periodic to match the active tracers.

You have options I don't recognize like ANA_PTOBC. Who's code are you using? It's not the default myroms.org code.

For these test purposes I think you could perhaps set initial values with ANA_PASSIVE and avoid creating the netcdf file until you get started.

I see you have AGE_MEAN. So you need to think carefully about initializing and forcing. Set a pair of age at initialization only will simply give you a patch of uniformly aging water.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

Re: ad passive tracer

#7 Unread post by thaitlhp »

wilkin wrote:Your roms_upwelling.in has 4 passive tracers, not one. You have set NPT=4.

Your roms_upwelling.in has also set "closed" as the options for open boundary conditions on passive tracers:

Code: Select all

   LBC(isTvar) ==   Per     Clo     Per     Clo \       ! temperature
                    Per     Clo     Per     Clo \
Clo     Clo     Clo     Clo \
Clo     Clo     Clo     Clo \
Clo     Clo     Clo     Clo \        ! salinity
		    Clo     Clo     Clo     Clo
when that should be periodic to match the active tracers.

You have options I don't recognize like ANA_PTOBC. Who's code are you using? It's not the default myroms.org code.

For these test purposes I think you could perhaps set initial values with ANA_PASSIVE and avoid creating the netcdf file until you get started.

I see you have AGE_MEAN. So you need to think carefully about initializing and forcing. Set a pair of age at initialization only will simply give you a patch of uniformly aging water.

Thank you very much for your suggestion! first i ran with NPT=1, after i tried with NPT=4,from

Code: Select all

  Per     Clo     Per     Clo 
to

Code: Select all

   Clo     Clo     Clo     Clo 
, all case get problem with mod_ncparam.F
about ANA_PTOBC Define for analytic passive tracer boundary conditions i read it in roms_manual.
where can i get ad passive tracer suggestion more clearly? forgive me, my English skill is still bad!I have been learning about ROMS for 1 year but My knowledge of ROMS is still small! thank all!
Last edited by thaitlhp on Sat May 25, 2019 9:54 am, edited 1 time in total.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: ad passive tracer

#8 Unread post by kate »

Which mod_ncparam.f90 error are you getting and what is on the line number in question? We want you to look at the file. Maybe you should also explain exactly what it is you want to do with your passive tracer.

thaitlhp
Posts: 15
Joined: Thu Apr 05, 2018 6:34 pm
Location: Ha Noi University Sience and Technology, Vietnam

Re: ad passive tracer

#9 Unread post by thaitlhp »

kate wrote:Which mod_ncparam.f90 error are you getting and what is on the line number in question? We want you to look at the file. Maybe you should also explain exactly what it is you want to do with your passive tracer.
Thank Kate!

Code: Select all

 
      SELECT CASE (Vinfo(1))
              CASE ('tracer2')
                varid=varid-1
                DO i=1,MT
                  varid=varid+1
                  idTTav(i)=varid
                  DO ng=1,Ngrids
                    Fscale(varid,ng)=scale
                    Iinfo(1,varid,ng)=gtype
                  END DO
                  WRITE (Vname(1,varid),'(a,a)')                        &
     &                  TRIM(ADJUSTL(Vname(1,idTvar(i)))), '_2'
                  WRITE (Vname(2,varid),'(a,1x,a)')                     &
     &                  'time-averaged squared',                        &
     &                  TRIM(ADJUSTL(Vname(2,idTvar(i))))
                  IF (TRIM(ADJUSTL(Vname(3,idTvar(i)))).eq.             &
     &                'nondimensional') THEN
                    WRITE (Vname(3,varid),'(a)')                        &
     &                    TRIM(ADJUSTL(Vname(3,idTvar(i))))
                  ELSE
                    WRITE (Vname(3,varid),'(a,a)')                      &
     &                    TRIM(ADJUSTL(Vname(3,idTvar(i)))), '2'
                  END IF
                  WRITE (Vname(4,varid),'(a,a)')                        &
     &                  TRIM(Vname(1,varid)), ', scalar, series'
                  WRITE (Vname(5,varid),'(a)')                          &
     &                  TRIM(ADJUSTL(Vinfo(5)))
                END DO
I want ad Radioactive function of Cs-137 with passive tracer. But first, i need understand how passive tracer work so i set default in ana_passive.h and looking for dye in output file.
:D :D :D
Attachments
ana_passive.h
(3.35 KiB) Downloaded 194 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: ad passive tracer

#10 Unread post by kate »

Are you by chance asking for Aout(idTTav)? I have never done so and can't guarantee that it would work.

heroh2o
Posts: 1
Joined: Tue Sep 20, 2011 2:27 pm
Location: TianJin University

Re: ad passive tracer

#11 Unread post by heroh2o »

I'm trying to add a passive tracer in upwelling test case, and I have the same problem.
If I edit passive tracer part of “varinfo.dat” file, it shows:

At line 1904 of file mod_ncparam.f90
Fortran runtime error: End of record

If I use the original “varinfo.dat” file, it shows:

Fortran runtime error: shape and target do not conform

Have you solved your problem?
I need some suggestion.Thanks!

Post Reply