Changeset 204 for LMDZ.3.3/trunk/libf/dyn3d/offlinenc.F

Timestamp:

Apr 13, 2001, 12:44:53 PM (23 years ago)

Author:

lmdz

Message:

Debogage du guidage et de la version debranchee et abandon de la version
debranchee non-netcdf FH/MAF
LF

File:

• LMDZ.3.3/trunk/libf/dyn3d/offlinenc.F (modified) (10 diffs)

LMDZ.3.3/trunk/libf/dyn3d/offlinenc.F

@@ +1 @@
+c
+c $Header
+c
       PROGRAM offlinenc
       USE ioipsl
@@ -66 +69 @@
         integer iiinj
       integer jour0,isplit,nsplit_dyn,nsplit_phy,nsplit
-      logical debut,lectstart,rnpb
+      logical debut,lectstart,rnpb,lafin
 
       EXTERNAL inidissip,iniconst,inifilr
@@ -190 +193 @@
      .              teta,q,masse,ps,phis, time_0)
          print*,'Lecture du start'
+c  On zappe le radon et le plomb.
+         q(:,:,:,1)=q(:,:,:,3)
+         q(:,:,:,2)=q(:,:,:,4)
+c Initialisation d'un traceur ` 1 pour tester l'impact du lessivage.
+         q(:,:,:,3)=1.
       else
          day_ini=0
@@ -203 +211 @@
 C      print*,'av iniconst'
 c   lecture du jour de demarrage
-c    premiere initialisation, eventuellement bidon
-      CALL iniconst
+c    premiere cnitialisation, eventuellement bidon
+      cALL iniconst
       CALL inigeom
 C
-      print*,'ENTREE DANS redecoupenc ou lectfluxnc0
+      print*,'ENTREE DANS redecoupenc ou lectfluxnc
      s        pour irec=0'
 
@@ -219 +227 @@
      s     frac_impa,frac_nucl,phis)
       else
-         call lectfluxnc0(0,masse,pbaru,pbarv,w,teta,phi,
+         call lectfluxnc(0,masse,pbaru,pbarv,w,teta,phi,
      s     nrec,avant,airefi,pphis,
      s     t,zmfu, zmfd, zen_u, zde_u,zen_d, zde_d, coefkz,
@@ -373 +381 @@
      s     frac_impa,frac_nucl,phis)
       else
-         call lectfluxnc0(irec,masse,pbaru,pbarv,w,teta,phi,
+         call lectfluxnc(irec,masse,pbaru,pbarv,w,teta,phi,
      s     nrec,avant,airefi,pphis,
      s     t,zmfu, zmfd, zen_u, zde_u,zen_d, zde_d, coefkz,
@@ -380 +388 @@
       endif
 
+
+      print*,'TESTPHYS: ON PREND LA PUIS ',1./float(nsplit*nsplit_phy)
+     s   ,'  DES FRAC A IT=',itau
+      frac_impa(:,:)=frac_impa(:,:)**(1./float(nsplit*nsplit_phy))
+      frac_nucl(:,:)=frac_nucl(:,:)**(1./float(nsplit*nsplit_phy))
 c    ...  ouverture du fichier de stockage netcdf ...
 C
@@ -391 +404 @@
          mode=1
 
-         CALL initdynav(dynhistave_file,day_ini,anne_ini,dtav,
-     .              t_ops, t_wrt, nq,mode, histaveid)
+c        CALL initdynav(dynhistave_file,day_ini,anne_ini,dtav,
+c    .              t_ops, t_wrt, nq,mode, histaveid)
 
          pi=2.*asin(1.)
@@ -832 +845 @@
              do iii=1,nsplit_phy
 C
+      lafin=.false. ! en attendant mieux.
       print*,'dtphys avant phytrac ',dtphys
-             call phytrac(rnpb,
-     I                   ecritphy, debutphy,
+      print*,'TESTPHYS: APPEL A PHYTRAC IT=',itau
+
+             call phytrac(rnpb,ecritphy,
+c    I                   ecritphy, debutphy,
+     I                   debutphy, lafin,
      I                   nq,
      I                   ngridmx,llm,dtphys,
@@ -855 +872 @@
              itauav=(itau-1)*nsplit+isplit
 c             itauav=itau*nsplit+isplit
-             CALL writedynav(histaveid, nq,mode, itauav,vcov ,
-     ,                   ucov,teta,pk,phi,q,masse,ps,phis)
+c            CALL writedynav(histaveid, nq,mode, itauav,vcov ,
+c    ,                   ucov,teta,pk,phi,q,masse,ps,phis)
 c            qmoy(:,:,:)=qmoy(:,:,:)+q(:,:,1,:)
          do iq=1,nq