Changeset 1067 for LMDZ4/trunk/libf/phylmd/cpl_mod.F90

Timestamp:

Dec 17, 2008, 2:30:13 PM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifications lie au premier niveau du modele pour la diffusion turbulent

du vent.

• Preparation pour un couplage des courrant oceaniques.

JG

File:

• LMDZ4/trunk/libf/phylmd/cpl_mod.F90 (modified) (12 diffs)

LMDZ4/trunk/libf/phylmd/cpl_mod.F90

@@ -1 @@
-!
-! $Header$
 !
 MODULE cpl_mod
@@ -67 +65 @@
   REAL, ALLOCATABLE, DIMENSION(:,:), SAVE   :: read_alb_sic ! albedo at sea ice
   !$OMP THREADPRIVATE(read_alb_sic)
+  REAL, ALLOCATABLE, DIMENSION(:,:), SAVE   :: read_u0, read_v0 ! ocean surface current
+  !$OMP THREADPRIVATE(read_u0,read_v0)
   
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: unity
@@ -177 +177 @@
     sum_error = sum_error + error
     ALLOCATE(read_alb_sic(iim, jj_nb), stat = error)
+    sum_error = sum_error + error
+
+    ALLOCATE(read_u0(iim, jj_nb), stat = error)
+    sum_error = sum_error + error
+    ALLOCATE(read_v0(iim, jj_nb), stat = error)
     sum_error = sum_error + error
 
@@ -272 +277 @@
 ! are stored in this module.
     USE surface_data
+    USE phys_state_var_mod, ONLY : rlon, rlat
+    USE Write_Field
 
     INCLUDE "indicesol.h"
@@ -296 +303 @@
     REAL, DIMENSION(iim,jj_nb,jpfldo2a)     :: tab_read_flds
     REAL, DIMENSION(klon,nbsrf)             :: pctsrf_old
+    REAL, DIMENSION(klon_mpi)               :: rlon_mpi, rlat_mpi
+    REAL, DIMENSION(iim, jj_nb)             :: tmp_lon, tmp_lat
+    REAL, DIMENSION(iim, jj_nb)             :: tmp_r0
 
 !*************************************************************************************
@@ -311 +321 @@
        time_sec=(itime-1)*dtime
 #ifdef CPP_COUPLE
-    time_sec=(itime-1)*dtime
 !$OMP MASTER
     CALL fromcpl(time_sec, tab_read_flds)
@@ -342 +351 @@
 !$OMP END MASTER
 
+       IF (cpl_current) THEN
+
+! Transform the longitudes and latitudes on 2D arrays
+          CALL gather_omp(rlon,rlon_mpi)
+          CALL gather_omp(rlat,rlat_mpi)
+!$OMP MASTER
+          CALL Grid1DTo2D_mpi(rlon_mpi,tmp_lon)
+          CALL Grid1DTo2D_mpi(rlat_mpi,tmp_lat)
+
+! Transform the currents from cartesian to spheric coordinates
+! tmp_r0 should be zero
+          CALL geo2atm(iim, jj_nb, tab_read_flds(:,:,5), tab_read_flds(:,:,6), tab_read_flds(:,:,7), &
+               tmp_lon, tmp_lat, &
+               read_u0(:,:), read_v0(:,:), tmp_r0(:,:))
+!$OMP END MASTER
+          CALL WriteField('read_u0',read_u0)
+          CALL WriteField('read_v0',read_v0)
+          CALL WriteField('read_r0',tmp_r0)
+       ELSE
+          read_u0(:,:) = 0.
+          read_v0(:,:) = 0.
+       ENDIF
+
 !*************************************************************************************
 !  Transform seaice fraction (read_sic : ocean-seaice mask) into global 
@@ -368 +400 @@
 !
 
-  SUBROUTINE cpl_receive_ocean_fields(knon, knindex, tsurf_new)
+  SUBROUTINE cpl_receive_ocean_fields(knon, knindex, tsurf_new, u0_new, v0_new)
 !
 ! This routine returns the field for the ocean that has been read from the coupler
@@ -384 +416 @@
 !*************************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT)      :: tsurf_new
+    REAL, DIMENSION(klon), INTENT(OUT)      :: u0_new
+    REAL, DIMENSION(klon), INTENT(OUT)      :: v0_new
 
 ! Local variables
@@ -396 +430 @@
     CALL cpl2gath(read_sst, tsurf_new, knon, knindex)
     CALL cpl2gath(read_sic, sic_new, knon, knindex)
+    CALL cpl2gath(read_u0, u0_new, knon, knindex)
+    CALL cpl2gath(read_v0, v0_new, knon, knindex)
 
 !*************************************************************************************
@@ -1120 +1156 @@
     ENDIF
 
-! Transform the wind from local atmospheric 2D coordinates to geocentric 
-! 3D coordinates 
+! Transform the wind from spherical atmospheric 2D coordinates to geocentric
+! cartesian 3D coordinates 
 !$OMP MASTER
     CALL atm2geo (iim, jj_nb, tmp_taux, tmp_tauy, tmp_lon, tmp_lat, &
@@ -1165 +1201 @@
     time_sec=(itime-1)*dtime
 #ifdef CPP_COUPLE
-    time_sec=(itime-1)*dtime
 !$OMP MASTER
     CALL intocpl(time_sec, lafin, tab_flds(:,:,:))