Changeset 14738


Ignore:
Timestamp:
2021-04-22T12:47:30+02:00 (2 weeks ago)
Author:
jwhile
Message:

Fix for gulf18

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/Gulf18_changes_for_review/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r14190 r14738  
    167167      LOGICAL :: lk_surft      ! Logical: T => Increments file contains surft variable  
    168168                               !               so only apply surft increments. 
     169      CHARACTER(50) :: cn_surftname   ! Name of the surft variable in the increments file 
    169170      !! 
    170171      NAMELIST/nam_asminc/ ln_bkgwri, ln_avgbkg, ln_balwri,                & 
     
    182183         &                 nitbkg, nitdin, nitiaustr, nitiaufin, niaufn,   & 
    183184         &                 ln_salfix, salfixmin, nn_divdmp, nitavgbkg,     & 
    184          &                 mld_choice, mld_choice_bgc, rn_maxchlinc 
     185         &                 mld_choice, mld_choice_bgc, rn_maxchlinc,       & 
     186         &                 cn_surftname 
     187          
    185188      !!---------------------------------------------------------------------- 
    186189 
     
    211214      niaufn    = 0 
    212215      nitavgbkg = 1 
     216      cn_surftname = "surft" 
    213217 
    214218      REWIND( numnam_ref )              ! Namelist nam_asminc in reference namelist : Assimilation increment 
     
    275279         WRITE(numout,*) '      Logical switch for applying pph increments               ln_pphinc = ', ln_pphinc 
    276280         WRITE(numout,*) '      Logical switch for applying po2 increments               ln_po2inc = ', ln_po2inc 
     281         WRITE(numout,*) '      Surface temperature variable name (if applicable)        cn_surftname = ', cn_surftname 
    277282      ENDIF 
    278283 
     
    520525 
    521526            !Test if the increments file contains the surft variable. 
    522             isurfstat = iom_varid( inum, 'bckinsurft', ldstop = .FALSE. ) 
     527            isurfstat = iom_varid( inum, 'bckin'//TRIM(cn_surftname), ldstop = .FALSE. ) 
    523528            IF ( isurfstat == -1 ) THEN 
    524529               lk_surft = .FALSE. 
     
    526531               lk_surft = .TRUE. 
    527532               CALL ctl_warn( ' Applying 2D temperature increment to bottom of ML: ', & 
    528             &                 ' bckinsurft found in increments file.' ) 
     533            &                 ' bckin'//TRIM(cn_surftname)//' found in increments file.' ) 
    529534            ENDIF              
    530535 
     
    553558                       
    554559               ALLOCATE( t_bkginc_2d(jpi,jpj) )  
    555                CALL iom_get( inum, jpdom_autoglo, 'bckinsurft', t_bkginc_2d, 1)  
     560               CALL iom_get( inum, jpdom_autoglo, 'bckin'//TRIM(cn_surftname), t_bkginc_2d, 1)  
    556561#if defined key_bdy                 
    557562               DO jk = 1,jpkm1  
Note: See TracChangeset for help on using the changeset viewer.