WC13 compilation error of Normalization

[dickson]

Hi,

I got the following errors during compilation of Normalization under the test case of WC13 when I add a CPP option "define NESTING" in the WC13.h. The compilation was successful using the default wc13 header file without adding the NESTING. May I ask for your help and how to fix ? Thank you

ROMS SVN 1151

cd /home/my_test/WC13/Normalization/Build_roms; /opt/intel/compilers/linux/mpi/intel64/bin/mpiifort -c -fp-model precise -heap-arrays -ip -O3 -traceback -check uninit set_contact.f90
set_contact.f90(608): error #6460: This is not a field name that is defined in the encompassing structure. [TL_VWEIGHT]
allocate ( Rcontact(cr) % tl_Vweight(2,N(dg),NpointsR(cr)) )
----------------------------------^
set_contact.f90(614): error #6460: This is not a field name that is defined in the encompassing structure. [AD_VWEIGHT]
allocate ( Rcontact(cr) % ad_Vweight(2,N(dg),NpointsR(cr)) )
----------------------------------^
set_contact.f90(915): error #6460: This is not a field name that is defined in the encompassing structure. [TL_BUSTR]
allocate ( COMPOSITE(cr) % tl_bustr(4,NpointsU(cr)) )
-------------------------------------^
set_contact.f90(916): error #6460: This is not a field name that is defined in the encompassing structure. [TL_BVSTR]
allocate ( COMPOSITE(cr) % tl_bvstr(4,NpointsV(cr)) )
-------------------------------------^
set_contact.f90(917): error #6460: This is not a field name that is defined in the encompassing structure. [TL_UBAR]
allocate ( COMPOSITE(cr) % tl_ubar(4,NpointsU(cr),2) )
-------------------------------------^
set_contact.f90(918): error #6460: This is not a field name that is defined in the encompassing structure. [TL_VBAR]
allocate ( COMPOSITE(cr) % tl_vbar(4,NpointsV(cr),2) )
-------------------------------------^
set_contact.f90(919): error #6460: This is not a field name that is defined in the encompassing structure. [TL_ZETA]
allocate ( COMPOSITE(cr) % tl_zeta(4,NpointsR(cr),2) )
-------------------------------------^
set_contact.f90(920): error #6460: This is not a field name that is defined in the encompassing structure. [TL_RZETA]
allocate ( COMPOSITE(cr) % tl_rzeta(4,NpointsR(cr)) )
-------------------------------------^
set_contact.f90(924): error #6460: This is not a field name that is defined in the encompassing structure. [AD_BUSTR]
allocate ( COMPOSITE(cr) % ad_bustr(4,NpointsU(cr)) )
-------------------------------------^
set_contact.f90(925): error #6460: This is not a field name that is defined in the encompassing structure. [AD_BVSTR]
allocate ( COMPOSITE(cr) % ad_bvstr(4,NpointsV(cr)) )
-------------------------------------^
set_contact.f90(926): error #6460: This is not a field name that is defined in the encompassing structure. [AD_UBAR]
allocate ( COMPOSITE(cr) % ad_ubar(4,NpointsU(cr),2) )
-------------------------------------^
set_contact.f90(927): error #6460: This is not a field name that is defined in the encompassing structure. [AD_VBAR]
allocate ( COMPOSITE(cr) % ad_vbar(4,NpointsV(cr),2) )
-------------------------------------^
set_contact.f90(928): error #6460: This is not a field name that is defined in the encompassing structure. [AD_ZETA]
allocate ( COMPOSITE(cr) % ad_zeta(4,NpointsR(cr),2) )
-------------------------------------^
set_contact.f90(929): error #6460: This is not a field name that is defined in the encompassing structure. [AD_RZETA]
allocate ( COMPOSITE(cr) % ad_rzeta(4,NpointsR(cr)) )
-------------------------------------^
set_contact.f90(947): error #6460: This is not a field name that is defined in the encompassing structure. [TL_DU_AVG1]
allocate ( COMPOSITE(cr) % tl_DU_avg1(4,NpointsU(cr)) )
-------------------------------------^
set_contact.f90(948): error #6460: This is not a field name that is defined in the encompassing structure. [TL_DV_AVG1]
allocate ( COMPOSITE(cr) % tl_DV_avg1(4,NpointsV(cr)) )
-------------------------------------^
set_contact.f90(949): error #6460: This is not a field name that is defined in the encompassing structure. [TL_ZT_AVG1]
allocate ( COMPOSITE(cr) % tl_Zt_avg1(4,NpointsR(cr)) )
-------------------------------------^
set_contact.f90(953): error #6460: This is not a field name that is defined in the encompassing structure. [TL_U]
allocate ( COMPOSITE(cr) % tl_u(4,N(dg),NpointsU(cr)) )
-------------------------------------^
set_contact.f90(954): error #6460: This is not a field name that is defined in the encompassing structure. [TL_V]
allocate ( COMPOSITE(cr) % tl_v(4,N(dg),NpointsV(cr)) )
-------------------------------------^
set_contact.f90(955): error #6460: This is not a field name that is defined in the encompassing structure. [TL_HUON]
allocate ( COMPOSITE(cr) % tl_Huon(4,N(dg),NpointsU(cr)) )
-------------------------------------^
set_contact.f90(956): error #6460: This is not a field name that is defined in the encompassing structure. [TL_HVOM]
allocate ( COMPOSITE(cr) % tl_Hvom(4,N(dg),NpointsV(cr)) )
-------------------------------------^
set_contact.f90(959): error #6460: This is not a field name that is defined in the encompassing structure. [TL_T]
allocate ( COMPOSITE(cr) % tl_t(4,N(dg),NpointsR(cr),NT(dg)) )
-------------------------------------^
set_contact.f90(961): error #6460: This is not a field name that is defined in the encompassing structure. [AD_DU_AVG1]
allocate ( COMPOSITE(cr) % ad_DU_avg1(4,NpointsU(cr)) )
-------------------------------------^
set_contact.f90(962): error #6460: This is not a field name that is defined in the encompassing structure. [AD_DV_AVG1]
allocate ( COMPOSITE(cr) % ad_DV_avg1(4,NpointsV(cr)) )
-------------------------------------^
set_contact.f90(963): error #6460: This is not a field name that is defined in the encompassing structure. [AD_ZT_AVG1]
allocate ( COMPOSITE(cr) % ad_Zt_avg1(4,NpointsR(cr)) )
-------------------------------------^
set_contact.f90(967): error #6460: This is not a field name that is defined in the encompassing structure. [AD_U]
allocate ( COMPOSITE(cr) % ad_u(4,N(dg),NpointsU(cr)) )
-------------------------------------^
set_contact.f90(968): error #6460: This is not a field name that is defined in the encompassing structure. [AD_V]
allocate ( COMPOSITE(cr) % ad_v(4,N(dg),NpointsV(cr)) )
-------------------------------------^
set_contact.f90(969): error #6460: This is not a field name that is defined in the encompassing structure. [AD_HUON]
allocate ( COMPOSITE(cr) % ad_Huon(4,N(dg),NpointsU(cr)) )
-------------------------------------^
set_contact.f90(970): error #6460: This is not a field name that is defined in the encompassing structure. [AD_HVOM]
allocate ( COMPOSITE(cr) % ad_Hvom(4,N(dg),NpointsV(cr)) )
-------------------------------------^
set_contact.f90(973): error #6460: This is not a field name that is defined in the encompassing structure. [AD_T]
allocate ( COMPOSITE(cr) % ad_t(4,N(dg),NpointsR(cr),NT(dg)) )

[arango]

We haven't released the Nested 4D-Var algorithms yet for you to activate the NESTING flag in any of the adjoint-based data assimilation algorithms. Nesting has nothing to do with the computation of the error covariance normalization coefficients. They are computed from the coarser nested grid instead. There is absolutely no reason for you to activate NESTING. It is only possible to activate nesting in the forward nonlinear trajectory using the 4D-Var split scheme, but I suspect this algorithm is very advanced for a regular user. It requires expertise, and we haven't set up a test case for users to see how it is done.