Changeset 330 for trunk/LMDZ.MARS/libf

Timestamp:

Oct 25, 2011, 10:10:00 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: changes related to g95 support. but the model does not compile apparently. notes taken in NOTES.txt but sorry g95 with LMD_MM_MARS no longer supported at least on 64bits machines. LMDZ.MARS: minor changes related to warning while trying to compile with g95.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• dimradmars.h (modified) (1 diff)
• meso_inc/meso_inc_les.F (modified) (1 diff)
• physiq.F (modified) (3 diffs)
• thermcell_dqupdown.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/dimradmars.h

@@ -49 +49 @@
 
 !  Warning : the "nir" thermal IR bands are not ordered by wavelength:
-!      iir=1 : central 15um CO2 bands \    
+!      iir=1 : central 15um CO2 bands     
 !      iir=2 : CO2 band wings    [long1co2-long2co2] MINUS central band
 !      iir=3 : 9 um band [long1ir - long1co2]

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

@@ -1 @@
-         if (callrichsl .eq. .false.) then
+         if (.not.callrichsl) then
        
              DO ig=1,ngrid

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -89 +89 @@
 c    pt(ngrid,nlayer)      Temperature (K)
 c    pq(ngrid,nlayer,nq)   Advected fields
-c    pudyn(ngrid,nlayer)    \ 
-c    pvdyn(ngrid,nlayer)     \ Dynamical temporal derivative for the
-c    ptdyn(ngrid,nlayer)     / corresponding variables
-c    pqdyn(ngrid,nlayer,nq) /
+c    pudyn(ngrid,nlayer)    | 
+c    pvdyn(ngrid,nlayer)    | Dynamical temporal derivative for the
+c    ptdyn(ngrid,nlayer)    | corresponding variables
+c    pqdyn(ngrid,nlayer,nq) |
 c    pw(ngrid,?)           vertical velocity
 c
@@ -98 +98 @@
 c   -------
 c
-c    pdu(ngrid,nlayermx)        \
-c    pdv(ngrid,nlayermx)         \  Temporal derivative of the corresponding
-c    pdt(ngrid,nlayermx)         /  variables due to physical processes.
-c    pdq(ngrid,nlayermx,nqmx)   /
-c    pdpsrf(ngrid)             /
+c    pdu(ngrid,nlayermx)       |
+c    pdv(ngrid,nlayermx)       |  Temporal derivative of the corresponding
+c    pdt(ngrid,nlayermx)       |  variables due to physical processes.
+c    pdq(ngrid,nlayermx,nqmx)  |
+c    pdpsrf(ngrid)             |
 c    tracerdyn                 call tracer in dynamical part of GCM ?
 
@@ -1773 +1773 @@
          if (calltherm .and. (z_out .gt. 0.)) then
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar)
+     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar,L_mo)
 
          zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))

trunk/LMDZ.MARS/libf/phymars/thermcell_dqupdown.F90

@@ -212 +212 @@
                  if ((q(ig,k) .gt. 0.) .and. ( q(ig,k) .gt. 1.e-15 )) then
                    !! only print if it is the thermal scheme which makes qa<0
-                   print*,'qa<0 created by thermals!!!',charvar,ig,k,fm0(ig,k),qa(ig,k-1),entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
+                   print*,'qa<0 created by thermals!!!',charvar,ig,k, &
+     &                     fm0(ig,k),qa(ig,k-1),     &
+     &                     entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
                    print*,'---------> Cancelling qa'
                  endif