Changeset 2100
- Timestamp:
- Jul 22, 2014, 6:33:56 PM (10 years ago)
- Location:
- LMDZ5/trunk/libf
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d/abort_gcm.F
r1907 r2100 28 28 C ierr = severity of situation ( = 0 normal ) 29 29 30 character(len=*) modname31 integer ierr32 character(len=*) message30 character(len=*), intent(in):: modname 31 integer, intent(in):: ierr 32 character(len=*), intent(in):: message 33 33 34 34 write(lunout,*) 'in abort_gcm' … … 54 54 stop 55 55 else 56 write(lunout,*) 'Houston, we have a problem ', ierr56 write(lunout,*) 'Houston, we have a problem, ierr = ', ierr 57 57 stop 1 58 58 endif -
LMDZ5/trunk/libf/dyn3dmem/abort_gcm.F
r1907 r2100 23 23 C ierr = severity of situation ( = 0 normal ) 24 24 25 character(len=*) modname25 character(len=*), intent(in):: modname 26 26 integer ierr, ierror_mpi 27 character(len=*) message27 character(len=*), intent(in):: message 28 28 29 29 write(lunout,*) 'in abort_gcm' … … 46 46 write(lunout,*) 'Everything is cool' 47 47 else 48 write(lunout,*) 'Houston, we have a problem ', ierr48 write(lunout,*) 'Houston, we have a problem, ierr = ', ierr 49 49 #ifdef CPP_MPI 50 50 C$OMP CRITICAL (MPI_ABORT_GCM) -
LMDZ5/trunk/libf/dyn3dpar/abort_gcm.F
r1907 r2100 27 27 C ierr = severity of situation ( = 0 normal ) 28 28 29 character(len=*) modname29 character(len=*), intent(in):: modname 30 30 integer ierr, ierror_mpi 31 character(len=*) message31 character(len=*), intent(in):: message 32 32 33 33 write(lunout,*) 'in abort_gcm' … … 53 53 write(lunout,*) 'Everything is cool' 54 54 else 55 write(lunout,*) 'Houston, we have a problem ', ierr55 write(lunout,*) 'Houston, we have a problem, ierr = ', ierr 56 56 #ifdef CPP_MPI 57 57 C$OMP CRITICAL (MPI_ABORT_GCM) -
LMDZ5/trunk/libf/dyn3dpar/guide_p_mod.F90
r2025 r2100 374 374 ENDIF 375 375 376 PRINT *,'---> on rentre dans guide_main'377 376 ! CALL AllGather_Field(ucov,ip1jmp1,llm) 378 377 ! CALL AllGather_Field(vcov,ip1jm,llm) -
LMDZ5/trunk/libf/phylmd/hgardfou.F90
r1992 r2100 13 13 include "iniprint.h" 14 14 REAL t(klon, klev), tsol(klon, nbsrf) 15 CHARACTER *(*)text15 CHARACTER(len=*), intent(in):: text 16 16 CHARACTER (LEN=20) :: modname = 'hgardfou' 17 CHARACTER (LEN=80) :: abort_message18 17 19 18 INTEGER i, k, nsrf … … 129 128 END DO 130 129 131 IF (.NOT. ok) THEN 132 abort_message = 'hgardfou s arrete ' // text 133 CALL abort_gcm(modname, abort_message, 1) 134 END IF 130 IF (.NOT. ok) CALL abort_gcm(modname, text, 1) 135 131 136 RETURN137 132 END SUBROUTINE hgardfou -
LMDZ5/trunk/libf/phylmd/physiq.F90
r2098 r2100 1929 1929 ntra = 1 1930 1930 1931 !========================================================================= ============1932 !ajout pour la parametrisation des poches froides: 1933 ! calcul de t_wake et t_undi: si pas de poches froides, t_wake=t_undi=t_seri1931 !========================================================================= 1932 !ajout pour la parametrisation des poches froides: calcul de 1933 !t_wake et t_undi: si pas de poches froides, t_wake=t_undi=t_seri 1934 1934 do k=1,klev 1935 1935 do i=1,klon … … 1952 1952 enddo 1953 1953 1954 !c-- Calcul de l'energie disponible ALE (J/kg) et de la puissance disponible ALP (W/m2) 1955 !c-- pour le soulevement des particules dans le modele convectif 1954 ! Calcul de l'energie disponible ALE (J/kg) et de la puissance 1955 ! disponible ALP (W/m2) pour le soulevement des particules dans 1956 ! le modele convectif 1956 1957 ! 1957 1958 do i = 1,klon … … 1981 1982 enddo 1982 1983 endif 1983 !combinaison avec ale et alp de couche limite: constantes si pas de couplage, valeurs calculees1984 !d ans le thermique sinon1984 !combinaison avec ale et alp de couche limite: constantes si pas 1985 !de couplage, valeurs calculees dans le thermique sinon 1985 1986 if (iflag_coupl.eq.0) then 1986 1987 if (debut.and.prt_level.gt.9) & … … 2004 2005 ! enddo 2005 2006 2006 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!2007 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2007 2008 ! Modif FH 2010/04/27. Sans doute temporaire. 2008 ! Deux options pour le alp_offset : constant si >?? 0 ou proportionnel ??a2009 ! w si <02010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!2009 ! Deux options pour le alp_offset : constant si >?? 0 ou 2010 ! proportionnel ??a w si <0 2011 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2011 2012 do i = 1,klon 2012 2013 ALE(i) = max(ale_wake(i),Ale_bl(i)) … … 2026 2027 endif 2027 2028 enddo 2028 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!2029 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2029 2030 2030 2031 … … 2046 2047 2047 2048 !fin calcul ale et alp 2048 !======================================================================= ==========================2049 !======================================================================= 2049 2050 2050 2051 … … 2174 2175 ! . d_u_con, d_v_con) 2175 2176 2176 !----------------------------------------------------------------------------------------- 2177 ! ajout des tendances de la diffusion turbulente 2178 CALL add_phys_tend(d_u_con,d_v_con,d_t_con,d_q_con,dql0,dqi0,paprs,'con') 2179 !----------------------------------------------------------------------------------------- 2177 CALL add_phys_tend(d_u_con, d_v_con, d_t_con, d_q_con, dql0, dqi0, paprs, & 2178 'convection') 2179 !---------------------------------------------------------------------------- 2180 2180 2181 2181 if (mydebug) then … … 2287 2287 ,wake_ddeltat,wake_ddeltaq) 2288 2288 ! 2289 !------------------------------------------------------------------------- ----------------2289 !------------------------------------------------------------------------- 2290 2290 ! ajout des tendances des poches froides 2291 2291 ! Faire rapidement disparaitre l'ancien dt_wake pour garder un d_t_wake … … 2294 2294 d_q_wake(:,:)=dq_wake(:,:)*dtime 2295 2295 CALL add_phys_tend(du0,dv0,d_t_wake,d_q_wake,dql0,dqi0,paprs,'wake') 2296 !------------------------------------------------------------------------ -----------------2296 !------------------------------------------------------------------------ 2297 2297 2298 2298 endif … … 2551 2551 endif 2552 2552 2553 !----------------------------------------------------------------------- ------------------2553 !----------------------------------------------------------------------- 2554 2554 ! ajout des tendances de l'ajustement sec ou des thermiques 2555 2555 CALL add_phys_tend(du0,dv0,d_t_ajsb,d_q_ajsb,dql0,dqi0,paprs,'ajsb') … … 2557 2557 d_q_ajs(:,:)=d_q_ajs(:,:)+d_q_ajsb(:,:) 2558 2558 2559 !--------------------------------------------------------------------- --------------------2559 !--------------------------------------------------------------------- 2560 2560 2561 2561 endif … … 2609 2609 WHERE (rain_lsc < 0) rain_lsc = 0. 2610 2610 WHERE (snow_lsc < 0) snow_lsc = 0. 2611 !----------------------------------------------------------------------------------------- 2612 ! ajout des tendances de la diffusion turbulente 2611 2613 2612 CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs,'lsc') 2614 !--------------------------------------------------------------------------- --------------2613 !--------------------------------------------------------------------------- 2615 2614 DO k = 1, klev 2616 2615 DO i = 1, klon … … 2740 2739 #else 2741 2740 2742 abort_message='You should compile with -rrtm if running with iflag_rrtm=1' 2741 abort_message = 'You should compile with -rrtm if running ' & 2742 // 'with iflag_rrtm=1' 2743 2743 call abort_gcm(modname,abort_message,1) 2744 2744 #endif
Note: See TracChangeset
for help on using the changeset viewer.