Changeset 3402 for branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2012-05-25T18:43:49+02:00 (12 years ago)

Author:

acc

Message:

Branch: dev_r3385_NOCS04_HAMF; #665. Stage 2 of 2012 development: suppression of emps array and introduction of sfx (salt flux) array with associated code to setup the options for embedding the seaice into the ocean

File:

• branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (5 diffs)

branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -107 +107 @@
       !! ** Action  : - time evolution of the CICE sea-ice model
       !!              - update all sbc variables below sea-ice:
-      !!                utau, vtau, qns , qsr, emp , emps
+      !!                utau, vtau, qns , qsr, emp , sfx 
       !!---------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt      ! ocean time step
@@ -480 +480 @@
 ! salinity calculation
       CALL cice2nemo(fresh_gbm,ztmp,'T', 1. )
-      emps(:,:)=emp(:,:)-ztmp(:,:)
+      sfx (:,:)=emp(:,:)-ztmp(:,:)
 ! Note the 1000.0 is to convert from kg salt to g salt (needed for PSU)
       CALL cice2nemo(fsalt_gbm,ztmp,'T', 1. )
@@ -486 +486 @@
          DO ji=1,jpi
             IF (sss_m(ji,jj).gt.0.0) THEN
-               emps(ji,jj)=emps(ji,jj)+ztmp(ji,jj)*1000.0/sss_m(ji,jj)
+               sfx (ji,jj)=sfx (ji,jj)+ztmp(ji,jj)*1000.0/sss_m(ji,jj)
             ENDIF
          ENDDO
@@ -494 +494 @@
 ! weight of the precip will affect the free surface even if it falls on the ice
 ! (same to the argument that freezing / melting of ice doesn't change the free surface) 
-! Sublimation from the ice is treated in a similar way (included in emp but not emps)  
+! Sublimation from the ice is treated in a similar way (included in emp but not sfx )  
 !
 ! This should not be done in the variable volume case
@@ -512 +512 @@
 
       CALL lbc_lnk( emp , 'T', 1. )
-      CALL lbc_lnk( emps , 'T', 1. )
+      CALL lbc_lnk( sfx  , 'T', 1. )
 
 ! Solar penetrative radiation and non solar surface heat flux