Changeset 14630 for utils/tools/DOMAINcfg/src/domzgr.F90

Timestamp:

2021-03-23T16:19:36+01:00 (3 years ago)

Author:

jchanut

Message:

AGFdomcfg: 1) restore use of analytical grids (i.e. case where jphgr_msh>0). That may be useful to set up test cases with AGRIF. 2) Add random topography over a flat bottom if nn_bathy = -1 (in place of a Gaussian bump). This illustrates well where the interface matching and update are done. #2638

File:

• utils/tools/DOMAINcfg/src/domzgr.F90 (modified) (5 diffs)

utils/tools/DOMAINcfg/src/domzgr.F90

@@ -562 +562 @@
       INTEGER , ALLOCATABLE, DIMENSION(:,:) ::   idta   ! global domain integer data
       REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zdta   ! global domain scalar data
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zrand 
       !!----------------------------------------------------------------------
       !
@@ -574 +575 @@
          ALLOCATE( idta(jpiglo,jpjglo), STAT=ierror )
          IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'zgr_bat: unable to allocate idta array' )
-         ALLOCATE( zdta(jpiglo,jpjglo), STAT=ierror )
+         ALLOCATE( zdta(jpiglo,jpjglo), zrand(jpiglo,jpjglo), STAT=ierror )
          IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'zgr_bat: unable to allocate zdta array' )
          !
@@ -599 +600 @@
          ELSE                                         ! bump centered in the basin
             IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) '         bathymetry field: flat basin with a bump'
+!            IF(lwp) WRITE(numout,*) '         bathymetry field: flat basin with a bump'
+            IF(lwp) WRITE(numout,*) '         bathymetry field: flat basin with random noise'
             ii_bump = jpiglo / 2                           ! i-index of the bump center
             ij_bump = jpjglo / 2                           ! j-index of the bump center
@@ -605 +607 @@
             h_bump  =  2700._wp                            ! bump height (meters)
             h_oce   = gdepw_1d(jpk)                        ! background ocean depth (meters)
-            IF(lwp) WRITE(numout,*) '            bump characteristics: '
-            IF(lwp) WRITE(numout,*) '               bump center (i,j)   = ', ii_bump, ii_bump
-            IF(lwp) WRITE(numout,*) '               bump height         = ', h_bump , ' meters'
-            IF(lwp) WRITE(numout,*) '               bump radius         = ', r_bump , ' index'
-            IF(lwp) WRITE(numout,*) '            background ocean depth = ', h_oce  , ' meters'
+!            IF(lwp) WRITE(numout,*) '            bump characteristics: '
+!            IF(lwp) WRITE(numout,*) '               bump center (i,j)   = ', ii_bump, ii_bump
+!            IF(lwp) WRITE(numout,*) '               bump height         = ', h_bump , ' meters'
+!            IF(lwp) WRITE(numout,*) '               bump radius         = ', r_bump , ' index'
+!            IF(lwp) WRITE(numout,*) '            background ocean depth = ', h_oce  , ' meters'
             !                                        
+            CALL RANDOM_SEED()
+            CALL RANDOM_NUMBER(zrand)
             DO jj = 1, jpjglo                              ! zdta :
                DO ji = 1, jpiglo
-                  zi = FLOAT( ji - ii_bump ) * ppe1_m / r_bump
-                  zj = FLOAT( jj - ij_bump ) * ppe2_m / r_bump
-                  zdta(ji,jj) = h_oce - h_bump * EXP( -( zi*zi + zj*zj ) )
+!                  zi = FLOAT( ji - ii_bump ) * ppe1_m / r_bump
+!                  zj = FLOAT( jj - ij_bump ) * ppe2_m / r_bump
+!                  zdta(ji,jj) = h_oce - h_bump * EXP( -( zi*zi + zj*zj ) )
+                  zdta(ji,jj) = h_oce + 0.1_wp *h_oce * (zrand(ji,jj)-1._wp) 
                END DO
             END DO
@@ -657 +662 @@
          misfdep(:,:)=1
          !
-         DEALLOCATE( idta, zdta )
+         DEALLOCATE( idta, zdta , zrand)
          !
          !                                            ! ================ !