Changeset 1482 for trunk/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2009-07-03T17:28:06+02:00 (15 years ago)

Author:

smasson

Message:

distribution of iom_put + cleaning of LIM2 outputs, see ticket:437

File:

• trunk/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (6 diffs)

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -25 +25 @@
    USE lbclnk           ! ocean lateral boundary condition
    USE in_out_manager   ! I/O manager
+   USE iom              ! 
    USE albedo           ! albedo parameters
    USE prtctl           ! Print control
@@ -92 +93 @@
       REAL(wp) ::   zsang, zmod, zfm
       REAL(wp), DIMENSION(jpi,jpj) ::   ztio_u, ztio_v   ! ocean stress below sea-ice
+      REAL(wp), DIMENSION(jpi,jpj) ::   zqnsoce          ! save qns before its modification by ice model
 
       !!---------------------------------------------------------------------
@@ -121 +123 @@
 !!gm re-verifies the non solar expression, especially over open ocen
 !!gm
+      zqnsoce(:,:) = qns(:,:)
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -185 +188 @@
          END DO
       END DO
-  
+
+      CALL iom_put( 'qns_io_cea', qns(:,:) - zqnsoce(:,:) * pfrld(:,:) )      
+      CALL iom_put( 'qsr_io_cea', fstric(:,:) * (1. - pfrld(:,:)) )
+
       !------------------------------------------!
       !      mass flux at the ocean surface      !
@@ -275 +281 @@
 
       !-----------------------------------------------!
-      !   Storing the transmitted variables           !
+      !   Coupling variables                          !
       !-----------------------------------------------!
-
-      fr_i(:,:) = 1.0 - frld(:,:)       ! sea-ice fraction
 
       IF ( lk_cpl ) THEN           
@@ -286 +290 @@
          CALL albedo_ice( tn_ice, reshape( hicif, (/jpi,jpj,1/) ), reshape( hsnif, (/jpi,jpj,1/) ), zalbp, zalb )
          alb_ice(:,:,1) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
+         CALL iom_put( "icealb_cea", alb_ice(:,:,1) * fr_i(:,:) )  ! ice albedo
       ENDIF