Changeset 6208

Timestamp:

07/21/22 12:04:05 (3 years ago)

Author:

aclsce

Message:

Fixed bug in perturbation protocol.

File:

• CONFIG_DEVT/IPSLCM6.5_work_ENSEMBLES/modeles/LMDZ/libf/dyn3dmem/dynetat0_loc.F90 (modified) (2 diffs)

CONFIG_DEVT/IPSLCM6.5_work_ENSEMBLES/modeles/LMDZ/libf/dyn3dmem/dynetat0_loc.F90

@@ -41 +41 @@
   CHARACTER(LEN=256) :: msg, var, modname
   INTEGER, PARAMETER :: length=100
-  INTEGER :: iq, fID, vID, idecal, ierr
+  INTEGER :: iq, fID, vID, idecal, ierr, ij
   REAL    :: time, tab_cntrl(length)               !--- RUN PARAMS TABLE
   REAL,             ALLOCATABLE :: vcov_glo(:,:),masse_glo(:,:),   ps_glo(:)
@@ -172 +172 @@
   CALL RANDOM_NUMBER(random_teta)
   CALL get_var2("teta",teta_glo)
+  DO ij = iip2,ip1jm,iip1
+    random_teta(ij,:)=random_teta(ij+iim,:)
+  ENDDO
+  random_teta(1:iip1,:)=0
+  random_teta(ip1jm+1:ip1jmp1,:)=0
   teta_glo=teta_glo*(1+(random_teta*2-1)*noise_teta)
   teta (ijb_u:ije_u,:)=teta_glo(ijb_u:ije_u,:);  DEALLOCATE(teta_glo) ; DEALLOCATE(random_teta)