Changeset 3802 for LMDZ6/trunk

Timestamp:

Jan 16, 2021, 7:57:49 PM (4 years ago)

Author:

lguez

Message:

Move existing constraint on ngroup from groupeun_loc to conf_gcm
(abort as soon as possible if input is bad) and add constraint on
ngroup and jjm.

Location:

LMDZ6/trunk/libf/dyn3dmem

Files:

• conf_gcm.F90 (modified) (1 diff)
• groupe_loc.F (modified) (1 diff)

LMDZ6/trunk/libf/dyn3dmem/conf_gcm.F90

@@ -345 +345 @@
   ngroup=3
   CALL getin('ngroup',ngroup)
+  if (mod(iim, 2**ngroup) /= 0) &
+       call abort_gcm("conf_gcm", 'iim must be multiple of 2**ngroup', 1)
+  if (2**ngroup > jjm + 1) &
+       call abort_gcm("conf_gcm", '2**ngroup must be <= jjm + 1', 1)
 
   ! mode_top_bound : fields towards which sponge relaxation will be done:

LMDZ6/trunk/libf/dyn3dmem/groupe_loc.F

@@ -40 +40 @@
 c$OMP THREADPRIVATE(firstcall)
 
-      data firstcall/.true./
       integer ijb,ije,jjb,jje
       
-      if (firstcall) then
-         if(mod(iim,2**ngroup).ne.0) stop'probleme du nombre ede point'
-         firstcall=.false.
-      endif
-
 c   Champs 1D