New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4718 for branches – NEMO

Changeset 4718 for branches


Ignore:
Timestamp:
2014-07-15T18:47:11+02:00 (10 years ago)
Author:
vancop
Message:

[sea ice tracer boundary conditions ORCA2_LIM2 PISCES]

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_4707_CNRS04_bgc_ice/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

    r4641 r4718  
    105105         !                                               
    106106         CALL wrk_alloc( jpi, jpj, zironice ) 
     107         zironice(:,:) = 0.0 
    107108         !                                               
    108          DO jj = 1, jpj 
    109             DO ji = 1, jpi 
    110                zdep    = rfact2 / fse3t(ji,jj,1) 
    111                zwflux  = fmmflx(ji,jj) / 1000._wp 
    112                zfminus = MIN( 0._wp, -zwflux ) * trn(ji,jj,1,jpfer) * zdep 
    113                zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep 
    114                zironice(ji,jj) =  zfplus + zfminus 
    115             END DO 
    116          END DO 
    117          ! 
    118          trn(:,:,1,jpfer) = trn(:,:,1,jpfer) + zironice(:,:)  
    119          !                                               
     109! MV NOW DONE IN trcsbc 
     110!        DO jj = 1, jpj 
     111!           DO ji = 1, jpi 
     112!              zdep    = rfact2 / fse3t(ji,jj,1) 
     113!              zwflux  = fmmflx(ji,jj) / 1000._wp 
     114!              zfminus = MIN( 0._wp, -zwflux ) * trn(ji,jj,1,jpfer) * zdep 
     115!              zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep 
     116!              zironice(ji,jj) =  zfplus + zfminus 
     117!           END DO 
     118!        END DO 
     119!        ! 
     120!        trn(:,:,1,jpfer) = trn(:,:,1,jpfer) + zironice(:,:)  
     121!        !                                               
    120122         IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc )   & 
    121123            &   CALL iom_put( "Ironice", zironice(:,:) * 1.e+3 * rfact2r * fse3t(:,:,1) * tmask(:,:,1) ) ! iron flux from ice 
    122124         CALL wrk_dealloc( jpi, jpj, zironice ) 
    123          !                                               
    124       ENDIF 
     125!        !                                               
     126!     ENDIF 
     127! END MV 
    125128 
    126129      ! Add the external input of nutrients from dust deposition 
Note: See TracChangeset for help on using the changeset viewer.