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 5624 for branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90 – NEMO

Ignore:
Timestamp:
2015-07-21T20:47:42+02:00 (9 years ago)
Author:
mathiot
Message:

UKMO_ISF : fix conservation issue based on the work of Jerome on runoff, simplification of trasbc (isf part only) and remove option to apply isf melting as volume flux or not

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r5621 r5624  
    266266               ! Add volume filter correction: compatibility with tracer advection scheme 
    267267               ! => time filter + conservation correction (only at the first level) 
    268                fse3t_b(:,:,1) = fse3t_b(:,:,1) - atfp * rdt * r1_rau0 * ( emp_b(:,:) - emp(:,:) & 
    269                               &                                          -rnf_b(:,:) + rnf(:,:) ) * tmask(:,:,1) 
    270                !DO jj = 1,jpj 
    271                !   DO ji = 1,jpi 
    272                !      jk = mikt(ji,jj) 
    273                !      fse3t_b(ji,jj,jk) = fse3t_b(ji,jj,jk) - atfp * rdt * r1_rau0  
    274                !                        &                          * ( (emp_b(ji,jj) - emp(ji,jj) ) & 
    275                !                        &                            - (rnf_b(ji,jj) - rnf(ji,jj) ) &     
    276                !                        &                            + (fwfisf_b(ji,jj) - fwfisf(ji,jj) ) ) * tmask(ji,jj,jk) ! ????? 
    277                !   END DO 
    278                !END DO 
     268               DO jj = 1,jpj 
     269                  DO ji = 1,jpi 
     270                     jk = mikt(ji,jj) 
     271                     fse3t_b(ji,jj,jk) = fse3t_b(ji,jj,jk) - atfp * rdt * r1_rau0                  &  
     272                                       &                          * ( (emp_b(ji,jj) - emp(ji,jj) ) & 
     273                                       &                            - (rnf_b(ji,jj) - rnf(ji,jj) ) &     
     274                                       &                            + (fwfisf_b(ji,jj) - fwfisf(ji,jj) ) ) * tmask(ji,jj,jk) 
     275                  END DO 
     276               END DO 
    279277            ENDIF 
    280278            ! 
Note: See TracChangeset for help on using the changeset viewer.