Ignore:
Timestamp:
2019-03-20T18:53:44+01:00 (19 months ago)
Author:
aumont
Message:

bug fixes corresponding to tickets #2262

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/TOP/PISCES/trcice_pisces.F90

    r10213 r10780  
    266266      zpisc(jpgon,4) =  1.05e-8_wp 
    267267      zpisc(jpgop,4) =  1.05e-8_wp 
    268   
     268! 
     269!     ln_ironice and tracers in seaice are redundant. Thus, if tracers in ice 
     270!     is activated, ln_ironice should be set to false 
     271!     ------------------------------------------------------------------------ 
     272      IF( nn_ice_tr >= 0 .AND. ln_ironice ) THEN 
     273         IF(lwp) THEN 
     274            WRITE(numout,*) '   ==>>>   ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 
     275            WRITE(numout,*) '           Specify your sea ice iron concentration in nampisice instead ' 
     276            WRITE(numout,*) '           ln_ironice is forced to .FALSE. ' 
     277         ENDIF 
     278         ln_ironice = .FALSE. 
     279      ENDIF 
     280!  
    269281      DO jn = jp_pcs0, jp_pcs1 
    270282         IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1)  ! Global case 
     
    310322            trc_i(:,:,jn) = trc_ice_prescr(jn) 
    311323         ENDIF 
    312         
    313324         !-- Baltic 
    314325         IF( cn_cfg == "orca" .OR. cn_cfg == "ORCA" ) THEN      
     
    318329                     trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn)  
    319330               END WHERE 
    320             ELSE                                 ! prescribed tracer concentration in ice 
    321                WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    322                    54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    323                      trc_i(:,:,jn) = trc_ice_prescr(jn) 
    324                END WHERE 
    325             ENDIF ! trc_ice_ratio 
     331            ENDIF 
    326332         ENDIF 
    327333      ! 
Note: See TracChangeset for help on using the changeset viewer.