Changeset 11809


Ignore:
Timestamp:
2019-10-26T14:15:48+02:00 (13 months ago)
Author:
smueller
Message:

Inclusion of river runoff (arrays rnf and rnf_b of module sbc_oce) in the nonlinear-model trajectory used by NEMOTAM; elimination of the use of an uninitialised array (array rnf_b of module sbc_oce) in the reconstruction of the vertical velocity field of the nonlinear-model trajectory used by NEMOTAM (fix for the bugs described in ticket #1876)

Location:
NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPATAM_SRC/TAM/trj_tam.F90

    r11774 r11809  
    113113      & hmlp1,    & 
    114114      & hmlp2 
     115  REAL(wp), ALLOCATABLE, DIMENSION(:,:), SAVE :: & 
     116      & rnfr1,    & 
     117      & rnfr2,    & 
     118      & rnf_br1,  & 
     119      & rnf_br2 
    115120 
    116121CONTAINS 
     
    290295                  & ) 
    291296#endif 
     297               ALLOCATE( & 
     298                  & rnfr1(jpi, jpj),       & 
     299                  & rnfr2(jpi, jpj),       & 
     300                  & rnf_br1(jpi, jpj),     & 
     301                  & rnf_br2(jpi, jpj)) 
    292302               ln_mem = .TRUE. 
    293303            ENDIF 
     
    351361#endif 
    352362#endif 
     363            rnfr2   (:,:)   = rnfr1   (:,:) 
     364            rnf_br2 (:,:)   = rnf_br1 (:,:) 
    353365 
    354366            istpr1 = INT( stpr1 ) 
     
    425437            CALL iom_get( inumtrj1, jpdom_autoglo, 'aeiw'  , aeiwr1  , inrc ) 
    426438#endif 
     439            CALL iom_get( inumtrj1, jpdom_autoglo, 'rnf'   , rnfr1   , inrc ) 
     440            CALL iom_get( inumtrj1, jpdom_autoglo, 'rnf_b' , rnf_br1 , inrc ) 
    427441            CALL iom_close( inumtrj1 ) 
    428442 
     
    484498#endif 
    485499#endif 
     500            rnfr1   (:,:)   = rnfr2   (:,:) 
     501            rnf_br1 (:,:)   = rnf_br2 (:,:) 
    486502 
    487503            istpr1 = INT( stpr1 ) 
     
    547563            CALL iom_get( inumtrj2, jpdom_autoglo, 'aeiw'  , aeiwr2  , inrc ) 
    548564#endif 
     565            CALL iom_get( inumtrj2, jpdom_autoglo, 'rnf'   , rnfr2   , inrc ) 
     566            CALL iom_get( inumtrj2, jpdom_autoglo, 'rnf_b' , rnf_br2 , inrc )             
    549567            CALL iom_close( inumtrj2 ) 
    550568 
     
    644662#endif 
    645663#endif 
     664      rnf(:,:)      = zwtr1 * rnfr1   (:,:)   + zwtr2 * rnfr2   (:,:) 
     665      rnf_b(:,:)    = zwtr1 * rnf_br1 (:,:)   + zwtr2 * rnf_br2 (:,:) 
    646666 
    647667      CALL ssh_wzv( kstp ) 
     
    927947               & ) 
    928948#endif 
     949            DEALLOCATE(    & 
     950               & rnfr1,    & 
     951               & rnfr2,    & 
     952               & rnf_br1,  & 
     953               & rnf_br2   & 
     954               & ) 
    929955    ENDIF 
    930956         END SUBROUTINE trj_deallocate 
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/TAM/tamtrj.F90

    r3604 r11809  
    211211         CALL iom_rstput( it, it, inum, 'aeiv'  , aeiv, ktype = ntype ) 
    212212         CALL iom_rstput( it, it, inum, 'aeiw'  , aeiw, ktype = ntype ) 
     213         CALL iom_rstput( it, it, inum, 'rnf'   , rnf,  ktype = ntype ) 
     214         CALL iom_rstput( it, it, inum, 'rnf_b' , rnf_b, ktype = ntype ) 
    213215 
    214216         CALL iom_close( inum ) 
Note: See TracChangeset for help on using the changeset viewer.