Dear all,
I have set up ORCA2_ICE_PISCES in NEMO v4.0.6 using gfortran compiler.
I met two significant errors for running model, which are the instability error and the p4zopt.f90 crash issue. I solved the problem temporarily by controlling the time interval (rn_rdt).
First, encountered the instability error, I reduced the time interval of model (rn_rdt).
Second, in “p4zopt.f90 crash issue”, I slightly increased the time step near the point.
However, in 2012, there was a serious trouble where two problems occurred at the same time on the same day. (I hope that the model runs by 2021.) I don’t know how to solve this problem.
Any suggestions and advices would be appreciated. Experimental setup and related files are below.
Thank you in advance.
Best regards,
Hwa-Jin Choi
1. Instability error
ocean.output_0052
===>>> : E R R O R
STOP
stp_ctl: |ssh| > 20 m or |U| > 10 m/s or S <= 0 or S >= 100 or NaN encounter in the tests
kt 394 |ssh| max 0.000 at i j 182 109 MPI rank 52
kt 394 |U| max 0.000 at i j k 182 109 1 MPI rank 52
kt 394 Sal min NaN at i j k 163 106 1 MPI rank 52
kt 394 Sal max NaN at i j k 163 106 1 MPI rank 52
=> output of last computed fields in output.abort* files
huge E-R-R-O-R : immediate stop
2. p4zopt.f90 crash issue
I receive the following message
Error termination. Backtrace:
At line 132 of file /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zopt.f90
Fortran runtime error: Index '0' of dimension 2 of array 'rkrgb' below lower bound of 1
Error termination. Backtrace:
#0 0x1529e9010d21 in ???
#1 0x1529e9011869 in ???
#2 0x1529e9011ee6 in ???
#3 0x55c6340a63a2 in __p4zopt_MOD_p4z_opt
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zopt.f90:132
#4 0x55c634adcbf2 in __p4zbio_MOD_p4z_bio
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zbio.f90:121
#5 0x55c6341b4d8d in __p4zsms_MOD_p4z_sms
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zsms.f90:163
#6 0x55c63458ac6b in __trcsms_pisces_MOD_trc_sms_pisces
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcsms_pisces.f90:81
#7 0x55c63458611a in __trcsms_MOD_trc_sms
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcsms.f90:92
#8 0x55c63326ab27 in __trcstp_MOD_trc_stp
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcstp.f90:144
#9 0x55c63307c9a9 in __step_MOD_stp
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/step.f90:237
#10 0x55c632fd8b7c in __nemogcm_MOD_nemo_gcm
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/nemogcm.f90:192
#11 0x55c632fd3f91 in nemo
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/WORK/nemo.f90:18
#12 0x55c632fd3fca in main
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/WORK/nemo.f90:11
Experimental Setup
- JRA55 reanalysis data was prescribed to surface boundary condition.
- The COARE 3.5 algorithm in bulk formulation was used to calculate fluxes from the atmosphere.
- Atmospheric pCO2 was changed as an annual observation data.
- The monthly satellite observation for chlorophyll was used.
→ The model had run from 1960 to 2012. I hope that the model runs by 2021.
My arch file:
%NCDF_INC -I%NCDF_HOME/include -I%HDF5_HOME/include
%NCDF_LIB -L%NCDF_HOME -lnetcdff -lnetcdf -lstdc++
%XIOS_INC -I%XIOS_HOME/inc
%XIOS_LIB -L%XIOS_HOME/lib -lxios -L/usr/lib/gcc/x86_64-linux-gnu/9 -lstdc++
%OASIS_INC -I%OASIS_HOME/build/lib/mct -I%OASIS_HOME/build/lib/psmile.MPI1
%OASIS_LIB -L%OASIS_HOME/lib -lpsmile.MPI1 -lmct -lmpeu -lscrip
%CPP /usr/bin/cpp-9 -Dkey_nosignedzero
%CPPFLAGS -P -traditional
%FC /media/cmlws/Data1/hjc/Tools/install/bin/mpif90
%FCFLAGS -fdefault-real-8 -funroll-all-loops -cpp -fcray-pointer -ffree-line-length-none -g -O0 -fbacktrace
%FFLAGS %FCFLAGS
%LD /media/cmlws/Data1/hjc/Tools/install/bin/mpif90
%LDFLAGS
%FPPFLAGS -P -C -traditional
%AR ar
%ARFLAGS -rs
%MK make
%USER_INC %XIOS_INC %OASIS_INC %NCDF_INC
%USER_LIB %XIOS_LIB %OASIS_LIB %NCDF_LIB
%CC gcc
%CFLAGS -O0 -fbacktrace