Changeset 5871

Timestamp:

06/25/21 13:59:41 (3 years ago)

Author:

ymipsl

Message:

Add possibility to add white noise in LMDZ restart

YM

File:

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

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

@@ -46 +46 @@
   REAL,             ALLOCATABLE :: ucov_glo(:,:),    q_glo(:,:), phis_glo(:)
   REAL,             ALLOCATABLE :: teta_glo(:,:)
+  REAL, ALLOCATABLE ::random_ucov(:,:),random_vcov(:,:),random_teta(:,:),random_ps(:)
+  REAL :: noise_ucov, noise_vcov, noise_ps, noise_teta
+  INTEGER :: seed_noise_ucov, seed_noise_vcov, seed_noise_ps, seed_noise_teta
+  INTEGER :: seed_size
+  INTEGER, ALLOCATABLE :: seed(:)
 !-------------------------------------------------------------------------------
   modname="dynetat0_loc"
@@ -125 +130 @@
   CALL err(NF90_GET_VAR(fID,vID,time),"get",var)
 
+  CALL RANDOM_SEED(size=seed_size)
+  ALLOCATE(seed(seed_size))
   ALLOCATE(phis_glo(ip1jmp1))
   CALL get_var1("phisinit",phis_glo)
@@ -130 +137 @@
 
   ALLOCATE(ucov_glo(ip1jmp1,llm))
+  ALLOCATE(random_ucov(ip1jmp1,llm))
   CALL get_var2("ucov",ucov_glo)
-  ucov (ijb_u:ije_u,:)=ucov_glo(ijb_u:ije_u,:);  DEALLOCATE(ucov_glo)
+  seed_noise_ucov=314159265
+  CALL getin("seed_noise_ucov",seed_noise_ucov)
+  seed=seed_noise_ucov
+  CALL RANDOM_SEED(put=seed)
+  noise_ucov = 0
+  CALL getin("noise_ucov",noise_ucov)
+  CALL RANDOM_NUMBER(random_ucov)
+  ucov_glo=ucov_glo*(1+(random_ucov*2-1)*noise_ucov)
+  ucov (ijb_u:ije_u,:)=ucov_glo(ijb_u:ije_u,:);  DEALLOCATE(ucov_glo) ; DEALLOCATE(random_ucov)
 
   ALLOCATE(vcov_glo(ip1jm,llm))
+  ALLOCATE(random_vcov(ip1jm,llm))
   CALL get_var2("vcov",vcov_glo)
-  vcov (ijb_v:ije_v,:)=vcov_glo(ijb_v:ije_v,:);  DEALLOCATE(vcov_glo)
+  seed_noise_vcov=314159265
+  CALL getin("seed_noise_ucov",seed_noise_vcov)
+  seed=seed_noise_vcov
+  CALL RANDOM_SEED(put=seed)
+  noise_vcov = 0
+  CALL getin("noise_ucov",noise_vcov)
+  CALL RANDOM_NUMBER(random_vcov)
+  vcov_glo=vcov_glo*(1+(random_vcov*2-1)*noise_vcov)
+  vcov (ijb_v:ije_v,:)=vcov_glo(ijb_v:ije_v,:);  DEALLOCATE(vcov_glo) ; DEALLOCATE(random_vcov)
 
   ALLOCATE(teta_glo(ip1jmp1,llm))
+  ALLOCATE(random_teta(ip1jmp1,llm))
+  seed_noise_teta=314159265
+  CALL getin("seed_noise_teta",seed_noise_teta)
+  seed=seed_noise_teta
+  CALL RANDOM_SEED(put=seed)
+  noise_teta = 0
+  CALL getin("noise_teta",noise_teta)
+  CALL RANDOM_NUMBER(random_teta)
   CALL get_var2("teta",teta_glo)
-  teta (ijb_u:ije_u,:)=teta_glo(ijb_u:ije_u,:);  DEALLOCATE(teta_glo)
+  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)
 
   ALLOCATE(masse_glo(ip1jmp1,llm))
@@ -146 +180 @@
   
   ALLOCATE(ps_glo(ip1jmp1))
+  ALLOCATE(random_ps(ip1jmp1))
+  seed_noise_ps=314159265
+  CALL getin("seed_noise_ps",seed_noise_ps)
+  seed=seed_noise_ps
+  CALL RANDOM_SEED(put=seed)
+  noise_ps = 0
+  CALL getin("noise_ps",noise_ps)
+  CALL RANDOM_NUMBER(random_ps)
   CALL get_var1("ps",ps_glo)
-  ps   (ijb_u:ije_u)  =   ps_glo(ijb_u:ije_u);   DEALLOCATE(ps_glo)
+  ps_glo=ps_glo*(1+(random_ps*2-1)*noise_ps)
+  ps   (ijb_u:ije_u)  =   ps_glo(ijb_u:ije_u);   DEALLOCATE(ps_glo) ; DEALLOCATE(random_ps)
 
 !--- Tracers