Changeset 5390 for branches/UKMO/dev_r5377_UKMO1_CICE_coupling_GSI7/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2015-06-10T10:48:39+02:00 (9 years ago)

Author:

jamrae

Message:

Changes required for GSI7 in coupled configuration.

File:

• branches/UKMO/dev_r5377_UKMO1_CICE_coupling_GSI7/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (6 diffs)

branches/UKMO/dev_r5377_UKMO1_CICE_coupling_GSI7/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -16 +16 @@
    USE domvvl
    USE eosbn2, only : eos_fzp ! Function to calculate freezing point of seawater
-   USE phycst, only : rcp, rau0, r1_rau0, rhosn, rhoic
+   USE phycst, only : rcp, rau0, r1_rau0, rhosn, rhoic, rt0
    USE in_out_manager  ! I/O manager
    USE iom, ONLY : iom_put,iom_use              ! I/O manager library !!Joakim edit
@@ -39 +39 @@
    USE ice_calendar, only: dt
 # if defined key_cice4
-   USE ice_state, only: aice,aicen,uvel,vvel,vsno,vsnon,vice,vicen
+   USE ice_state, only: aice,aicen,uvel,vvel,vsno,vsnon,vice,vicen, &
+                        nt_Tsfc
    USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow,  &
                 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_gbm,     &
@@ -45 +46 @@
                 flatn_f,fsurfn_f,fcondtopn_f,                    &
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
-                swvdr,swvdf,swidr,swidf,Tf
+                swvdr,swvdf,swidr,swidf,Tf              !!!   &
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!! LINE BELOW COMMENTED OUT TEMPORARILY UNTIL CICE MULTILAYERS  !!!
+!!! BRANCH CAN BE INCLUDED AT CICE5.                             !!!
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!              ,keffn_top,Tn_top
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
+
    USE ice_therm_vertical, only: calc_Tsfc
 #else
@@ -56 +65 @@
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
                 swvdr,swvdf,swidr,swidf,Tf
-   USE ice_therm_shared, only: calc_Tsfc
+   USE ice_therm_shared, only: calc_Tsfc, heat_capacity
    USE ice_shortwave, only: apeffn
 #endif
@@ -453 +462 @@
       CALL nemo2cice(sss_m,sss,'T', 1. )
 
+! Sea ice surface skin temperature
+      DO jpl=1,ncat
+        CALL nemo2cice(tsfc_ice(:,:,jpl), trcrn(:,:,nt_tsfc,jpl,:),'T',1.)
+      ENDDO 
+
 ! x comp and y comp of surface ocean current
 ! U point to F point
@@ -727 +741 @@
       ENDDO
 #endif
+
+
+! If using multilayers thermodynamics in CICE then get top layer temperature
+! and effective conductivity       
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!!! LINES COMMENTED OUT UNTIL CICE MULTILAYERS BRANCH IS IMPLEMENTED 
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
+!!      IF (heat_capacity) THEN
+!!         DO jpl = 1,ncat
+!!            CALL cice2nemo(Tn_top(:,:,jpl,:),tn_ice(:,:,jpl),'T', 1. )
+!!            CALL cice2nemo(keffn_top(:,:,jpl,:),kn_ice(:,:,jpl),'T', 1. )
+!!         ENDDO
+!!! Convert surface temperature to Kelvin
+!!         tn_ice(:,:,:)=tn_ice(:,:,:)+rt0
+!!      ELSE
+!!         tn_ice(:,:,:) = 0.0
+!!         kn_ice(:,:,:) = 0.0
+!!      ENDIF       
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
 
       !