Changeset 1017 for codes/icosagcm/trunk/src/initial

Timestamp:

02/28/20 10:25:21 (4 years ago)

Author:

ymipsl

Message:

Add possibilty to set a white noise at restart. The amplitude of white noise is fixed by the folowing parameter that can be used in run.def :

• etat0_ps_white_noise
• etat0_theta_rhodz_white_noise
• etat0_u_white_noise
• etat0_q_white_noise

The value set the relative amplitude of the random perturbation.
ex :

etat0_ps_white_noise=1e-15 -> ps(ij)=ps(ij)*(1+1e-15*ran) where ran inside [0,1[

default perturbation is 0

YM

File:

• codes/icosagcm/trunk/src/initial/etat0.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/initial/etat0.f90

@@ -95 +95 @@
     
     REAL(rstd),POINTER :: ps(:), mass(:,:)
+    REAL    :: etat0_ps_white_noise
+    REAL    :: etat0_theta_rhodz_white_noise
+    REAL    :: etat0_u_white_noise
+    REAL    :: etat0_q_white_noise
     LOGICAL :: autoinit_mass, collocated
     INTEGER :: ind
@@ -178 +182 @@
        END DO
     END IF
+    
+    etat0_ps_white_noise=0.
+    CALL getin("etat0_ps_white_noise",etat0_ps_white_noise)
+    CALL add_white_noise(f_ps, etat0_ps_white_noise)
+
+    etat0_theta_rhodz_white_noise=0.
+    CALL getin("etat0_theta_rhodz_white_noise",etat0_theta_rhodz_white_noise)
+    CALL add_white_noise(f_theta_rhodz, etat0_theta_rhodz_white_noise)
  
+    etat0_u_white_noise=0.
+    CALL getin("etat0_u_white_noise",etat0_u_white_noise)
+    CALL add_white_noise(f_u, etat0_u_white_noise)
+
+    etat0_q_white_noise=0.
+    CALL getin("etat0_q_white_noise",etat0_q_white_noise)
+    CALL add_white_noise(f_q, etat0_q_white_noise)
+    
   END SUBROUTINE etat0
 
+
+  SUBROUTINE add_white_noise(field, factor)
+  USE icosa
+  IMPLICIT NONE
+    TYPE(t_field),POINTER :: field(:)   ! INOUT
+    REAL,INTENT(IN) :: factor    
+    
+    INTEGER,ALLOCATABLE :: seed(:)
+    REAL,ALLOCATABLE :: random2d(:)
+    REAL,ALLOCATABLE :: random3d(:,:)
+    REAL,ALLOCATABLE :: random4d(:,:,:)
+    REAL,POINTER :: field2d(:)
+    REAL,POINTER :: field3d(:,:)
+    REAL,POINTER :: field4d(:,:,:)
+    INTEGER :: ind
+    INTEGER :: m
+    
+    CALL RANDOM_SEED(SIZE=m)
+    ALLOCATE(seed(m))
+    
+    DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind) .OR. .NOT. is_omp_level_master) CYCLE
+      seed=domloc_glo_ind(ind) ! to be reproducible
+      CALL RANDOM_SEED(PUT=seed)
+      IF (field(ind)%ndim==2) THEN
+        field2d=field(ind)
+        ALLOCATE(random2d(size(field2d,1)))
+        CALL RANDOM_NUMBER(random2d)
+        field2d=field2d*(1.+random2d*factor)
+        DEALLOCATE(random2d)
+      ELSE IF (field(ind)%ndim==3) THEN
+        field3d=field(ind)
+        ALLOCATE(random3d(size(field3d,1),size(field3d,2)))
+        CALL RANDOM_NUMBER(random3d)
+        field3d=field3d*(1.+random3d*factor)
+        DEALLOCATE(random3d)
+      ELSE IF (field(ind)%ndim==4) THEN
+        field4d=field(ind)
+        ALLOCATE(random4d(size(field4d,1),size(field4d,2),size(field4d,3)))
+        CALL RANDOM_NUMBER(random4d)
+        field4d=field4d*(1.+random4d*factor)
+        DEALLOCATE(random4d)
+      ENDIF
+    ENDDO
+    
+  END SUBROUTINE
+    
+    
   SUBROUTINE etat0_collocated(f_phis,f_ps,f_mass,f_theta_rhodz,f_u, f_geopot,f_W, f_q)
     USE theta2theta_rhodz_mod