Changeset 1013 for trunk/NEMO/OFF_SRC/dtadyn.F90

Timestamp:

2008-05-30T12:15:06+02:00 (16 years ago)

Author:

cetlod

Message:

adding wind speed module variable, see ticket 172

File:

• trunk/NEMO/OFF_SRC/dtadyn.F90 (modified) (6 diffs)

trunk/NEMO/OFF_SRC/dtadyn.F90

@@ -38 +38 @@
       jptaux = 1 , & ! indice in flux for taux
       jptauy = 2 , & ! indice in flux for tauy
-      jpemp = 3  , & ! indice in flux for E-P
-      jpice = 4  , & ! indice in flux for ice concentration
-      jpqsr = 5      ! indice in flux for shortwave heat flux 
+      jpwind = 3 , & ! indice in flux for wind speed
+      jpemp = 4  , & ! indice in flux for E-P
+      jpice = 5  , & ! indice in flux for ice concentration
+      jpqsr = 6      ! indice in flux for shortwave heat flux 
 
    LOGICAL , PUBLIC :: &
@@ -561 +562 @@
        ! keep needed fluxes
        !
+                    wndm(:,:) = flx(:,:,jpwind)
                     freeze(:,:) = flx(:,:,jpice)
                     emp(:,:) = flx(:,:,jpemp)
@@ -625 +627 @@
        ! keep needed fluxes
        !
+                  wndm(:,:) = flx(:,:,jpwind)
                   freeze(:,:) = flx(:,:,jpice)
                   emp(:,:)    = flx(:,:,jpemp)
@@ -680 +683 @@
 
       REAL(wp), DIMENSION(jpi,jpj) :: &
-         zemp, zqsr, zmld, zice 
+         zemp, zqsr, zmld, zice, zwspd
 #if defined key_trcbbl_dif   ||   defined key_trcbbl_adv
       REAL(wp), DIMENSION(jpi,jpj) :: &
@@ -756 +759 @@
       CALL iom_get ( numfl_t, jpdom_data, 'soshfldo', zqsr (:,:  ), jkenr )
       CALL iom_get ( numfl_t, jpdom_data, 'soicecov', zice (:,:  ), jkenr )
+      CALL iom_get ( numfl_t, jpdom_data, 'sowindsp', zwspd(:,:  ), jkenr )
 
       ! file grid-U
@@ -887 +891 @@
       ! flux :
       !
+      flxdta(:,:,jpwind,2) = zwspd(:,:) * tmask(:,:,1)
       flxdta(:,:,jpice,2)  = MIN( 1., zice(:,:) ) * tmask(:,:,1)
       flxdta(:,:,jpemp,2)  = zemp(:,:) * tmask(:,:,1)