Changeset 3559 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/OPA_SRC/SBC/cyclone.F90

Timestamp:

2012-11-15T11:18:04+01:00 (11 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB minor changes for Tropical Cyclones module

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/OPA_SRC/SBC/cyclone.F90 (modified) (10 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/OPA_SRC/SBC/cyclone.F90

@@ -21 +21 @@
    USE in_out_manager  ! I/O manager
    USE geo2ocean       ! tools for projection on ORCA grid
+   USE wrk_nemo        ! work arrays
    USE lib_mpp       
 
@@ -46 +47 @@
 CONTAINS
 
-   SUBROUTINE wnd_cyc( kt, pwnd_i, pwnd_j, ptmask_tc )
+   SUBROUTINE wnd_cyc( kt, pwnd_i, pwnd_j )
       !!----------------------------------------------------------------------
       !!                    ***  ROUTINE wnd_cyc  ***
@@ -58 +59 @@
       REAL(wp), INTENT(out), DIMENSION(jpi,jpj) ::   pwnd_i   ! wind speed i-components at T-point ORCA direction
       REAL(wp), INTENT(out), DIMENSION(jpi,jpj) ::   pwnd_j   ! wind speed j-components at T-point ORCA direction
-      REAL(wp), INTENT(out), DIMENSION(jpi,jpj) ::   ptmask_tc ! mask = 1 where TC are added
-      ! add Manu !
-      
+      ! 
       !!
       INTEGER  ::   ji, jj , jtc        ! loop arguments
@@ -80 +79 @@
       REAL(wp) ::   zrmw                ! mean radius of Max wind of a tropical cyclone (Willoughby 2004) [m]
       REAL(wp) ::   zwnd_r, zwnd_t      ! radial and tangential components of the wind
-      REAL(wp), DIMENSION(jpi,jpj) ::   zwnd_x, zwnd_y   ! zonal and meridional components of the wind
-
+      REAL(wp) ::   zvmax               ! timestep interpolated vmax
+      REAL(wp) ::   zrlon, zrlat        ! temporary 
+      REAL(wp), DIMENSION(:,:), POINTER ::   zwnd_x, zwnd_y   ! zonal and meridional components of the wind
       REAL(wp), DIMENSION(14,5)    ::   ztct                ! tropical cyclone track data at kt
-      REAL(wp)                     ::   zvmax               ! timestep interpolated vmax
-      REAL(wp)                     ::   zrlon, zrlat        ! temporary 
-      !!
+      !
       CHARACTER(len=100) ::  cn_dir            ! Root directory for location of files
       TYPE(FLD_N), DIMENSION(1) ::   slf_i     ! array of namelist informations on the TC position
@@ -92 +90 @@
       !!--------------------------------------------------------------------
 
+      CALL wrk_alloc( jpi,jpj, zwnd_x, zwnd_y )
 
       !                                         ! ====================== !
@@ -113 +112 @@
          ENDIF
          ALLOCATE( sf(1)%fnow(14,5,1) )
-         ALLOCATE( sf(1)%fdta(14,5,2,1) )
+         ALLOCATE( sf(1)%fdta(14,5,1,2) )
          slf_i(1) = sn_tc
          !
@@ -135 +134 @@
       zwnd_y(:,:) = 0.e0 
       
-      ptmask_tc(:,:) = 0.e0  
-      ! add Manu !
-  
       DO jtc = 1, 14
          !
@@ -202 +198 @@
                   zwnd_x(ji,jj) = zwnd_x(ji,jj) - zhemi * COS(ztheta)*zwnd_t + SIN(ztheta)*zwnd_r
                   zwnd_y(ji,jj) = zwnd_y(ji,jj) + zhemi * SIN(ztheta)*zwnd_t + COS(ztheta)*zwnd_r
-
-                  ptmask_tc(ji,jj) = 1. !MAX( 1., ptmask_tc(ji,jj) ) 
-                  ! add Manu !
-  
-
                  ENDIF
                END DO
@@ -273 +264 @@
                   zwnd_y(ji,jj) = zwnd_y(ji,jj) + zhemi * SIN(ztheta)*zwnd_t + COS(ztheta)*zwnd_r
                   
-                  ptmask_tc(ji,jj) = 1. !MAX( 1., ptmask_tc(ji,jj) ) 
-                  ! add Manu !
-  
                  ENDIF
                END DO
@@ -286 +274 @@
       CALL rot_rep ( zwnd_x, zwnd_y, 'T', 'en->j', pwnd_j ) !rotation of components on ORCA grid
 
+      CALL wrk_dealloc( jpi,jpj, zwnd_x, zwnd_y )
 
    END SUBROUTINE wnd_cyc